Large Eddy Simulation Model Research Articles

Large-eddy simulation of cultivated fields with Cyperus malaccensis Lam

ABSTRACT In this study, the honami phenomenon in Cyperus malaccensis Lam. was reproduced to establish a simulation technique that can accurately reproduce the waves of it and predict the wind flow over it and through it, which is related to its growth. The analysis code was constructed on the basis of a large-eddy simulation (LES). The LES model is based on the nonrotating coherent-structure Smagorinsky model (NRCSM), which is an expansion of the Smagorinsky model. In addition, this code is solved for a bundle of reaped rice plants, which are based on the vegetation movement, using the Newmark-beta method with an LES to analyze the movement of vegetation and air in the field. As a result, vortex structures and vegetation oscillations were visualized. We observed not only the streamwise vortex structures near the tips of C. malaccensis Lam. but also characteristic structures that sink into the vegetation simultaneously with vegetation shaking. Consequently, we found that the vibration of C. malaccensis Lam. is influenced by air inflow within the planting.

A small-eddy-dissipation mechanism for turbulence modeling and application to wall-bounded flows

A small-eddy-dissipation (SED) mechanism is proposed in the present study for the development of turbulence models. According to the SED mechanism, a turbulence model introduces artificial dissipation to filter out small eddies so that a lower resolution mesh can be used in the simulation. In addition, the artificial dissipation should be applied outside the energy-containing range so that the large-scale motions are not affected. A small-eddy-dissipation mixing length (SED-ML) model is developed based on the SED mechanism to calculate wall-bounded flows. A local Reynolds number ys+=|det(∇u)|1/3/s is introduced in the SED-ML model to distinguish laminar flows from turbulent flows. Therefore, in addition to fully turbulent flows, the SED-ML model can also be used to calculate weakly turbulent or laminar flows. To demonstrate the performance of the SED-ML model, turbulent channel flows with the Reynolds numbers Reτ up to 4200 are simulated. The numerical results are extensively compared with the large eddy simulation (LES) results using the classical subgrid-scale models. The numerical results show that the SED-ML model predicts the statistical results with a good accuracy, while requiring a lower mesh resolution than the classical LES models. The accuracy of the calculated statistical results can be further improved by the parameter extension. The friction coefficient f for channel flows can be extended directly from the reference solution according to the relation ∂f/∂ϕ/f=−9, where ϕ is a dissipative strength coefficient. The energy and dissipation spectra confirm that the SED-ML model introduces significant dissipation at high wavenumbers. This feature is beneficial in maintaining the high accuracy of the simulation results and reducing the computational cost. The numerical study also shows that a sharp filtering of the turbulent kinetic energy in the wavenumber space is essential for the development of the turbulence model with high performance.

Characteristics of instantaneous intersecting airflow from main nozzle and various relay nozzles in an air jet loom

Due to the low insertion speed and poor weft motion stability with single round-hole relay nozzle in the air jet loom, this study is aims to investigate the characteristics of transient intersecting airflow from a main nozzle and tri-elliptical hole relay nozzles. The turbulent intersecting airflow was simulated based on a large eddy simulations (LES). The correctness of the LES model was evaluated by comparing the time-averaged axial velocity with experimental and the Reynolds-averaged Navier–Stokes (RANS) results. The LES simulation has shown that only one potential-core arose its height about 4 mm for single round hole relay nozzle, whereas three potential-core arouse its height about 4 mm for tri-elliptical hole relay nozzles. It is possible to enhance the mean axial velocity of the intersecting airflow by up to 10% and decrease fluctuation magnitude by 25% using a new type of tri-elliptical center array relay nozzle.

Filter stabilization for the mildly compressible Euler equations with application to atmosphere dynamics simulations

We present a filter stabilization technique for the mildly compressible Euler equations that relies on a linear or nonlinear indicator function to identify the regions of the domain where artificial viscosity is needed and determine its amount. For the realization of this technique, we adopt a three step algorithm called Evolve-Filter-Relax (EFR), which at every time step evolves the solution (i.e., solves the Euler equations on a coarse mesh), then filters the computed solution, and finally performs a relaxation step to combine the filtered and non-filtered solutions. We show that the EFR algorithm is equivalent to an eddy-viscosity model in Large Eddy Simulation. Three indicator functions are considered: a constant function (leading to a linear filter), a function proportional to the norm of the velocity gradient (recovering a Smagorinsky-like model), and a function based on approximate deconvolution operators. Through well-known benchmarks for atmospheric flow, we show that the deconvolution-based filter yields stable solutions that are much less dissipative than the linear filter and the Smagorinsky-like model and we highlight the efficiency of the EFR algorithm.

