Atmospheric turbulence, which produces chaotic motions in the planetary boundary layer, can inhibit mixing between fast-reacting species produced or released at different locations. This segregation process modifies the effective rate at which reactions occur between these species, and are not appropriately accounted for in coarse resolution models since these models assume complete mixing of tracers within each gridbox. Here we present a few examples of LES-based simulations applied to chemically reactive species in a forested area with high emissions of biogenic hydrocarbons, an urban area rich in anthropogenic emissions; and a maritime area with high emissions of reduced sulfur species.

The present paper proposes an accurate characterization of a NACA 23012 airfoil at near stall conditions at a Reynolds number $Re=3\cdot 10^5$. In light of the unavoidable limits of which experiments suffer near the stall regime both in terms of effective two--dimensionality and data portability across different research groups, the present characterization is performed through high fidelity numerical simulations. Taking advantage of the local discontinuous Galerkin (LDG) LES solver, implemented in the open source library \textit{FEMilaro}, the airfoil behavior is investigated ranging from $\alpha=5^{\circ}$ to $17^{\circ}$, with a peculiar focus for $\alpha \ge 10^{\circ}$. Accuracy of the found results is ensured both from the realistic and accurate reconstructed flow physics and by means of a proper comparison with existing experimental data. Then, a reliable numerical baseline is obtained for any future investigation involving a NACA 23012 airfoil at $Re=3 \cdot 10^5$. The future study of the actual research group will be a parallel BVI analysis, but results are completely extendable as they are to any other research.

Aim at a better understanding of the flow field around a fully appended Joubert BB2 submarine model and complementing the experimental investigations of the wake of the hydroplanes and sail, the Large eddy simulation (LES) with the dynamic Smagorinsky model is conducted. Three sets of grids with a maximum grid number of up to 228 million are designed, to perform the LES simulation for the Joubert BB2 at 10° yaw conditions, with a freestream Reynolds (Re) number based on the local freestream velocity and the hull length of ReL=2.2×107. And comparison of the wake of the cruciform appendage is made with experiments, to verify the computational accuracy and examine the influence of the spatial resolution. A satisfactory result shows more favorable agreement with experimental measurements as the improvement of spatial resolution, with the relative error of the vortice centers well within 7.2% in the most refined grid arrangement. And under this grid arrangement, the evolution characteristics of three co-rotating vortices originating from the cruciform appendage are further described in detail at straight ahead and 10° yaw conditions.

We have developed an unsteady and non-linear wind synopsis simulator called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain) in order to simulate the airflow on a microscale, i.e., a few tens of km or less. In RIAM-COMPACT, the large-eddy simulation (LES) has been adopted for turbulence modeling. LES is a technique in which the structures of relatively large eddies are directly simulated and smaller eddies are modeled using a sub-grid scale model. In the present study, we have conducted the numerical wind diagnoses for Taikoyama Wind Farm nacelle separation accident in Japan. The simulation results suggest that all six wind turbines at the Taikoyama Wind Farm are subject to significant influence from separated flow (terrain-induced turbulence) which is generated due to the topographic irregularities in the vicinity of the wind turbines. A proposal has been also made on reconstruction of the wind farm.

Ventilation in cities is crucial for the well being of their inhabitants. Therefore, local governments require air ventilation assessments (AVAs) prior to the construction of new buildings. In a standard AVA, however, only neutral stratification is considered, although diabatic and particularly unstable conditions may be observed more frequently in nature. The results presented here indicate significant changes in ventilation within most of the area of Kowloon City, Hong Kong, included in the study. A new definition for calculating ventilation was introduced, and used to compare the influence of buildings on ventilation under conditions of neutral and unstable stratification. The overall ventilation increased due to enhanced vertical mixing. In the vicinity of exposed buildings, however, ventilation was weaker for unstable stratification than for neutral stratification. The influence on ventilation by building parameters, such as the plan area index, was altered when unstable stratification was considered. Consequently, differences in stratification were shown to have marked effects on ventilation estimates, which should be taken into consideration in future AVAs.

