Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

High-Resolution Micro-Siting Technique for Large Scale Wind Farm Outside of Japan Using LES Turbulence Model

The spatial distribution of wind speed varies greatly over steep complex terrain, thus the selection of an optimal site in such terrain for wind turbine construction requires great care. We have developed a numerical model for simulating unsteady flows called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COM putational Prediction of Airflow over Complex Terrain), which is based on the LES (Large-Eddy Simulation) technique. The computational domain of RIAM-COMPACT can be varied from several meters to several kilometers, and the model is able to predict airflow over complex terrain with high accuracy. The present paper discusses the application of RIAM-COMPACT to the micro-siting of wind turbines at sites outside Japan. The results from two case studies will be presented.

Large eddy simulation of aircraft wake vortex with self-adaptive grid method

A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.

Numerical Study on Flow-induced Noise for a Steam Stop-valve Using Large Eddy Simulation

The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.

Large Eddy Simulation and Study of the Urban Boundary Layer

Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.

Large eddy simulations of a circular orifice jet with and without a cross-sectional exit plate

The effect of a cross-sectional exit plane on the downstream mixing characteristics of a circular turbulent jet is in- vestigated using large eddy simulation （LES）. The turbulent jet is issued from an orifice-type nozzle at an exit Reynolds number of 5 ×104. Both instantaneous and statistical velocity fields of the jet are provided. Results show that the rates of the mean velocity decay and jet spread are both higher in the case with the exit plate than without it. The existence of the plate is found to increase the downstream entrainment rate by about 10% on average over the axial range of 8-30de （exit diameter）. Also, the presence of the plate enables the formation of vortex rings to occur further downstream by 0.5-1 .Ode. A physical insight into the near-field jet is provided to explain the importance of the boundary conditions in the evolution of a turbulent jet. In addition, a method of using the decay of the centreline velocity and the half-width of the jet to calculate the entrainment rate is proposed.

A LARGE EDDY SIMULATION OF THE NEAR WAKE OF A CIRCULAR CYLINDER

Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors.

Large eddy simulation of turbulent premixed piloted flame using artificial thickened flame model coupled with tabulated chemistry

A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.

Application of large eddy simulation in the performance study of wave blocker

The configuration and aerodynamic performance of the inlet system are important aspects in the process of installing a gas turbine on a naval vessel. Under the requirements, large eddy simulation (LES) is used to simulate the three-dimensional fluid flow in the wave blocker of a marine inlet filter. The Smagorinsky-Lilly sub-grid model was used to model motions of small-scale structures. During numerical simulation, the SIMPLE algorithm was applied. The central-differencing spatial discretization scheme and the second order accuracy finite difference for the temporal discretization were used. Simulation gives satisfactory distribution of the vorticity fields and turbulent kinetic energy.Compared with the k-ε turbulent model, the results of LES are better for the distribution of parameters. The results of experimental study in a small-scale wind tunnel indicate that numerical calculation has higher accuracy. Therefore, the methods used are worthy of reference and introduction for the design of an inlet system.

Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method

Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.

Investigation of Wall Pressure Fluctuations in a Turbulent Boundary Layer by Large Eddy Simulation

Large eddy simulation (LES) was used to investigate the space-time field of the low Mach number, fully developed turbulent boundary layer on a smooth, rigid flat plate. The wall-pressure field simulated by LES was analyzed to obtain the pressure statistics, including the wall-pressure root-mean square, skewness and flatness factors, which show the wall pressure distribution was not Gaussian. The profile of the auto-power spectral density and the contour of the streamwise wavenumber-frequency spectral density of wall-pressure were plotted. The "convection ridge" can be observed clearly and the convection velocity can be calculated from the location of the convection peak.

Large eddy simulation of water flow over series of dunes

Large eddy simulation was used to investigate the spatial development of open channel flow over a series of dunes. The three-dimensional filtered Navier-Stokes （N-S） equations were numerically solved with the fractional-step method in sigma coordinates. The subgrid-scale turbulent stress was modeled with a dynamic coherent eddy viscosity model proposed by the authors The computed velocity profiles are in good agreement with the available experimental results. The mean velocity and the turbulent Reynolds stress affected by a series of dune-shaped structures were compared and analyzed. The variation of turbulence statistics along the flow direction affected by the wavy bottom roughness has been studied. The turbulent boundary layer in a complex geographic environment can be simulated well with the proposed large eddy simulation （LES） model.

Large eddy simulation of the rotation effect on the ocean turbulence kinetic energy budget in the surface mixed layer