A super-grid approach for LES combustion closure using the Linear Eddy Model

LES–LEM is a simulation approach for turbulent combustion in which the stochastic Linear Eddy Model (LEM) is used for sub-grid mixing and combustion closure in Large-Eddy Simulation (LES). LEM resolves, along a one-dimensional line, all spatial and temporal scales, provides on-the-fly local turbulent flame statistics, captures finite rate chemistry effects and directly incorporates turbulence-chemistry interaction. However, the approach is computationally expensive as it requires advancing an LEM-line in each LES cell. This paper introduces a novel turbulent combustion closure model for LES using LEM to address this issue. It involves coarse-graining the LES mesh to generate a coarse- level ‘super-grid’ comprised of cell-clusters. Each cell-cluster, instead of each LES cell, then contains a single LEM domain. This domain advances the combined advection–reaction–diffusion solution and also provides suitably conditioned statistics for thermochemical scalars such as species mass fractions. Local LES-filtered thermochemical states are then obtained by probability-density-function (PDF) weighted integration of binned conditionally averaged scalars, akin to standard presumed PDF approaches for reactive LES but with physics-based determination of the full thermochemical state for particular values of the conditioning variables. The proposed method is termed ‘super-grid LEM’ or ‘SG-LEM’. The paper describes LEM reaction–diffusion advancement, the LEM representation of turbulent advection, a novel splicing algorithm (a key feature of LES–LEM) formulated for the super-grid approach, a wall treatment, and a thermochemical LES closure procedure. To validate the proposed model, a pressure-based solver was developed using the OpenFOAM library and tested on a premixed ethylene flame stabilised over a backward facing step, a setup for which some DNS data is available. SG-LEM provides high resolution flame structures, temperature and mass fractions suitable for LES thermochemical closure. Additionally, it provides reaction-rate data at the coarse level, a unique feature compared to other mapping-type closure methods. Quantitative comparisons are made between the proposed model and time-averaged DNS data, focussing on velocity, temperature and species mass fraction. Results show good agreement downstream of the step. Furthermore, comparison with an equivalent Partially-Stirred Reactor (PaSR) simulation demonstrates the superior predictive capability of SG-LEM. Additionally, the paper briefly examines the sensitivity of the model to coarse-graining parameters and finally, explores computational efficiency highlighting the substantial speedup achieved when compared to the standard LES–LEM approach with potentially significant speedup relative to PaSR closure for the intensely turbulent regimes of principal interest.

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

A parameter-free LES model for anisotropic mesh adaptivity

Balancing accuracy and computational cost is a challenge in the modelling of turbulent flows. A widely used method for turbulence modelling is large–eddy simulation (LES). LES allows one to describe large scale flow features at a reasonable computational cost compared to the more accurate direct numerical simulation (DNS), making it a popular choice for engineering applications. One strategy to balance accuracy and cost with LES is through the use of mesh adaptivity, which allows the degrees of freedom in a problem to be reduced by changing spatial discretisation. However, mesh adaptivity can affect accuracy when using the standard Smagorinsky LES model with an implicit filter, considering that the parameter Cs is highly dependent on the filter width, which depends on mesh resolution. This work is aimed to develop an LES model that does not require any user–defined parameters and is suitable for mesh adaptivity with implicit filter. In this study we introduce a parameter–free LES model incorporating an anisotropic eddy–viscosity formulation combined with anisotropic mesh adaptivity. In our model, the parameter Cs in the eddy–viscosity formulation of the Smagorinsky model, is replaced by a function that evaluates the relative location of turbulence fluctuations in each element with respect to the turbulence spectrum inertial range. The anisotropic formulation of the eddy–viscosity allows for the application of an appropriate filter width in different directions, improving accuracy. Additionally, the mesh adaptivity algorithm assesses the local turbulence fluctuations via local Reynolds number and vortex identification criteria. This assessment leads to the refinement of regions with higher turbulence fluctuations down to the smallest scale limit in the inertial range in the corresponding direction, and also leads to the coarsening of regions without turbulence fluctuation up to largest scale limit in the inertial range. This method is tested using a flow past a sphere test case. The results are compared both qualitatively and quantitatively to results obtained with the standard Smagorinsky model, and demonstrate the better performance of our method with lower computational cost. This allows us to simulate large Reynolds number cases and compare our quantitative results to experimental results found in the literature, emphasising that our method produces good results at reasonable computational cost.