The shear stress generated by wall turbulence is the main cause of wall friction resistance in turbulent flow through pipes. This paper investigates the impact of inserting rest periods (regions of constant Reynolds number) within the pulsating operating cycle of velocity on the fully turbulent flow at large time-averaged Reynolds numbers, using the Large Eddy Simulation (LES) method. The study aims to explore the effect of increasing rest periods within pulsations on drag resistance. The dimensionless shear stress and drag reduction rates during different time periods of rest were analyzed. Numerical simulation results indicate that the pulsating velocity operation mode does not necessarily lead to drag reduction; it may even result in increased resistance. Inserting rest periods within the pulsation cycle can achieve drag reduction effects, with the maximum drag reduction rate reaching 21.8\%. Comparisons with experimental data and Direct Numerical Simulation (DNS) from the literature validate the feasibility of using the LES method for pipe pulsating operation modes.

This research deals with the processing of wood using a circular saw, which is associated with several risks and factors that affect the wood splitting process. One of them is the wear of the teeth of the saw blade and the change in the size of the cutting forces. The aim of this study is to compare real measured values and theoretical values when processing Norway Spruce wood, using two types of saw blades SK and HSS with a diameter of 600 mm and the number of teeth 54 and 56. Experimental measurements were carried out on a special test device with the possibility of moving the tested wood into the cutting devices and a simple saw blade replacement system. The results of the contributions show significant differences, based on defined factors, which are supported by FEM analysis of cutting forces in the process of cutting wood - when changing the design of the saw blade. The results of the experiment show that the cutting conditions (cutting speed, feed speed) do not have a significant effect on the size of the cutting forces. The theoretical and real stress values for HSS discs ranged from 14 to 22 MPa. Significantly larger differences were recorded for discs with SK, where the real values of the maximum stresses were around 14 MPa and the theoretical ones around 25 MPa. When cutting spruce with a disc with SK slices, the maximum values of the the-oretical cutting forces were in the range of 76 to 94 N and the maximum values of the measured cutting forces were in the range of 40 to 44 N.

In the present study, field observation wind data from the time of the wind turbine blade damage accident on Shiratakiyama Wind Farm were analyzed in detail. In parallel, high-resolution LES turbulence simulations were performed in order to examine the model’s ability to numerically reproduce terrain-induced turbulence. The comparison of the observed and simulated time series (1 second average values) from a 10 minute period from the time of the accident led to the conclusion that the settings of the horizontal grid resolution and time increment are important to numerically reproduce the terrain-induced turbulence that caused the wind turbine blade damage accident on Shiratakiyama Wind Farm. A spectral analysis of the same set of observed and simulated data revealed that the simulated data reproduced the energy cascade of the actual terrain-induced turbulence well.

Computational fluid dynamics is used to study the impact of the support structure of a tidal turbine on performance and the downstream wake characteristics. A high-fidelity computational model of a dual rotor, contra-rotating tidal turbine in a large channel domain is presented, with turbulence modelled using large eddy simulation. Actuator lines represent the turbine blades, permitting the analysis of transient flow features and turbine diagnostics. The following four cases are considered: the flow in an unexploited, empty channel; flow in a channel containing the rotors; flow in a channel containing the support structure; and flow in a channel with both rotors and support structure. The results indicate that the support structure contributes significantly to the behaviour of the turbine and to turbulence levels downstream, even when the rotors are upstream. This implies that inclusion of the turbine structure, or some parametrisation thereof, is a prerequisite for the realistic prediction of turbine performance and reliability, particularly for array layouts where wake effects become significant.

This paper proposes a procedure for predicting the actual wind speed for flow over complex terrain with CFD. It converts a time-series of wind speed data acquired from field observations into a time-series of actual scalar wind speed by using non-dimensional wind speed parameters which are determined beforehand with the use of CFD output. The accuracy and reproducibility of the prediction procedure were examined by simulating the flow with CFD with the use of high resolution (5 m) surface elevation data for the Noma Wind Park in Kagoshima Prefecture, Japan. The errors of the predicted average monthly wind speeds relative to the observed values were less than approximately 20%.

