Frontogenesis and Frontolysis of a Cold Filament Driven by the Cross-Filament Wind and Wave Fields Simulated by a Large Eddy Simulation

The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.

Shallow Convection Dataset Simulated by Three Different Large Eddy Models

Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes,yet its parameterization in numerical models remains a great challenge,partly due to the lack of high-resolution observations.This study describes a large eddy simulation(LES)dataset for four shallow convection cases that differ primarily in inversion strength,which can be used as a surrogate for real data.To reduce the uncertainty in LES modeling,three different large eddy models were used,including SAM(System for Atmospheric Modeling),WRF(Weather Research and Forecasting model),and UCLA-LES.Results show that the different models generally exhibit similar behavior for each shallow convection case,despite some differences in the details of the convective structure.In addition to grid-averaged fields,conditionally sampled variables,such as in-cloud moisture and vertical velocity,are also provided,which are indispensable for calculation of the entrainment/detrainment rate.Considering the essentiality of the entraining/detraining process in the parameterization of cumulus convection,the dataset presented in this study is potentially useful for validation and improvement of the parameterization of shallow convection.

Large Eddy Simulations of Rayleigh-Bénard Convection in a Cubic Cavity

This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.

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.

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.

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 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 a planar jet flow with nanoparticle coagulation

Coagulation and growth of nanoparticles subject to large coherent structures in a planar jet has been explored by using large eddy simulation. The particle field is obtained by employing a moment method to approximate the nanoparticle general dynamic equa- tion. An incompressible fluid containing particles of 1 nm in diameter is projected into a particle-free ambient. The results show that the coherent structures dominate the evolution of the nanoparticle number intensity, diameter and polydispersity distributions as the jet develops. In addition, the coherent structures act to increase the diffusion of particles, and the vortex rolling-up makes the particles distributing more irregularly while the vortex pairing causes particle distributions to become uniform. As the jet travels downstream, the time-averaged particle number concentration becomes lower in the jet core and higher in the outskirts, whereas the time- averaged particle mass over the entire flow field maintains unaltered, and the time-averaged particle diameter and geometric standard deviations grow and reach their maximum on the interface of the jet region and the ambient.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

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 vel

A large eddy simulation of flows around an underwater vehicle model using an immersed boundary method

A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.

Simulation of Surface Wave with Large Eddy Simulation in σ-Coordinate System

A three dimensional numerical model in the σ coordinate system is developed to study the problem of waves. Turbulence effects are modeled by a subgrid scale (SGS) model with the concept of large eddy simulation (LES). The σ coordinate transformation is introduced to map the irregular physical domain of the wavy free surface and uneven bottom onto the regular computational domain of the shape of rectangular prism. The operator splitting method, which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified Navier Stokes Equation. The model is used to simulate the propagation of solitary wave and wave passing over a submerged breakwater. Numerical results are compared with available analytical solutions and experimental data in terms of velocity profiles, free surface displacement, and energy conservation. Good agreement is obtained. The method is proved to be of high accuracy and efficiency in simulating surface wave propagation and wave structure interaction. It is suitable for the large and irregular physical domain, and requiring the non uniform grid system. The present work provides a foundation for further studies of random waves, wave structure interaction, wave discharge interaction, etc.

Large eddy simulation of vertical turbulent jets under JONSWAP waves

The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.

LARGE-EDDY SIMULATION OF STRATIFIED CHANNEL FLOW

An investigation of large-eddy simulation(LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smagorinsky eddy-viscosity model and Yoshizawa eddy-viscosity model were used to describe the unresolved subgrid scale (SGS) fluctuations respectively After some numerical testing, the latter was further simplified so that if can be used in the dynamic model closure. A LES code was developed for parallel computations by using the parallel technique, and aas run on the Dawn-1000 parallel computer. To demonstrate the viability and accuracy of the code, our results are compared with and found in good agreement with available LES results.

Large eddy simulation of turbulence in ocean surface boundary layer at Zhangzi Island offshore

This study uses a large eddy simulation （LES） model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island （39°N, 122°E） during July 2009 are used to drive the LES model. The LES results capture a clear diurnal cycle in the oceanic turbulence boundary layer. The process of the heat penetration and heat distribution characteristics are analyzed through the heat flux results from the LES and their differences between two diurnal cycles are discussed as well. Energy balance and other dynamics are investigated which show that the tide-induced shear production is the main source of the turbulence energy that balanced dissipation. Momentum flux near the surface shows better agreement with atmospheric data computed by the eddy correlation method than those computed by bulk formula.

Large Eddy Simulation for Wave Breaking in the Surf Zone

In this paper, the large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads to the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, Such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale model is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this paper has been used to study the propagation of the solitary wave in constant water depth and the shoaling of the non-breaking solitary wave on a beach. The model is based on large eddy simulation, and to track free-surface movements, the Tokyo University Modified Marker and Cell (TUMMAC) method is employed. In order to ensure the accuracy of each component of this wave mathematical model, several steps have been taken to verify calculated solutions; with either analytical solutions or experimental data. For non-breaking waves, very accurate results are obtained for a solitary wave propagating over a constant depth and on a beach. Application of the model to cnoidal wave breaking in the surf zone shows that the model results are in good agreement with analytical solution and experimental data. From the present model results, it can be seen that the turbulent eddy viscosity increases from the bottom to the water surface in surf zone. In the eddy viscosity curve, there is a turn-point obviously, dividing water depth into two parts, in the upper part, the eddy viscosity becomes very large near the wave breaking position.

Large Eddy Simulation of the Effects of Plasma Actuation Strength on Film Cooling Efficiency

In this article, numerical investigation of the effects of different plasma actuation strengths on the film cooling flow characteristics has been conducted using large eddy simulation （LES）. For this numerical research, the plasma actuator is placed downstream of the trailing edge of the film cooling hole and a phenomenological model is employed to provide the electric field generated by it, resulting in the body forces. Our results show that as the plasma actuation strength grows larger, under the downward effect of the plasma actuation, the jet trajectory near the cooling hole stays closer to the wall and the recirculation region observably reduces in size. Meanwhile, the momentum injection effect of the plasma actuation also actively alters the distributions of the velocity components downstream of the cooling hole. Consequently, the influence of the plasma actuation strength on the Reynolds stress downstream of the cooling hole is remarkable. Furthermore, the plasma actuation weakens the strength of the kidney shaped vortex and prevents the jet from lifting off the wall. Therefore, with the increase of the strength of the plasma actuation, the coolant core stays closer to the wall and tends to split into two distinct regions. So the centerline film cooling efficiency is enhanced, and it is increased by 55% at most when the plasma actuation strength is 10.

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.

Large Eddy Simulations of Mixing Time in a Stirred Tank

Large eddy simulations (LES) of mixing process in a stirred tank of 0.476m diameter with a 3-narrow blade hydrofoil CBY impeller were reported. The turbulent flow field and mixing time were calculated using LES with Sma-gorinsky-Lilly subgrid scale model. The impeller rotation was modeled using the sliding mesh technique. Better agree-ment of power demand and mixing time was obtained between the experimental and the LES prediction than that by the traditional Reynolds-averaged Navier-Stokes (RANS) approach. The curve of tracer response predicted by LES was in good agreement with the experimental. The results show that LES is a reliable tool to investigate the unsteady and quasi-periodic behavior of the turbulent flow in stirred tanks.