A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.

Study of numerical errors in direct numerical simulation and large eddy simulation

By comparing the energy spectrum and total kinetic energy, the effects of numerical errors （which arise from aliasing and discretization errors）, subgrid-scale （SGS） models, and their interactions on direct numerical simulation （DNS） and large eddy simulation （LES） are investigated, The decaying isotropic turbulence is chosen as the test case. To simulate complex geometries, both the spectral method and Pade compact difference schemes are studied. The truncated Navier-Stokes （TNS） equation model with Pade discrete filter is adopted as the SGS model. It is found that the discretization error plays a key role in DNS. Low order difference schemes may be unsuitable. However, for LES, it is found that the SGS model can represent the effect of small scales to large scales and dump the numerical errors. Therefore, reasonable results can also be obtained with a low order discretization scheme.

Large Eddy Simulation of Stratified Turbulent Channel Flow with a Dynamic Subgrid Scale Model

In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.

Application of LES in the study of inlet flow field in marine gas turbine

The structure and aerodynamics performance of gas turbine inlet system is an important part of technology of gas turbine installed on naval vessels. The design and numerical simulations of gas turbine inlet system are conducted and reliable foundation for design and manufacture of marine gas turbine inlet system of high performance is provided. Numerical simulations and experiments of two inlet system models of gas turbine are conducted with satisfactory results and are of significance to the actual application of the inlet system.

THE EXAMINATION OF TURBULENCE MODELING WITH LES DATABASE

In this paper the database obtained from LES is used to examine the algebraicturbulence model in Demuren and Rodi’s work. The results show that the prediction ofnormal Reynolas stresses and turbulence energy by means of turbulence modeling isbetter than that of shear Reynolde stresses. The comparison shows the LES methodcan be used to examine turbulence modelling.

A hybrid LES (Large Eddy Simulation)/assumed sub-grid PDF (Probability Density Function) model for supersonic turbulent combustion

A hybrid LES (Large Eddy Simulation)/assumed sub-grid PDF (Probability Density Function) closure model has been devel-oped for supersonic turbulent combustion. Scalar transport equations for all species in a given chemical kinetic mechanism were solved, which are necessary in the supersonic combustion where the non-equilibrium chemistry is essentially involved. The clipped Gaussian PDF of temperature and multivariate ? PDF of composition were used to close the sub-grid chemical sources that appear in the conservation equations. The sub-grid variances of temperature and composition were constructed based on scale similarity approach. A semi-implicit approach based on the PDF model was proposed to tackle the resulting numerical stiffness associated with finite rate chemistry. The model was applied to simulate a supersonic, coaxial H2-air burner, where both the mean and rms (root mean square) results were compared with the experimental data. In general, good agree-ments were achieved, which indicated that the present sub-grid PDF method could work well in simulating supersonic turbu-lent combustion. Moreover, the calculation showed that the sub-grid fluctuations of temperature and major species in the combustion region were of the order of 10%-20% of their rms, while the sub-grid fluctuation of hydroxyl might be as high as 40%-50% of its rms.

Investigation on the width-to-depth ratio effect on turbulent flows in a sharp meandering channel with periodic boundaries using large eddy simulations

As one of the most common river patterns in nature,meandering river has very complex flow structures in its curved channel bends,including secondary flow structure and primary flow velocity redistributions.To date,most of the studies have been carried out on the flow structures in channel bends with unavoidable influences from inlet and outlet boundaries,while a streamwise periodic boundary can overcome this shortcoming elegantly.In this paper,large eddy simulations(LES)are employed to investigate the complex flow structures in periodically continuous sharp sine-generated bends.The influence of width-to-depth ratios and dimensionless curvature radiuses are studied.The results highlight two additional vortex structures beyond the commonly known secondary currents:The recirculation zone(RZ)and the inner bank cell(IBC).The width-to-depth ratio shows the determining effect on the recirculation zone.The size of recirculation zone is usually bigger in sine-generated-curve(SGC)channel with large width-to-depth ratios.The biggest recirculation zones appear between the zero-curvature section and the apex section.The inner bank cell only forms in SGC channels with small width-to-depth ratios and low curvature.For SGC channel with large width-to-depth ratios,only one circulation cell is observed near the inner bank.The spatial variations of turbulent features are also revealed by statistical analysis based on the LES sampling data.Results highlight remarkable effect of width-to-depth ratio and dimensionless curvature radius on the turbulent kinetic energy(TKE)and bed shear stress in SGC channels.