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.

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.

In this paper, for the Reynolds-averaged Navier-Stokes equations, a self-closed turbulence model without any adjustable parameter is formulated and a simple model is proposed. The validation rule for self-closed turbulence model is rigorously derived from the Reynolds-averaged Navier-Stokes equation. The rule is not effected by turbulence modelling on the Reynolds stresses. As an application, two dimensional turbulent boundary layers, and plate turbulent boundary layers have been formulated.

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.

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.

The problem of simulating wakes in a stratified oceanic environment with active background turbulence is considered. Anisotropic RANS turbulence models are tested against laboratory and eddy-resolving models of the problem. An important aspect of our work is to acknowledge that the environment is not quiescent; therefore, additional sources are included in the models to provide a non-zero background turbulence. The RANS models are found to reproduce some key features from the eddy-resolving and laboratory descriptions of the problem. Tests using the freestream sources show the intuitive result that background turbulence causes more rapid wake growth and decay.

An open source computational fluid dynamics (CFD) solver has been incorporated into the WindNinja modeling framework widely used by wildland fire managers as well as researchers and practitioners in other fields, such as wind energy, wind erosion, and search and rescue. Here we describe incorporation of the CFD solver and evaluate its performance compared to the conservation of mass (COM) solver in WindNinja and previously published large-eddy simulations (LES) for three field campaigns conducted over isolated terrain obstacles of varying terrain complexity: Askervein Hill, Bolund Hill, and Big Southern Butte. We also compare the effects of two important model settings in the CFD solver and provide guidance on model sensitivity to these settings. Additionally, we investigate the computational mesh and difficulties regarding terrain representation. Two important findings from this work are: (1) the choice of discretization scheme for advection has a significantly larger effect on the simulated winds than the choice of turbulence model and (2) CFD solver predictions are significantly better than the COM solver predictions at windward and lee side observation locations, but no difference was found in predicted speed-up at ridgetop locations between the two solvers.

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.

The study of turbulent heat transfer in liquid metal flows has gained interest because of applications in several industrial fields. The common assumption of similarity between the dynamical and thermal turbulence, namely the Reynolds analogy, has been proven to be not valid for these fluids. Many methods have been proposed in order to overcome the difficulties encountered in a proper definition of the turbulent heat flux, such as global or local correlations for the turbulent Prandtl number or four parameter turbulence models. In this work we assess a four parameter logarithmic turbulence model for liquid metals based on RANS approach. Several simulation results considering fluids with Pr = 0.01 and Pr = 0.025 are reported in order to show the validity of this approach. The Kays turbulence model is also assessed and compared with integral heat transfer correlations for a wide range of Peclet numbers.

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.

The paper presents the prediction of hydrodynamic characteristics of different rudder -plane configurations behind a submarine in deep water through the Reynolds-Averaged Navier-Stokes method in Ansys Fluent Solver. The hydrodynamic force and moment acting on the Cross-plane, X-plane, and Y-plane rudder configurations are analyzed for estimating hydrodynamic coefficientsfor predicting the advantages of each configuration to control the ship. The obtained results show that X-force of the X-plane rudder configuration is larger than the force acting on the Cross-plane rudder and Y-plane rudder configurations. However, Y-force, and Z-force of the X-plane rudder configuration are significantly greater than the Y-force of the left configurations. The same tendency can be seen in the moment of the X-plane rudder about the y-axis and z-axis. The roll moment induced by the Y-plane and Cross-plane rudder configurations is significantly larger than the case of the X-plane rudder configuration.

This study revisits the Reynolds-averaged Navier--Stokes equations (RANS) and finds that the existing literature is erroneous regarding the primary unknowns and the number of independent unknowns in the RANS. The literature claims that the Reynolds stress tensor has six independent unknowns, but in fact the six unknowns can be reduced to three that are functions of the three velocity fluctuation components, because the Reynolds stress tensor is simply an integration of a second-order dyadic tensor of flow velocity fluctuations rather than a general symmetric tensor. This difficult situation is resolved by returning to the time of Reynolds in 1895 and revisiting Reynolds' averaging formulation of turbulence. The study of turbulence modeling could focus on the velocity fluctuations instead of on the Reynolds stress. An advantage of modeling the velocity fluctuations is, from both physical and experimental perspectives, that the velocity fluctuation components are observable whereas the Reynolds stress tensor is not.

