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 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.

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.

Water jet propulsion is widely used in various military and civilian fields due to its advantages of simple structure and high propulsion efficiency. The process of mooring involves utilizing specially designed equipment to secure a ship at a designated berth. During the process of water jet propulsion, the single propeller operates within a complex and turbulent three-dimensional flow. Hence, studying the coupling between the water jet propeller and the hull is critical to comprehending the characteristics of the device and the distribution of the flow field in detail.

Firstly, we conducted computational fluid dynamics (CFD)-based self-propulsion calculations to evaluate the interaction between the hull and the propeller. We subsequently analyzed the propeller's performance and the forces acting on the hull to understand how the presence or absence of the hull influenced the water jet propeller. Finally, we performed calculations and analysis of the cavitation characteristics of the coupling between the hull and the water jet propeller, considering different rotational speeds and water depths at the bottom of the pool.

The study demonstrated that the presence of the hull boundary layer under the hull-propeller coupling condition led to reduced uniformity of propeller inlet flow and lower efficiency of the propulsion pump. However, it also increased the bias toward low-flow conditions. Additionally, increasing the impeller speed led to a gradual increase in the cavitation volume within the water jet propeller, resulting in a gradual decrease in the propeller's performance.

This research provides the technical support required for effective design and operation of water jet propulsion systems. This paper involves studying and analyzing the performance and flow field of the coupling between the hull and propeller under mooring conditions with a specified hull model.

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.

Temperature is a critical factor in the fused filament fabrication (FFF) process, which affects the flow behavior and adhesion of the melted filament and the mechanical properties of the final object. Therefore, modeling and predicting temperature in FFF is crucial for achieving high-quality prints, repeatability, process control and failure prediction. This study aims to investigate the melt deposition and temperature profile in FFF both numerically and experimentally using different Acrylonitrile Butadiene Styrene single-strand specimens. The process parameters, including layer thickness, nozzle temperature and build platform temperature, were varied.

COMSOL Multiphysics software was used to perform numerical simulations of fluid flow and heat transfer for the printed strands. The polymer melt/air interface was tracked using the coupling of continuity equation, equation of motion and the level set equation, and the heat transfer equation was used to simulate the temperature distribution in the deposited strand.

The numerical results show that increasing the nozzle temperature or layer thickness leads to an increase in temperature at points close to the nozzle, but the bed temperature is the main determinant of the overall layer temperature in low-thickness strands. The experimental temperature profile of the deposited strand was measured using an infrared (IR) thermal imager to validate the numerical results. The comparison between simulation and observed temperature at different points showed that the numerical model accurately predicts heat transfer in the three-dimensional (3D) printing of a single-strand under different conditions. Finally, a parametric analysis was performed to investigate the effect of selected parameters on the thermal history of the printed strand.

The numerical results show that increasing the nozzle temperature or layer thickness leads to an increase in temperature at points close to the nozzle, but the bed temperature is the main determinant of the overall layer temperature in low-thickness strands. The experimental temperature profile of the deposited strand was measured using an IR thermal imager to validate the numerical results. The comparison between simulation and observed temperature at different points showed that the numerical model accurately predicts heat transfer in the 3D printing of a single-strand under different conditions. Finally, a parametric analysis was performed to investigate the effect of selected parameters on the thermal history of the printed strand.

The purpose of this study is an assessment of the existing roughness models to simulate the flow in the narrow gap between corotating and rough disks. A specific configuration of the flow through the gap, which forms a minichannel with variable cross sections and rotating walls, makes it a complex problem and, therefore, worth discussing in more detail.

Two roughness models were examined, the first one was based on the wall function modification by application of the shift in the dimensionless velocity profile, and the second one was based on the correction of turbulence parameters at the wall, proposed by Aupoix. Due to the lack of data to validate that specific case, the approach to deal with was selected after a systematic study of reported test cases. It started with a zero-pressure-gradient boundary layer in the flow over a flat plate, continued with flow through minichannels with stationary walls, and finally, focused on the flow between corotating discs, pertaining each time to smooth and rough surfaces.

The limitations of the roughness models were highlighted, which make the models not reliable in the application to minichannel flows. It concerns turbulence models, near-wall discretization and roughness approaches. Aupoix’s method to account for roughness was selected, and the influence of minichannel height, mass flow rate, fluid properties and roughness height on the velocity profile between corotating discs in both smooth and rough cases was discussed.

The originality of this study is the evaluation and validation of different methods to account for the roughness in rotating mini channels, where the protrusions can cover a substantial part of the channel. Flow behavior and performance of different turbulence models were analyzed as well.