This study explores the feasibility of using an oscillating plate downstream of a cylindrical body to produce mechanical energy from a Von Kármán vortex street. The study aims to quantify the impact of the plate length, its separation from the cylinder, and a machine damping factor on the power coefficient and the blade’s displacement to identify the optimal configuration. This preliminary assessment assumes that the plate oscillation is small enough to avoid changes in the vortex dynamics. This assumption allows the construction of a surrogate model using CFD to evaluate the effect of plate length and separation from the cylinder on the fluctuating lift forces over the plate. Later, the surrogate model, combined with varying machine damping factors, facilitates the description of the device’s dynamics through the numerical integration of an angular momentum equation. The results showed that a plate with 0.52D length, 5.548D separation from the cylinder, and a damping factor of 0.013 achieved a power coefficient of 0.147 and a perpendicular displacement of 0.266D. These results demonstrate a substantial improvement in the performance of bladeless generators.

In numerical simulations, achieving high accuracy without significantly increasing computational costs is often challenging. To address this, this paper presents an improved finite volume weighted essentially non-oscillatory (WENO) scheme tailored for applicability in computational fluid dynamics (CFD) and implemented in the flow solver NUAA-Turbo for simulating turbomachinery flows using both RANS and RANS/LES coupling. Firstly, the new WENO scheme is validated against classic numerical test cases to assess its stability and reliability in handling discontinuities, the Double problem, and Raleigh-Taylor (RT) instability issues. Compared to the original format, this enhanced finite volume WENO scheme demonstrates superior stability near discontinuities and resolves flow features with the same grid resolution more effectively. Next, for engineering applications related to turbomachinery, such as compressor and turbine characteristics, computations are performed using RANS, and the results obtained using WENO-ZQ3 and WENO-JS3 are compared. Finally, the new fifth-order WENO scheme is applied to RANS/LES coupled simulations of turbine wakes and film cooling. The results show that the enhanced finite volume WENO scheme offers improved stability and accuracy in engineering applications, allowing for high-precision calculations with fewer grid points to capture more detailed flow physics, thereby reducing computational costs in aerodynamic applications.

To conduct high-precision and high-resolution numerical simulation of complex flow structures in turbomachinery, a high-order finite volume weighted essentially non-oscillatory (WENO) scheme for the large eddy simulation (LES) is improved and embedded into the three-dimensional viscous unsteady CFD solver NUAA-Turbo. Firstly, the spectral characteristics and unsteady convergence of the improved WENO scheme are studied. Compared with the classical WENO scheme, the improved WENO scheme has better dissipation and dispersion characteristics and faster convergence speed. Then, the coefficient CW of the Wall-Adapting Local Eddy-viscosity (WALE) model is calibrated by decaying homogeneous isotropic turbulence (DHIT) test case. Finally, the scheme is applied to conduct high-precision LES calculations of turbulent boundary layer flow, supersonic compression corner flow, and low-pressure turbine cascade separation flow. The numerical simulation results show that this WENO scheme has excellent shock wave capture ability and turbulence resolution and can capture more flow field details with less mesh size, greatly reducing the computational cost and laying a foundation for large-scale engineering application of LES.

Due to the accuracy of numerical calculation of fluid flow inside a hydrocyclone can be obtained using Computational Fluid Dynamics (CFD), highly modified super computers are used to simulate the fluid flow and track particle motion inside a hydrocyclone. This paper deals with the numerical study using three multiphase models viz. Volume of fluid, Mixture and Eulerian model. The dimensions of the hydrocyclone taken into consideration for numerical analysis is same as considered by Rajamani. Validation of axial and tangential velocities at different strategically decided axial stations, RMS axial and tangential velocity profiles of the hydrocyclone is done using Reynolds Stress Model (RSM). The hydrocyclone model has been designed in Creo 3.0 using the same dimensions which later was imported to CFD for meshing. Fine hexagonal mesh numbering up to 5 lacs were constructed to obtain optimum results. Fluid flow was allowed to be developed in ANSYS FLUENT 16.2. Entire simulation took 96 hours to generate results and track particle movements inside the hydrocyclone. The particle tracking has been done using three multiphase model. The first being the volume of fluid was used for validation purposes and the comparison of the Mixture and Eulerian model are the basic focus of this research work. Conclusive results indicate that usage of different multiphase model does not result in variation in particle motion. The slight variation in grade efficiency values are hardly noticeable. The Mixture model and Eulerian model predict lower separation efficiency as compared with Volume of fluid multiphase model.