Consumer products contain the chemical substances that threaten human health. The modeling methods and experimental methods have been used to estimate the inhalation exposure concentration by the consumer products. The model and measurement methods have the spatial property problem and time/cost consuming problem, respectively. For solving the problems due to the conventional methodology, this study performed the feasibility of applying CFD for evaluation of inhalation exposure by comparing the experiment results and the zero-dimensional results with CFD results. To calculate the aerosol concentration, the CFD was performed by combined the 3D Reynolds averaged Navier Stoke’s equation and discrete phased model using ANSYS FLUENT. As a result of comparing the three methodologies performed under the same simulation/experimental conditions, we found the zero-dimensional spray model shows approximately 5 times underestimated inhalation exposure concentration when compared with the CFD results and measurement results in near field. Also, the results of the measured concentration of aerosols at five locations and the CFD results at the same location were compared to show the possibility of evaluating inhalation exposure at various locations using CFD instead of experimental method. The CFD results according to measurement positions can predict rationally the measurement results with low error. In conclusion, in the field of exposure science, a guideline for exposure evaluation using CFD was found that complements the shortcomings of the conventional methodology, the zero-dimensional spray model and measurement method.

This paper deals with the characteristics of the cyclone separator from the Lagrangian perspective to design important dependent variables, develops a neural network model for predicting the separation performance parameter, and compares the predictive performance between the traditional surrogate model and the neural network model. In order to design the important parameters of the cyclone separator based on the particle separation theory, the force acting until the particles are separated was calculated using the Lagrangian-based CFD methodology. As a result, it was proved that the centrifugal force and drag acting on the critical diameter having a separation efficiency of 50% were similar, and the particle separation phenomenon in the cyclone occurred from the critical diameter, and it was set as an important dependent variable. For developing a critical diameter prediction model based on machine learning and multiple regression methods, Unsteady-RANS analyzes according to shape dimensions were performed. The input design variables for predicting the critical diameter were selected as four geometry parameters that affect the turbulent flow inside the cyclone. As a result of comparing the model prediction performances, the ML model showed the 32.5 % of improvement rate of R2 compared to the traditional MLR considering the nonlinear relationship between the cyclone design variable and the critical diameter. The proposed techniques have proven to be fast and practical tools for cyclone design.

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.

Steady-state Reynolds-Averaged Navier-Stokes (RANS) simulations are performed for a leading-edge inflatable wing for airborne wind energy applications. Expanding on previous work where only the inflatable leading edge tube was considered, eight additional inflatable strut tubes that support the wing canopy are now included. The shape of the wing is considered to be constant. The influence of the strut tubes on the aerodynamic performance of the wing and the local flow field is assessed, considering flow configurations with and without side-slip. The simulations show that the aerodynamic performance of the wing decreases with increasing side-slip component of the inflow. On the other hand, the chordwise struts have little influence on the integral lift and drag of the wing, irrespective of the side-slip component. The overall flow characteristics are in good agreement with previous studies. In particular, it is confirmed that at a low Reynolds number of Re=10^5, a laminar separation bubble exists on the suction side of this hypothetical rigid wing shape with perfectly smooth surface. The destruction of this bubble at low angles of attack impacts negatively on the aerodynamic performance.

Based on the insight gained from single jet analysis performed earlier [1], CFD analysis on turbulent jets impinging on one another at an angle was performed. Multiple impingement angles were considered for this study to gain better understanding of the parameters affecting resultant jet growth and velocity distribution. From the study of single jet, it was concluded that the SST k-ω turbulence model was the ideal turbulence model capable of accurately predicting the flow physics of the jets exiting a fully developed pipe at low Reynolds number. Hence, for the study of impinging jets, SST k-ω turbulence model was used to study the velocity and jet growth characteristics. It became evident that the mesh alignment with the velocity vector at exit of the pipe domain plays a crucial role in the accuracy of the results. The parameter used to evaluate this condition was identified as False Diffusion and was observed to affect the TKE parameter significantly. Methods to reduce false diffusion are also discussed in this article. Based on the mesh obtained from the grid sensitivity study, jets impinging at 30, 45 and 60 degrees at Reynolds number of 7500 were numerically analyzed. It was observed that the profile of the resultant jet closely matched with the prediction of elliptical profile predicted by past researchers, [2] and [3]. Also, it was seen that higher jet growth was predicted in case of jets impinging at a higher impingement angle.

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.

RANS simulations have been broadly used to investigate turbulence in the oceans and atmosphere. Within these environments there are a multitude of tracers undergoing reactions (eg. phytoplankton growth, chemical reactions). The distribution of these reactive tracers is strongly influenced by turbulent mixing. With a 50 member ensemble of two dimensional Rayleigh-Taylor induced turbulent mixing, we show that the dynamics of a reactive tracer growing according to Fisher’s equation are poorly captured by the mean. A fluctuation-dependent sink introduced by Reynolds averaging Fisher’s equation transfers concentration from the mean to the fluctuations. We compare the dynamics of the reactive tracer with those of a passive tracer. The reaction causes the concentration of reactive tracer to increase thereby increasing Fickian diffusion and allowing the reactive tracer to diffuse into turbulent structures that the passive tracer cannot reach. A positive feedback between turbulent mixing and fluctuation growth is identified. We show that eddy viscosity and diffusivity parameterizations fail to capture the bulk trends of the system and identify a need for negative eddy diffusivities. One must therefore be cautious when interpreting RANS results for reactive tracers.

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.