Degradation analysis of horizontal steam generator tube bundles through crack growth due to two-phase flow induced vibration

A correct understanding of vibration-based degradation is crucial from the standpoint of maintenance for Steam Generators (SG) as crucial mechanical equipment in nuclear power plants. This study has established a novel approach to developing a model for investigating tube bundle degradation according to crack growth caused by two-phase Flow-Induced Vibration (FIV). An important step in the approach is to calculate the two-phase flow field parameters between the SG tube bundles in various zones using the porous media model to determine the velocity and vapor volume fraction. Afterward, to determine the vibration properties of the tube bundles, the Fluid-Solid Interaction (FSI) analysis is performed in eighteen thermal-hydraulic zones. Tube bundle degradation based on crack growth using the sixteen most probable initial cracks and within each SG thermal-hydraulic zone is performed to calculate useful lifetime. Large Eddy Simulation (LES) model, Paris law, and Wiener process model are considered to model the turbulent crossflow around the tube bundles, simulation of elliptical crack growth due to the vibration characteristics, and estimation of SG tube bundles degradation, respectively. The analysis shows that the tube deforms most noticeably in the zone with the highest velocity. As a result, cracks propagate more quickly in the tube with a higher height. In all simulations based on different initial crack sizes, it was observed that zone 16 experiences the greatest deformation and, subsequently, the fastest degradation, with a velocity and vapor volume fraction of 0.5 m/s and 0.4, respectively.

Large eddy simulation model for two-way coupled particle-laden turbulent flows

Large eddy simulation model for two-way coupled particle-laden turbulent flows

LES model of flash-boiling and pressure recovery phenomena during release from large-scale pressurised liquid hydrogen storage tank

The paper presents large eddy simulations (LES) model of two-phase (liquid-gas) hydrogen behaviour during pressurised tank venting. The results of simulations are compared against published experiment on flash-boiling of liquid hydrogen (LH2) in a large-scale tank of 7.622 m height, 2.3 m internal diameter, filled with LH2 to 6.248 m (80% tank filling) at pressure 321.8 kPa gauge. The volume-of-fluid model is applied for tracking phases and Lee's model is used to simulate evaporation and condensation. The effect of various Lee's model evaporation-condensation coefficients on the process dynamics is assessed. The simulations are compared against experimental pressure dynamics monitored at the venting pipe and experimental temperature measured by thermocouples in gaseous and liquid phases during the first 15 s of venting. The LES model reasonably predicts physical phenomena during depressurisation of LH2 storage tank, including the pressure recovery phenomenon, LH2 level and temperature dynamics in the tank.

Ensemble data assimilation-based mixed subgrid-scale model for large-eddy simulations

An ensemble Kalman filter (EnKF)-based mixed model (EnKF-MM) is proposed for the subgrid-scale (SGS) closure in the large-eddy simulation (LES) of turbulence. The model coefficients are determined through the EnKF-based data assimilation technique. The direct numerical simulation (DNS) results are filtered to obtain the benchmark data for the LES. Reconstructing the correct kinetic energy spectrum of the filtered DNS (fDNS) data has been adopted as the target for the EnKF to optimize the coefficient of the functional part in the mixed model. The proposed EnKF-MM framework is subsequently tested in the LES of both the incompressible homogeneous isotropic turbulence and turbulent mixing layer. The performance of the LES is comprehensively examined through the predictions of the flow statistics including the velocity spectrum, the probability density functions (PDFs) of the SGS stress, the PDF of the strain rate, and the PDF of the SGS energy flux. The structure functions, the evolution of turbulent kinetic energy, the mean flow, the Reynolds stress profile, and the iso-surface of the Q-criterion are also examined to evaluate the spatial–temporal predictions by different SGS models. The results of the EnKF-MM framework are consistently more satisfying compared to the traditional SGS models, including the dynamic Smagorinsky model, the dynamic mixed model, and the velocity gradient model, demonstrating its great potential in the optimization of SGS models for the LES of turbulence.

Numerical modeling of scour and erosion processes around spur dike