In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation field was applied for the solution of the problems of local strong wind areas in railway fields, and the mechanism of wind generation was discussed. At the same time, the affectivity of the application of computational fluid dynamic approach to railway field was discussed. The problem of local wind that occurs on the railway line in winter was taken up in this research. A computational simulation for the prediction of wind conditions by LES was implemented and it was clarified that the local strong wind area is mainly caused by separated flows originating from the small‐scale terrain positioned at its upstream (at approximately 180.0 m above sea level). Meanwhile, the effects of the size of calculation area and spatial grid resolution on the result of calculation and the effect of atmospheric stability were also discussed. It was clarified that when the air flow characteristic of the separated flow originating from the small‐scale terrain (at altitude of approximately 180.0 m) targeted in this research is reproduced at high accuracy by computational simulation of wind conditions, approximately 10.0 m of spatial resolution of computational grid in horizontal direction is required. As a result of the computational simulation of wind conditions of stably stratified flow (Fr = 1.0), lee waves were excited at the downstream of the terrain over time. As a result, the reverse‐flow region lying behind the terrain that had been observed at a neutral time was inhibited. Consequently, local strong wind area was generated at the downstream of the terrain and the strong wind area passing through the observation mast was observed. By investigating the speed increasing rate of local strong wind area induced at the time of stable stratification, it was found that the wind was approximately 1.2 times stronger than what was generated at a neutral time.

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in e.g. the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. The paper consists of three distinct parts – an introduction to and review of the field of computational rhinology, a review of published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of the paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

An approach based on the OpenFOAM library has been developed to solve a high-speed, multicomponent mixture of a reacting, compressible flow. This work presents comprehensive validation of the newly developed solver called compressibleCentralReactingFoam with different supersonic flows, including shocks, expansion waves, and turbulence-combustion interaction. The comparisons of the simulation results with experimental and computational data confirm the fidelity of this solver for problems involving multicomponent high-speed reactive flows. The gas dynamics of turbulence-chemistry interaction are modeled using a partially-stirred reactor formulation and provide promising results to better understand the complex physics involved in supersonic combustors. A time-scale analysis based on local Damköhler numbers reveals different regimes of turbulent combustion. In the core of the jet flow, the Damköhler number is relatively high indicating that the reaction time scale is smaller than the turbulent mixing time scale. This means that the combustion is controlled by turbulent mixing. In the shear layer, where the heat release rate and the scalar dissipation rate have the highest value, the flame is stabilized due to finite rate chemistry with small Damköhler numbers and a high fraction of fine structure. This solver allows three-dimensional gas dynamic simulation of high-speed multicomponent reactive flows relevant to practical combustion applications.

Mathematical modeling is the most important tool for constructing the theory of different kinds of two-phase flows. The review is devoted to the analysis of the advent of the approach of mathematical modeling of two-phase flows, where solid particles act mainly as the dispersed phase. The main problems and features of the study of gas-solid two-phase flows are excluded. The main characteristics of gas flows with solid particles are given and the classification of two-phase flows is developed on their basis. The Lagrangian and Euler approaches to modeling the motion of a dispersed phase (particles) are described. Much attention is paid to the consideration of numerical simulation methods that describe descriptions of turbulent gas flow at different hierarchical levels (RANS, LES, DNS), different levels of description of interphase interactions (one-way coupling (OWC), two-way coupling (TWC) and four-way coupling (FWC)), as well as at different levels of interface resolution (partial-point (PP) and particle-resolved (PR)). Excluded are examples of some studies carried out on the basis of the identified approaches, as well as their exclusion for mathematical modeling of various classes of gas-solid two-phase flows.