AbstractSpur dikes are popular river training structures used for flow diversion, flood control, bank protection, and improving navigation. This paper aims to study the numerical modeling of scour depth in the vicinity of a rectangular spur dike using the renormalization group (RNG) turbulence model and sedimentary van Rijn model with nested mesh configuration using FLOW‐3D. This study is also carried out to analyze the three turbulence models such as RNG, 𝑘−𝜀, and Large‐Eddy simulation (LES), for estimating the scour depth and evaluating the temporal variation of maximum scour depth. It is observed that scour depth is overpredicted with the LES model configuration used in the present study, as compared to 𝑘−𝜀 and RNG models. Furthermore, the scour depth increases by increasing the bed load coefficient near to the spur dike. Simulated results show a good agreement within a 5% error with experimental results regarding the net change of sediment elevation and scour depth.

Numerical investigation on dynamic behavior of premixed hydrogen/air flame propagation in a closed tube

In this paper, the large eddy simulation (LES) model of OpenFOAM was used to study the flame propagation characteristics of hydrogen-air premixed flame in a closed tube with an aspect ratio of 12. The numerical simulation reproduces the characteristic flame shape evolution observed in the experiment, such as classical tulip flame, T-shaped flame, and distorted tulip deformation. The predicted flame front velocity, overpressure, and its fluctuation are all in good agreement with the measurements. The simulation demonstrates that after the classic tulip flame is formed, a series of waves cause an obvious adverse pressure gradient at the front of the flame. This adverse pressure gradient drives reverse flow which pulls the backward tip further into the combustion zone, forming an elongated tulip flame. A small amount of unburnt gas trapped in the cusp rapidly consumes to form a featured T-shaped flame. Then the T-shaped flame transforms into the second finger-shaped flame, and finally the tulip distortion occurs. The vortex and the reverse flow near the flame front generated by the interaction between the pressure field and the flow field accounts primarily for the formation of distorted tulip flame.

Dynamic Mixed Modeling in Large Eddy Simulation Using the Concept of a Subgrid Activity Sensor

Following the relative success of mixed models in the Large Eddy Simulation of complex turbulent flow configurations, an alternative formulation is suggested here which incorporates the concept of a local subgrid activity sensor. The general idea of mixed models is to combine the advantages of structural models (superior alignment properties), usually of the scale similarity type, and functional models (superior stability), usually of the eddy viscosity type, while avoiding their disadvantages. However, the key question is the mathematical realization of this combination, and the formulation in this work accounts for the local level of underresolution of the flow. The justification and evaluation of the newly proposed mixed model is based on a priori and a posteriori analysis of homogeneous isotropic turbulence and laminar–turbulent transition in the Taylor–Green vortex, respectively. The suggested model shows a robust and accurate behavior for the cases investigated. In particular, it outperforms the separate structural and functional base models as well as the simulation without an explicit subgrid-scale model.

Theoretical analysis of the spray characteristics of ternary blends of diesel, biodiesel, and long-chain alcohols inside a combustion chamber

The diesel and alcohol blends are potential alternatives to biogenic fuels to reduce biodiesel consumption, which might help reduce toxic emissions. Due to the high-quality requirements for engines, the fuel fundaments in the fluid properties and its impact on the macroscopic spray characteristics. Herein, a study proposed the experimental determination of density, viscosity, and surface tension for ternary mixtures of diesel (D) + biodiesel (B) + 1-hexanol (H) and + 1-octanol (O) at 298.15 K. These properties comply with the limits established by the EN 14214 and ASTM D975 standards. Thus, to extend the impact of these results, a theoretical analysis via CFD is carried out in ANSYS FLUENT ® by considering a macroscopic and isothermal system, which involved the independence tests of meshing and turbulence modeling. The numerical results evaluated the spray tip penetration, spray angle, spray morphology, eddy turbulence, and relation of spray volume and surface spray area at the injection pressure of 90 MPa using the Large Eddy Simulation (LES) model. The behavior of the spray tip penetration and spray angle agree with experimental values reported in the literature for diesel fuel. This condition validates our simulations for biodiesel and ternary blends. Therefore, the use of alcohols and biodiesel as new additives to diesel improves the macroscopic spray characteristics, and some ternary mixtures showed similar behavior to conventional diesel: spray tip penetration (70D5B25O, 70D25B5H, and 70D25B5O) and spray angle (70D25B5H, 70D15B15H, 80D5B15O, and 80D15B5). The mixture of 70D25B5H agrees with both spray characteristics; therefore, this mixture could be considered an alternative fuel to conventional diesel. Hence, this analysis presents the essential elements to visualize the macroscopic spray characteristics and the fluid properties. This knowledge serves as a foundation for elucidating the intricate dynamics of spray during combustion processes and the development of complementary studies aimed at optimizing engine performance.

Application of a Local Dynamic Model of Large Eddy Simulation to a Marine Propeller Wake

With recent development of computer technology, the use of large eddy simulation method to solve industrial problems is gaining acceptance. From a theoretical and applied perspective, a local dynamic k-equation subgrid-scale model is applied to study the flow over a marine propeller. This local dynamic SGS model of LES have already been used in simple flows such as classical Taylor-Green vortex flow to investigate its robustness and superior than other dynamic SGS models. In this paper it will be applied to a more complex flows, i.e., simulation of a marine propeller wake, to further evaluate its ability range. 42 intentionally selected numerical experiments were conducted. The results of this local dynamic model of LES shows some superior than the dynamic Smagorinsky model, and well captures the wake evolution mechanism of a propeller, although which actually depends on its geometry.

Large Eddy Simulation of Bubble Column Bubbly Flow Considering Subgrid‐Scale Turbulent Diffusion Effects and Bubble Oscillation

AbstractThrough Euler/Euler large eddy simulation (LES) modeling, it is demonstrated that turbulent dispersion of bubbles can effectively indicate the impact of turbulent eddies on the bubble dynamics, i.e., the bubble oscillation behavior. This finding builds on previous work using the Euler/Lagrange LES modeling approach and leads to a significant improvement in predicting bubble lateral dispersion. Spatially filtered terms were proposed for the subgrid‐scale (SGS) turbulent dispersion and added mass stress force models, with a modification made to the SGS eddy viscosity to reflect bubble turbulent dispersion and oscillations. The proposed model substantially improves the prediction of bubble volume fraction distribution, bubble and liquid phase velocity profiles, the turbulent kinetic energy spectrum, and mass transfer.

Implementation of an embedded LES model with parameter assessment for predicting surface pressure and surrounding flow of an isolated building

Aiming at reducing computational costs, an embedded large eddy simulation (ELES) model is introduced to simulate the wind-induced response of an isolated tall building immersed in a turbulent atmospheric boundary layer (ABL). The unsteady LES with high grid resolution is computed at the building-located area, while the low-resolution RANS model is applied to other areas. The novelty of the present model is using the random flow generation (RFG) technique to provide turbulent fluctuations, applying a mixed condition to tangential interfaces, and the good pressure coupling between separated domains. The simulated fluctuating surface pressure and the surrounding wind field of the building are analyzed, as well as the effects of fraction factor and boundary locations of the inner LES domain on simulations are also studied. Comparison with experimental measurements and pure LES result verifies the good performance of present ELES and further reveals the weakened effect of the near-wall lateral boundaries and the distortion of the neighboring RANS flow field with an insufficient downstream boundary distance. The comparison shows that over 75% of the computing time is saved using present ELES compared with pure LES, demonstrating ELES's good prospects in the urban wind environment simulation.

Direct numerical simulation of an unsteady wall-bounded turbulent flow configuration for the assessment of large-eddy simulation models

A new benchmark case for the evaluation of direct numerical simulation (DNS) and large-eddy simulation (LES) models and methods is presented in this study. The known Taylor–Green vortex is modified by replacing the periodic boundary conditions in one direction with a no-slip boundary. A passive scalar is added and transported from the wall into the fluid. The addition of walls allows for the study of transient-instationary flows in a simple geometry with clean boundary and initial conditions, which is a key requirement for the assessment of LES modeling strategies. The added scalar mimics heat transfer through the wall. The case features reasonable computational cost for highly-resolved LES and DNS calculations. Simulations of the wall-bounded Taylor–Green vortex are easy to setup and do not require additional modeling. The proposed modification of the case is compared to the default Taylor–Green vortex and the difference in flow-physics is discussed. A detailed convergence study with four meshes, each of them refined by a factor of 2, has been conducted. The results reveal that converged second-order statistics can be obtained up to a dimensionless time of t/t_0 = 20. Beyond that, due to the unsteady chaotic nature of the flow, some uncertainties remain. The results show that the case features challenging (near-wall) flow dynamics, which cannot be covered using the default Taylor–Green vortex and hence, justify the proposed case as a useful benchmark.

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.