Gyrokinetic large eddy simulations
Morel, P.; Navarro, A. Banon; Albrecht-Marc, M.; Carati, D.; Merz, F.; Goerler, T.; Jenko, F.
2011-07-15
The large eddy simulation approach is adapted to the study of plasma microturbulence in a fully three-dimensional gyrokinetic system. Ion temperature gradient driven turbulence is studied with the GENE code for both a standard resolution and a reduced resolution with a model for the sub-grid scale turbulence. A simple dissipative model for representing the effect of the sub-grid scales on the resolved scales is proposed and tested. Once calibrated, the model appears to be able to reproduce most of the features of the free energy spectra for various values of the ion temperature gradient.
Large Eddy Simulation Study for Fluid Disintegration and Mixing
NASA Technical Reports Server (NTRS)
Bellan, Josette; Taskinoglu, Ezgi
2011-01-01
A new modeling approach is based on the concept of large eddy simulation (LES) within which the large scales are computed and the small scales are modeled. The new approach is expected to retain the fidelity of the physics while also being computationally efficient. Typically, only models for the small-scale fluxes of momentum, species, and enthalpy are used to reintroduce in the simulation the physics lost because the computation only resolves the large scales. These models are called subgrid (SGS) models because they operate at a scale smaller than the LES grid. In a previous study of thermodynamically supercritical fluid disintegration and mixing, additional small-scale terms, one in the momentum and one in the energy conservation equations, were identified as requiring modeling. These additional terms were due to the tight coupling between dynamics and real-gas thermodynamics. It was inferred that if these terms would not be modeled, the high density-gradient magnitude regions, experimentally identified as a characteristic feature of these flows, would not be accurately predicted without the additional term in the momentum equation; these high density-gradient magnitude regions were experimentally shown to redistribute turbulence in the flow. And it was also inferred that without the additional term in the energy equation, the heat flux magnitude could not be accurately predicted; the heat flux to the wall of combustion devices is a crucial quantity that determined necessary wall material properties. The present work involves situations where only the term in the momentum equation is important. Without this additional term in the momentum equation, neither the SGS-flux constant-coefficient Smagorinsky model nor the SGS-flux constant-coefficient Gradient model could reproduce in LES the pressure field or the high density-gradient magnitude regions; the SGS-flux constant- coefficient Scale-Similarity model was the most successful in this endeavor although not
Applied large eddy simulation.
Tucker, Paul G; Lardeau, Sylvain
2009-07-28
Large eddy simulation (LES) is now seen more and more as a viable alternative to current industrial practice, usually based on problem-specific Reynolds-averaged Navier-Stokes (RANS) methods. Access to detailed flow physics is attractive to industry, especially in an environment in which computer modelling is bound to play an ever increasing role. However, the improvement in accuracy and flow detail has substantial cost. This has so far prevented wider industrial use of LES. The purpose of the applied LES discussion meeting was to address questions regarding what is achievable and what is not, given the current technology and knowledge, for an industrial practitioner who is interested in using LES. The use of LES was explored in an application-centred context between diverse fields. The general flow-governing equation form was explored along with various LES models. The errors occurring in LES were analysed. Also, the hybridization of RANS and LES was considered. The importance of modelling relative to boundary conditions, problem definition and other more mundane aspects were examined. It was to an extent concluded that for LES to make most rapid industrial impact, pragmatic hybrid use of LES, implicit LES and RANS elements will probably be needed. Added to this further, highly industrial sector model parametrizations will be required with clear thought on the key target design parameter(s). The combination of good numerical modelling expertise, a sound understanding of turbulence, along with artistry, pragmatism and the use of recent developments in computer science should dramatically add impetus to the industrial uptake of LES. In the light of the numerous technical challenges that remain it appears that for some time to come LES will have echoes of the high levels of technical knowledge required for safe use of RANS but with much greater fidelity. PMID:19531503
Mesoscale Ocean Large Eddy Simulations
NASA Astrophysics Data System (ADS)
Pearson, Brodie; Fox-Kemper, Baylor; Bachman, Scott; Bryan, Frank
2015-11-01
The highest resolution global climate models (GCMs) can now resolve the largest scales of mesoscale dynamics in the ocean. This has the potential to increase the fidelity of GCMs. However, the effects of the smallest, unresolved, scales of mesoscale dynamics must still be parametrized. One such family of parametrizations are mesoscale ocean large eddy simulations (MOLES), but the effects of including MOLES in a GCM are not well understood. In this presentation, several MOLES schemes are implemented in a mesoscale-resolving GCM (CESM), and the resulting flow is compared with that produced by more traditional sub-grid parametrizations. Large eddy simulation (LES) is used to simulate flows where the largest scales of turbulent motion are resolved, but the smallest scales are not resolved. LES has traditionally been used to study 3D turbulence, but recently it has also been applied to idealized 2D and quasi-geostrophic (QG) turbulence. The MOLES presented here are based on 2D and QG LES schemes.
A study of differentiation errors in large-eddy simulations based on the EDQNM theory
Berland, J. Bogey, C.; Bailly, C.
2008-09-10
This paper is concerned with the investigation of numerical errors in large-eddy simulations by means of two-point turbulence modeling. Based on the eddy-damped quasi-normal Markovian (EDQNM) theory, a stochastic model is developed in order to predict the time evolution of the kinetic energy spectrum obtained by a large-eddy simulation (LES), including the effects of the numerics. Using this framework, the influence of the accuracy of the approximate space differencing schemes on LES quality is studied, for decaying homogeneous isotropic incompressible turbulence, with Reynolds numbers Re{sub {lambda}} based on the transverse Taylor scale equal to 780, 2500 and 8000. The results show that the discretization of the filtered Navier-Stokes equations leads to differentiation and aliasing errors. Error spectra are also presented, and indicate that the numerical errors are mainly originating from the approximate differentiation. In addition, increasing the order of accuracy of the differencing schemes or using algorithms optimized in the Fourier space is found to widen the range of well-resolved scales. Unfortunately, for all the schemes, the smaller scales with wavenumbers close to the grid cut-off wavenumber, are badly calculated and generate differentiation errors over the whole energy spectrum. The eventual use of explicit filtering to remove spurious motions with short wavelength is finally shown to significantly improve LES accuracy.
Temporal Large-Eddy Simulation
NASA Technical Reports Server (NTRS)
Pruett, C. D.; Thomas, B. C.
2004-01-01
In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.
A posterirori study of models for large eddy simulations of drop-laden flows
NASA Technical Reports Server (NTRS)
Leboissetier, A.; Okong'o, N. A.; Bellan, J.
2003-01-01
Large Eddy Simulation (LES) is conducted of a three-dimensional temporal mixing layer whose stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in Eulerian frame for two perfect gas species (carrier gas and vapor emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame.
Large Eddy Simulation study of fully developed thermal wind-turbine array boundary layers
NASA Astrophysics Data System (ADS)
Meneveau, Charles; Calaf, Marc; Parlange, Marc B.
2010-05-01
It is well known that when wind turbines are deployed in large arrays, their efficiency decreases due to complex interactions among themselves and with the atmospheric boundary layer (ABL). For wind farms whose length exceeds the height of the ABL by over an order of magnitude, a "fully developed" flow regime can be established. In this asymptotic regime, changes in the stream-wise direction can be neglected and the relevant exchanges occur in the vertical direction. Such a fully developed wind-turbine array boundary layer (WTABL) has recently been studied using Large Eddy Simulations (LES) under neutral stability conditions (Calaf et al. Physics of Fluids 22, 2010). Related wind-tunnel experiments on the WTABL are reported in Cal et al., J. Renewable and Sustainable Energy 2, 2010). The simulations showed the existence of two log-laws, one above and one below the wind turbine region. These results confirm basic assumptions made in prior work by Frandsen (J. Wind Eng. Ind. Aerodyn. 39, 1992) and Frandsen et al. (Wind Energy 9, 2006), and have enabled the development of more accurate parameterizations of the effective roughness scale for a wind farm. Now, a suite of Large Eddy Simulations, in which wind turbines are also modeled using the classical "drag disk" concept are performed but for non-neutral conditions. The aim is to study the effects of different thermal ABL stratifications, and thus to better understand the efficiency and characteristics of large wind farms and the associated land-atmosphere interactions for realistic atmospheric flow regimes. Such studies help to unravel the physics involved in extensive aggregations of wind turbines, allowing us to design better wind farm arrangements. By considering various turbine loading factors, surface roughness values and different atmospheric stratifications, it is possible to analyze the influence of these on the induced surface roughness, and the sensible heat roughness length. These last two can be used to
Large-eddy simulation of flow past urban-like surfaces: A model validation study
NASA Astrophysics Data System (ADS)
Cheng, Wai Chi; Porté-Agel, Fernando
2013-04-01
Accurate prediction of atmospheric boundary layer (ABL) flow and its interaction with urban surfaces is critical for understanding the transport of momentum and scalars within and above cities. This, in turn, is essential for predicting the local climate and pollutant dispersion patterns in urban areas. Large-eddy simulation (LES) explicitly resolves the large-scale turbulent eddy motions and, therefore, can potentially provide improved understanding and prediction of flows inside and above urban canopies. This study focuses on developing and validating an LES framework to simulate flow past urban-like surfaces. In particular, large-eddy simulations were performed of flow past an infinite long two-dimensional (2D) building and an array of 3D cubic buildings. An immersed boundary (IB) method was employed to simulate both 2D and 3D buildings. Four subgrid-scale (SGS) models, including (i) the traditional Smagorinsky model, (ii) the Lagrangian dynamic model, (iii) the Lagrangian scale-dependent dynamic model, and (iv) the modulated gradient model, were evaluated using the 2D building case. The simulated velocity streamlines and the vertical profiles of the mean velocities and variances were compared with experimental results. The modulated gradient model shows the best overall agreement with the experimental results among the four SGS models. In particular, the flow recirculation, the reattachment position and the vertical profiles are accurately reproduced with a grid resolution of (Nx)x(Ny)x(Nz) =160x40x160 ((nx)x(nz) =13x16 covering the block). After validating the LES framework with the 2D building case, it was further applied to simulate a boundary-layer flow past a 3D building array. A regular aligned building array with seven rows of cubic buildings was simulated. The building spacings in the streamwise and spanwise directions were both equal to the building height. A developed turbulent boundary-layer flow was used as the incoming flow. The results were
Large Eddy Simulation of a Turbulent Jet
NASA Technical Reports Server (NTRS)
Webb, A. T.; Mansour, Nagi N.
2001-01-01
Here we present the results of a Large Eddy Simulation of a non-buoyant jet issuing from a circular orifice in a wall, and developing in neutral surroundings. The effects of the subgrid scales on the large eddies have been modeled with the dynamic large eddy simulation model applied to the fully 3D domain in spherical coordinates. The simulation captures the unsteady motions of the large-scales within the jet as well as the laminar motions in the entrainment region surrounding the jet. The computed time-averaged statistics (mean velocity, concentration, and turbulence parameters) compare well with laboratory data without invoking an empirical entrainment coefficient as employed by line integral models. The use of the large eddy simulation technique allows examination of unsteady and inhomogeneous features such as the evolution of eddies and the details of the entrainment process.
Underwater Oil Plume Intrusion from Deepwater Blowouts - A Large-Eddy Simulation Study
NASA Astrophysics Data System (ADS)
Yang, D.; Chen, B.; Chamecki, M.; Meneveau, C. V.
2015-12-01
The interaction of buoyancy-driven hydrocarbon plumes with the stably stratified deep-ocean environment plays a crucial role in the formation of underwater oil intrusions. As gas bubbles and oil droplets are released from an underwater oil well blowout, they induce a strong buoyancy flux that lifts entrained sea water to form an upward plume. Towards higher elevations, the stratification-induced negative buoyancy increases and eventually exceeds the gas/oil-induced buoyancy, causing the plume to decelerate and a large fraction of entrained sea water to peel off from the rising plume to form a fountain-like downward outer plume. During this peeling process, weakly buoyant particles (e.g. small oil droplets) are trapped and fall together with the detrained fluid, and then migrate horizontally at the equilibrium buoyancy depth, forming underwater oil intrusion layers. In this study, the complex plume dynamics and oil intrusion are studied using a large-eddy simulation (LES) model. The LES model captures the essential characteristics of the plume structure and the peeling/intrusion processes, and yields good agreement with prior laboratory experiments. Applying to the Deepwater Horizon oil well blowout condition, the LES model shows considerable underwater trapping and intrusion of oil droplets under various conditions, with the trapping rate significantly affected by the diameter of the oil droplet. This study is supported by Gulf of Mexico Research Initiative RFP-II research grant.
Large eddy simulation study of spanwise spacing effects on secondary flows in turbulent channel flow
NASA Astrophysics Data System (ADS)
Aliakbarimiyanmahaleh, Mohammad; Anderson, William
2015-11-01
The structure of turbulent flow over a complex topography composed of streamwise-aligned rows of cones with varying spanwise spacing, s is studied with large-eddy simulation (LES). Similar to the experimental study of Vanderwel and Ganapathisubramani, 2015: J. Fluid Mech., we investigate the relationship between secondary flow and s, for 0 . 25 <= s / δ <= 5 . For cases with s / δ > 2 , domain-scale rollers freely exist. These had previously been called ``turbulent secondary flows'' (Willingham et al., 2014: Phys. Fluids; Barros and Christensen, 2014: J. Fluid Mech.; Anderson et al., 2015: J. Fluid Mech.), but closer inspection of the statistics indicates these are a turbulent tertiary flow: they only remain ``anchored'' to the conical roughness elements for s / δ > 2 . For s / δ < 2 , turbulent tertiary flows are prevented from occupying the domain by virtue of proximity to adjacent, counter-rotating tertiary flows. Turbulent secondary flows are associated with the conical roughness elements. These turbulent secondary flows emanate from individual conical topographic elements and set the roughness sublayer depth. The turbulent secondary flows remain intact for large and small spacing. For s / δ < 1 , a mean tertiary flow is not present. This work was supported by the Air Force Office of Sci. Research, Young Inv. Program (PM: Dr. R. Ponnoppan and Ms. E. Montomery) under Grant # FA9550-14-1-0394. Computational resources were provided by the Texas Adv. Comp. Center at the Univ. of Texas.
Influence of atmospheric stability on wind-turbine wakes: A large-eddy simulation study
NASA Astrophysics Data System (ADS)
Abkar, Mahdi; Porté-Agel, Fernando
2014-05-01
In this study, large-eddy simulation is combined with a turbine model to investigate the influence of atmospheric stability on wind-turbine wakes. In the simulations, subgrid-scale turbulent fluxes are parameterized using tuning-free Lagrangian scale-dependent dynamic models. These models optimize the local value of the model coefficients based on the dynamics of the resolved scales. The turbine-induced forces are parameterized with an actuator-disk model with rotation. In this technique, blade-element theory is used to calculate the lift and drag forces acting on the blades. Emphasis is placed on the structure and characteristics of wind-turbine wakes in the cases where the incident flows to the turbine have the same mean velocity at the hub height but different stability conditions. The simulation results show that atmospheric stability has a significant effect on the spatial distribution of the mean velocity deficit and turbulent fluxes in the wake region. In particular, the magnitude of the velocity deficit increases with increasing stability in the atmosphere. In addition, the locations of the maximum turbulence intensity and turbulent stresses are closer to the turbine in convective boundary layer compared with neutral and stable ones. Detailed analysis of the resolved turbulent kinetic energy (TKE) budget inside the wake reveals also that the thermal stratification of the incoming wind considerably affects the magnitude and spatial distribution of the turbulent production, transport term and dissipation rate (transfer of energy to the subgrid scales). It is also shown that the near-wake region can be extended to a farther distance downstream in stable condition compared with neutral and unstable counterparts. In order to isolate the effect of atmospheric stability, additional simulations of neutrally-stratified atmospheric boundary layers are performed with the same turbulence intensity at hub height as convective and stable ones. The results show that the
Amplification of turbulent exchange over wide Arctic leads: Large-eddy simulation study
NASA Astrophysics Data System (ADS)
Esau, I. N.
2007-04-01
Leads (narrow openings in the sea ice cover) are perhaps the most pronounced examples of heat islands naturally occurring on Earth. Large air-water temperature differences induce strong turbulent convection. In addition, large ice-water temperature differences induce more regular, breeze-like circulation at ice edges. Both the turbulent convection and the breeze result in intensive turbulent heat exchange between the ocean and the atmosphere. This study describes a series of turbulence-resolving experiments with the Large Eddy Simulation Nansen Center Improved Code (LESNIC). The numerical experiments quantify the turbulent heat exchange over leads of different widths. Contrary to the expected gradual decrease of the surface heat flux per unit area of open water, a strong amplification of the heat flux has been discovered for certain leads. This amplification results from a positive feedback between the horizontal entrainment of cool air in breeze and the turbulent heat exchange. Gradual reduction of the turbulent exchange for wider leads is thought to be due to development of self-organized structures in the convection. Pressure anomalies induced by the convective overturning could be comparable with the pressure anomalies due to the surface temperature difference. Their superposition limits the penetration of the cold breeze into the lead area. Without the horizontal entrainment, the near-surface temperature rises, reducing the average turbulent fluxes. In addition, the structures use the available kinetic energy to drive convective overturning. It also reduces the near-surface velocity and therefore fluxes. The maximum heat flux over open water was obtained for 2 km to 4 km leads. The maximum flux exceeds five-fold the flux in the homogeneous convection case. The revealed flux enhancement may have significant impact on the Arctic climate and more generally on the climate of urban areas and other heat islands. Therefore direct confirmation of the results from
Large-Eddy Simulation and Multigrid Methods
Falgout,R D; Naegle,S; Wittum,G
2001-06-18
A method to simulate turbulent flows with Large-Eddy Simulation on unstructured grids is presented. Two kinds of dynamic models are used to model the unresolved scales of motion and are compared with each other on different grids. Thereby the behavior of the models is shown and additionally the feature of adaptive grid refinement is investigated. Furthermore the parallelization aspect is addressed.
A Numerical Study of Self-Similarity in a Turbulent Plane Wake Using Large-Eddy Simulation
NASA Technical Reports Server (NTRS)
Ghosal, Sandip; Rogers, Michael M.
1996-01-01
Turbulent wakes are known to develop self-similarly sufficiently far downstream from obstacles that generate them. It has long been assumed that the spreading rate of the wake in the self-similar regime is independent of the details of the body generating the wake, being dependent only on the total drag (or momentum deficit). This assumption seems to be in contradiction with some recent experiments. In this study we attempt to complement these experimental investigations through a numerical study of a time-developing wake. A numerical study has the advantage of eliminating many of the uncontrolled factors present in experiments and allowing precise control of initial conditions. Large-eddy simulations employing the recently developed dynamic localization model are used to extend previous results from direct numerical simulations. The large-eddy simulation results are compared to the direct numerical simulation database, wherever such comparisons are feasible, as a check of the method. Like the experiments, the large-eddy simulations suggest that non-unique self-similar states, characterized by different spreading rates and turbulent statistics, are possible and that they can be maintained for significant time periods. The study also demonstrates the predictive capability of the dynamic localization subgrid model.
Study of tip clearance flow in a turbomachinery cascade using large eddy simulation
NASA Astrophysics Data System (ADS)
You, Donghyun
In liquid handling systems like pumps and ducted propulsors, low pressure events in the vicinity and downstream of the rotor tip gap can induce tip-leakage cavitation which leads to noise, vibration, performance loss, and erosions of blade and casing wall. In order to analyze the dynamics of the tip-clearance flow and determine the underlying mechanism for the low pressure events, a newly developed large-eddy simulation (LES) solver which combines an immersed-boundary method with a generalized curvilinear structured grid has been employed. An analysis of the LES results has been performed to understand the mean flow field, turbulence characteristics, vortex dynamics, and pressure fluctuations in the turbomachinery cascade with tip gap. In the cascade passage, the tip-leakage jet, which is generated by the pressure difference between the pressure and suction sides of the blade tip, is found to produce highly enhanced vorticity magnitude and significant levels of turbulent kinetic energy. Based on the understanding of the flow field, a guideline for reducing viscous loss in the cascade is provided. Analyses of the energy spectra and space-time correlations of the velocity fluctuations suggest that the tip-leakage vortex is subject to pitchwise wandering motion. The largest pressure drop and most intense pressure fluctuations due to the formation of the tip-leakage vortex are found at the location where the strongest portion of the tip-leakage vortex is found. Present study suggests that the tip-leakage vortex needs to be controlled in its origin to reduce cavitation in the present configuration. The effects of tip-gap size on the end-wall vortical structures and on the velocity and pressure fields have been investigated. The present analysis indicates that the mechanism for the generation of the vorticity and turbulent kinetic energy is mostly unchanged by the tip-gap size variation. However, larger tip-gap sizes are found to be more inductive to tip
Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation
NASA Technical Reports Server (NTRS)
Swanson, R. C.; Rumsey, Christopher L.; Rubinstein, Robert; Balakumar, Ponnampalam; Zang, Thomas A.
2012-01-01
Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.
Study of microphysical and radiative properties of contrail cirrus using large-eddy simulations
NASA Astrophysics Data System (ADS)
Paoli, Roberto; Thouron, Odile; Cariolle, Daniel
2013-11-01
Contrails are ice clouds that form by condensation of water vapor exhaust from aircraft engines and develop further in the wake as they are entrained by the airplane trailing vortices. When contrails spread to form contrail cirrus, they can persist for hours resulting in additional (artificial) cloud cover that adds to the cover due to natural cirrus. This talk presents recent results from large-eddy situations (LES) of contrail cirrus dispersion that are carried out using the atmospheric model Méso-NH. The objective is to investigate whether and how the ambient conditions and the microphysical and optical properties of ice crystals (e.g. shape, albedo), affect the three-dimensional structure and the overall microphysical and radiative characteristics of the contrail. The analysis is carried out by changing the radiative properties of the atmosphere (e.g. day/night conditions) for a given level of atmospheric turbulence. The turbulent field is generated by means of a stochastic forcing technique that reproduces the atmospheric conditions encountered in the upper troposphere. In addition to helping understanding the physics of contrails, the LES data retrieved from this study may provide useful inputs to the parameterization of contrail cirrus into global or climate models.
Large eddy simulation of turbulent cavitating flows
NASA Astrophysics Data System (ADS)
Gnanaskandan, A.; Mahesh, K.
2015-12-01
Large Eddy Simulation is employed to study two turbulent cavitating flows: over a cylinder and a wedge. A homogeneous mixture model is used to treat the mixture of water and water vapor as a compressible fluid. The governing equations are solved using a novel predictor- corrector method. The subgrid terms are modeled using the Dynamic Smagorinsky model. Cavitating flow over a cylinder at Reynolds number (Re) = 3900 and cavitation number (σ) = 1.0 is simulated and the wake characteristics are compared to the single phase results at the same Reynolds number. It is observed that cavitation suppresses turbulence in the near wake and delays three dimensional breakdown of the vortices. Next, cavitating flow over a wedge at Re = 200, 000 and σ = 2.0 is presented. The mean void fraction profiles obtained are compared to experiment and good agreement is obtained. Cavity auto-oscillation is observed, where the sheet cavity breaks up into a cloud cavity periodically. The results suggest LES as an attractive approach for predicting turbulent cavitating flows.
Large eddy simulation of longitudinal stationary vortices
NASA Astrophysics Data System (ADS)
Sreedhar, Madhu; Ragab, Saad
1994-07-01
The response of longitudinal stationary vortices when subjected to random perturbations is investigated using temporal large-eddy simulation. Simulations are obtained for high Reynolds numbers and at a low subsonic Mach number. The subgrid-scale stress tensor is modeled using the dynamic eddy-viscosity model. The generation of large-scale structures due to centrifugal instability and their subsequent breakdown to turbulence is studied. The following events are observed. Initially, ring-shaped structures appear around the vortex core. These structures are counter-rotating vortices similar to the donut-shaped structures observed in a Taylor-Couette flow between rotating cylinders. These structures subsequently interact with the vortex core resulting in a rapid decay of the vortex. The turbulent kinetic energy increases rapidly until saturation, and then a period of slow decay prevails. During the period of maximum turbulent kinetic energy, the normalized mean circulation profile exhibits a logarithmic region, in agreement with the universal inner profile of Hoffman and Joubert [J. Fluid Mech. 16, 395 (1963)].
Second-moment budgets in cloud topped boundary layers: A large-eddy simulation study
NASA Astrophysics Data System (ADS)
Heinze, Rieke; Mironov, Dmitrii; Raasch, Siegfried
2015-06-01
A detailed analysis of second-order moment budgets for cloud topped boundary layers (CTBLs) is performed using high-resolution large-eddy simulation (LES). Two CTBLs are simulated—one with trade wind shallow cumuli, and the other with nocturnal marine stratocumuli. Approximations to the ensemble-mean budgets of the Reynolds-stress components, of the fluxes of two quasi-conservative scalars, and of the scalar variances and covariance are computed by averaging the LES data over horizontal planes and over several hundred time steps. Importantly, the subgrid scale contributions to the budget terms are accounted for. Analysis of the LES-based second-moment budgets reveals, among other things, a paramount importance of the pressure scrambling terms in the Reynolds-stress and scalar-flux budgets. The pressure-strain correlation tends to evenly redistribute kinetic energy between the components, leading to the growth of horizontal-velocity variances at the expense of the vertical-velocity variance which is produced by buoyancy over most of both CTBLs. The pressure gradient-scalar covariances are the major sink terms in the budgets of scalar fluxes. The third-order transport proves to be of secondary importance in the scalar-flux budgets. However, it plays a key role in maintaining budgets of TKE and of the scalar variances and covariance. Results from the second-moment budget analysis suggest that the accuracy of description of the CTBL structure within the second-order closure framework strongly depends on the fidelity of parameterizations of the pressure scrambling terms in the flux budgets and of the third-order transport terms in the variance budgets. This article was corrected on 26 JUN 2015. See the end of the full text for details.
Large eddy simulation in the ocean
NASA Astrophysics Data System (ADS)
Scotti, Alberto
2010-12-01
Large eddy simulation (LES) is a relative newcomer to oceanography. In this review, both applications of traditional LES to oceanic flows and new oceanic LES still in an early stage of development are discussed. The survey covers LES applied to boundary layer flows, traditionally an area where LES has provided considerable insight into the physics of the flow, as well as more innovative applications, where new SGS closure schemes need to be developed. The merging of LES with large-scale models is also briefly reviewed.
Large eddy simulation of powered Fontan hemodynamics.
Delorme, Y; Anupindi, K; Kerlo, A E; Shetty, D; Rodefeld, M; Chen, J; Frankel, S
2013-01-18
Children born with univentricular heart disease typically must undergo three open heart surgeries within the first 2-3 years of life to eventually establish the Fontan circulation. In that case the single working ventricle pumps oxygenated blood to the body and blood returns to the lungs flowing passively through the Total Cavopulmonary Connection (TCPC) rather than being actively pumped by a subpulmonary ventricle. The TCPC is a direct surgical connection between the superior and inferior vena cava and the left and right pulmonary arteries. We have postulated that a mechanical pump inserted into this circulation providing a 3-5 mmHg pressure augmentation will reestablish bi-ventricular physiology serving as a bridge-to-recovery, bridge-to-transplant or destination therapy as a "biventricular Fontan" circulation. The Viscous Impeller Pump (VIP) has been proposed by our group as such an assist device. It is situated in the center of the 4-way TCPC intersection and spins pulling blood from the vena cavae and pushing it into the pulmonary arteries. We hypothesized that Large Eddy Simulation (LES) using high-order numerical methods are needed to capture unsteady powered and unpowered Fontan hemodynamics. Inclusion of a mechanical pump into the CFD further complicates matters due to the need to account for rotating machinery. In this study, we focus on predictions from an in-house high-order LES code (WenoHemo(TM)) for unpowered and VIP-powered idealized TCPC hemodynamics with quantitative comparisons to Stereoscopic Particle Imaging Velocimetry (SPIV) measurements. Results are presented for both instantaneous flow structures and statistical data. Simulations show good qualitative and quantitative agreement with measured data. PMID:23177085
Large Eddy Simulation of Powered Fontan Hemodynamics
Delorme, Y.; Anupindi, K.; Kerlo, A.E.; Shetty, D.; Rodefeld, M.; Chen, J.; Frankel, S.
2012-01-01
Children born with univentricular heart disease typically must undergo three open heart surgeries within the first 2–3 years of life to eventually establish the Fontan circulation. In that case the single working ventricle pumps oxygenated blood to the body and blood returns to the lungs flowing passively through the Total Cavopulmonary Connection (TCPC) rather than being actively pumped by a subpulmonary ventricle. The TCPC is a direct surgical connection between the superior and inferior vena cava and the left and right pulmonary arteries. We have postulated that a mechanical pump inserted into this circulation providing a 3–5 mmHg pressure augmentation will reestablish bi-ventricular physiology serving as a bridge-to-recovery, bridge-to-transplant or destination therapy as a “biventricular Fontan” circulation. The Viscous Impeller Pump (VIP) has been proposed by our group as such an assist device. It is situated in the center of the 4-way TCPC intersection and spins pulling blood from the vena cavae and pushing it into the pulmonary arteries. We hypothesized that Large Eddy Simulation (LES) using high-order numerical methods are needed to capture unsteady powered and unpowered Fontan hemodynamics. Inclusion of a mechanical pump into the CFD further complicates matters due to the need to account for rotating machinery. In this study, we focus on predictions from an in-house high-order LES code (WenoHemo™) for unpowered and VIP-powered idealized TCPC hemodynamics with quantitative comparisons to Stereoscopic Particle Imaging Velocimetry (SPIV) measurements. Results are presented for both instantaneous flow structures and statistical data. Simulations show good qualitative and quantitative agreement with measured data. PMID:23177085
A Large-eddy Simulation Study of Vertical Axis Wind Turbine Wakes in the Atmospheric Boundary Layer
NASA Astrophysics Data System (ADS)
Shamsoddin, Sina; Porté-Agel, Fernando
2016-04-01
Vertical axis wind turbines (VAWTs) offer some advantages over their horizontal axis counterparts, and are being considered as a viable alternative to conventional horizontal axis wind turbines (HAWTs). Nevertheless, a relative shortage of scientific, academic and technical investigations of VAWTs is observed in the wind energy community with respect to HAWTs. Having this in mind, in this work, we aim to study the wake of a single VAWT, placed in the atmospheric boundary layer, using large-eddy simulation (LES) coupled with actuator line model (ALM). It is noteworthy that this is the first time that such a study is being performed. To do this, for a typical 1 MW VAWT design, first, the variation of power coefficient with both the chord length of the blades and the tip-speed ratio is analyzed using LES-ALM, and an optimum combination of chord length and tip-speed ratio is obtained. Subsequently, the wake of a VAWT with these optimum specifications is thoroughly examined by showing different relevant mean and turbulent wake flow statistics. Keywords: vertical axis wind turbine (VAWT); VAWT wake; Atmospheric Boundary Layer (ABL); large eddy simulation (LES); actuator line model (ALM); turbulence.
Large eddy simulations in 2030 and beyond.
Piomelli, U
2014-08-13
Since its introduction, in the early 1970s, large eddy simulations (LES) have advanced considerably, and their application is transitioning from the academic environment to industry. Several landmark developments can be identified over the past 40 years, such as the wall-resolved simulations of wall-bounded flows, the development of advanced models for the unresolved scales that adapt to the local flow conditions and the hybridization of LES with the solution of the Reynolds-averaged Navier-Stokes equations. Thanks to these advancements, LES is now in widespread use in the academic community and is an option available in most commercial flow-solvers. This paper will try to predict what algorithmic and modelling advancements are needed to make it even more robust and inexpensive, and which areas show the most promise. PMID:25024415
Large eddy simulations in 2030 and beyond
Piomelli, U
2014-01-01
Since its introduction, in the early 1970s, large eddy simulations (LES) have advanced considerably, and their application is transitioning from the academic environment to industry. Several landmark developments can be identified over the past 40 years, such as the wall-resolved simulations of wall-bounded flows, the development of advanced models for the unresolved scales that adapt to the local flow conditions and the hybridization of LES with the solution of the Reynolds-averaged Navier–Stokes equations. Thanks to these advancements, LES is now in widespread use in the academic community and is an option available in most commercial flow-solvers. This paper will try to predict what algorithmic and modelling advancements are needed to make it even more robust and inexpensive, and which areas show the most promise. PMID:25024415
Autonomic Closure for Large Eddy Simulation
NASA Astrophysics Data System (ADS)
King, Ryan; Hamlington, Peter; Dahm, Werner J. A.
2015-11-01
A new autonomic subgrid-scale closure has been developed for large eddy simulation (LES). The approach poses a supervised learning problem that captures nonlinear, nonlocal, and nonequilibrium turbulence effects without specifying a predefined turbulence model. By solving a regularized optimization problem on test filter scale quantities, the autonomic approach identifies a nonparametric function that represents the best local relation between subgrid stresses and resolved state variables. The optimized function is then applied at the grid scale to determine unknown LES subgrid stresses by invoking scale similarity in the inertial range. A priori tests of the autonomic approach on homogeneous isotropic turbulence show that the new approach is amenable to powerful optimization and machine learning methods and is successful for a wide range of filter scales in the inertial range. In these a priori tests, the autonomic closure substantially improves upon the dynamic Smagorinsky model in capturing the instantaneous, statistical, and energy transfer properties of the subgrid stress field.
Large eddy simulation applications in gas turbines.
Menzies, Kevin
2009-07-28
The gas turbine presents significant challenges to any computational fluid dynamics techniques. The combination of a wide range of flow phenomena with complex geometry is difficult to model in the context of Reynolds-averaged Navier-Stokes (RANS) solvers. We review the potential for large eddy simulation (LES) in modelling the flow in the different components of the gas turbine during a practical engineering design cycle. We show that while LES has demonstrated considerable promise for reliable prediction of many flows in the engine that are difficult for RANS it is not a panacea and considerable application challenges remain. However, for many flows, especially those dominated by shear layer mixing such as in combustion chambers and exhausts, LES has demonstrated a clear superiority over RANS for moderately complex geometries although at significantly higher cost which will remain an issue in making the calculations relevant within the design cycle. PMID:19531505
Statistical Ensemble of Large Eddy Simulations
NASA Technical Reports Server (NTRS)
Carati, Daniele; Rogers, Michael M.; Wray, Alan A.; Mansour, Nagi N. (Technical Monitor)
2001-01-01
A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.
Large Eddy Simulation of Cirrus Clouds
NASA Technical Reports Server (NTRS)
Wu, Ting; Cotton, William R.
1999-01-01
The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent
Parallel Optimization with Large Eddy Simulations
NASA Astrophysics Data System (ADS)
Talnikar, Chaitanya; Blonigan, Patrick; Bodart, Julien; Wang, Qiqi; Alex Gorodetsky Collaboration; Jasper Snoek Collaboration
2014-11-01
For design optimization results to be useful, the model used must be trustworthy. For turbulent flows, Large Eddy Simulations (LES) can capture separation and other phenomena that traditional models such as RANS struggle with. However, optimization with LES can be challenging because of noisy objective function evaluations. This noise is a consequence of the sampling error of turbulent statistics, or long time averaged quantities of interest, such as the drag of an airfoil or heat transfer to a turbine blade. The sampling error causes the objective function to vary noisily with respect to design parameters for finite time simulations. Furthermore, the noise decays very slowly as computational time increases. Therefore, robustness with noisy objective functions is a crucial prerequisite to optimization candidates for LES. One way of dealing with noisy objective functions is to filter the noise using a surrogate model. Bayesian optimization, which uses Gaussian processes as surrogates, has shown promise in optimizing expensive objective functions. The following talk presents a new approach for optimization with LES incorporating these ideas. Applications to flow control of a turbulent channel and the design of a turbine blade trailing edge are also discussed.
Large eddy simulation of trailing edge noise
NASA Astrophysics Data System (ADS)
Keller, Jacob; Nitzkorski, Zane; Mahesh, Krishnan
2015-11-01
Noise generation is an important engineering constraint to many marine vehicles. A significant portion of the noise comes from propellers and rotors, specifically due to flow interactions at the trailing edge. Large eddy simulation is used to investigate the noise produced by a turbulent 45 degree beveled trailing edge and a NACA 0012 airfoil. A porous surface Ffowcs-Williams and Hawkings acoustic analogy is combined with a dynamic endcapping method to compute the sound. This methodology allows for the impact of incident flow noise versus the total noise to be assessed. LES results for the 45 degree beveled trailing edge are compared to experiment at M = 0 . 1 and Rec = 1 . 9 e 6 . The effect of boundary layer thickness on sound production is investigated by computing using both the experimental boundary layer thickness and a thinner boundary layer. Direct numerical simulation results of the NACA 0012 are compared to available data at M = 0 . 4 and Rec = 5 . 0 e 4 for both the hydrodynamic field and the acoustic field. Sound intensities and directivities are investigated and compared. Finally, some of the physical mechanisms of far-field noise generation, common to the two configurations, are discussed. Supported by Office of Naval research.
Large-Eddy Simulation of Aeroacoustic Applications
NASA Technical Reports Server (NTRS)
Pruett, C. David; Sochacki, James S.
1999-01-01
This report summarizes work accomplished under a one-year NASA grant from NASA Langley Research Center (LaRC). The effort culminates three years of NASA-supported research under three consecutive one-year grants. The period of support was April 6, 1998, through April 5, 1999. By request, the grant period was extended at no-cost until October 6, 1999. Its predecessors have been directed toward adapting the numerical tool of large-eddy simulation (LES) to aeroacoustic applications, with particular focus on noise suppression in subsonic round jets. In LES, the filtered Navier-Stokes equations are solved numerically on a relatively coarse computational grid. Residual stresses, generated by scales of motion too small to be resolved on the coarse grid, are modeled. Although most LES incorporate spatial filtering, time-domain filtering affords certain conceptual and computational advantages, particularly for aeroacoustic applications. Consequently, this work has focused on the development of subgrid-scale (SGS) models that incorporate time-domain filters.
Large Eddy Simulation of Transitional Boundary Layer
NASA Astrophysics Data System (ADS)
Sayadi, Taraneh; Moin, Parviz
2009-11-01
A sixth order compact finite difference code is employed to investigate compressible Large Eddy Simulation (LES) of subharmonic transition of a spatially developing zero pressure gradient boundary layer, at Ma = 0.2. The computational domain extends from Rex= 10^5, where laminar blowing and suction excites the most unstable fundamental and sub-harmonic modes, to fully turbulent stage at Rex= 10.1x10^5. Numerical sponges are used in the neighborhood of external boundaries to provide non-reflective conditions. Our interest lies in the performance of the dynamic subgrid scale (SGS) model [1] in the transition process. It is observed that in early stages of transition the eddy viscosity is much smaller than the physical viscosity. As a result the amplitudes of selected harmonics are in very good agreement with the experimental data [2]. The model's contribution gradually increases during the last stages of transition process and the dynamic eddy viscosity becomes fully active and dominant in the turbulent region. Consistent with this trend the skin friction coefficient versus Rex diverges from its laminar profile and converges to the turbulent profile after an overshoot. 1. Moin P. et. al. Phys Fluids A, 3(11), 2746-2757, 1991. 2. Kachanov Yu. S. et. al. JFM, 138, 209-247, 1983.
Churchfield, Matthew J; Li, Ye; Moriarty, Patrick J
2013-02-28
This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally periodic precursor simulation is performed to create turbulent flow data. Then those data are used as inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modelled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. We found that staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement. For example, using a larger precursor domain would better capture elongated turbulent structures, and including salinity and temperature equations would account for density stratification and its effect on turbulence. Additionally, the wall shear stress modelling could be improved, and more array configurations could be examined. PMID:23319713
Churchfield, M. J.; Li, Y.; Moriarty, P. J.
2012-07-01
This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally-periodic precursor simulation is performed to create turbulent flow data. Then that data is used as inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modeled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. Staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement, such as the use of a larger precursor domain to better capture elongated turbulent structures, the inclusion of salinity and temperature equations to account for density stratification and its effect on turbulence, improved wall shear stress modelling, and the examination of more array configurations.
Churchfield, M. J.; Li, Y.; Moriarty, P. J.
2011-07-01
This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally-periodic precursor simulation is performed to create turbulent flow data. Then that data is used to determine the inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modeled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. Staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement, such as the use of a larger precursor domain to better capture elongated turbulent structures, the inclusion of salinity and temperature equations to account for density stratification and its effect on turbulence, improved wall shear stress modeling, and the examination of more array configurations.
NASA Astrophysics Data System (ADS)
Cohen, Raymond; Iaccarino, Gianluca
2005-11-01
Previously published exprimental data of the flow around two circular cylinders arranged in tandem have shown that for small spacings between the cylinders, the shear layer from the upstream cylinder reattaches to the downstream cylinder, hence creating a recirculation region in between the two cylinders. The experimental data was obtained at Re=65,000 and it was found that beyond a critical spacing (L/D ˜ 4.0), the upstream shear layer ceases to attach to the downstream cylinder, resulting in a dramatic change in the flow mechanisms. Previous numerical studies using two-dimensional RANS and URANS were unsatisfactory at predicting the length of the recirculation region of the upstream cylinder and consequently badly predicted the hydrodynamic forces between the two cylinders. In this study, Large Eddy Simulation with a dynamic Smagorinsky subgrid-scale model was used to investigate the flow around two circular cylinders arranged in tandem. Results from high Reynolds numbers simulations will be presented and practical considerations in using LES in such a flow configuration will be discussed.
Large eddy simulations of laminar separation bubble
NASA Astrophysics Data System (ADS)
Cadieux, Francois
The flow over blades and airfoils at moderate angles of attack and Reynolds numbers ranging from ten thousand to a few hundred thousands undergoes separation due to the adverse pressure gradient generated by surface curvature. In many cases, the separated shear layer then transitions to turbulence and reattaches, closing off a recirculation region -- the laminar separation bubble. To avoid body-fitted mesh generation problems and numerical issues, an equivalent problem for flow over a flat plate is formulated by imposing boundary conditions that lead to a pressure distribution and Reynolds number that are similar to those on airfoils. Spalart & Strelet (2000) tested a number of Reynolds-averaged Navier-Stokes (RANS) turbulence models for a laminar separation bubble flow over a flat plate. Although results with the Spalart-Allmaras turbulence model were encouraging, none of the turbulence models tested reliably recovered time-averaged direct numerical simulation (DNS) results. The purpose of this work is to assess whether large eddy simulation (LES) can more accurately and reliably recover DNS results using drastically reduced resolution -- on the order of 1% of DNS resolution which is commonly achievable for LES of turbulent channel flows. LES of a laminar separation bubble flow over a flat plate are performed using a compressible sixth-order finite-difference code and two incompressible pseudo-spectral Navier-Stokes solvers at resolutions corresponding to approximately 3% and 1% of the chosen DNS benchmark by Spalart & Strelet (2000). The finite-difference solver is found to be dissipative due to the use of a stability-enhancing filter. Its numerical dissipation is quantified and found to be comparable to the average eddy viscosity of the dynamic Smagorinsky model, making it difficult to separate the effects of filtering versus those of explicit subgrid-scale modeling. The negligible numerical dissipation of the pseudo-spectral solvers allows an unambiguous
Wall Modeled Large Eddy Simulation of Airfoil Trailing Edge Noise
NASA Astrophysics Data System (ADS)
Kocheemoolayil, Joseph; Lele, Sanjiva
2014-11-01
Large eddy simulation (LES) of airfoil trailing edge noise has largely been restricted to low Reynolds numbers due to prohibitive computational cost. Wall modeled LES (WMLES) is a computationally cheaper alternative that makes full-scale Reynolds numbers relevant to large wind turbines accessible. A systematic investigation of trailing edge noise prediction using WMLES is conducted. Detailed comparisons are made with experimental data. The stress boundary condition from a wall model does not constrain the fluctuating velocity to vanish at the wall. This limitation has profound implications for trailing edge noise prediction. The simulation over-predicts the intensity of fluctuating wall pressure and far-field noise. An improved wall model formulation that minimizes the over-prediction of fluctuating wall pressure is proposed and carefully validated. The flow configurations chosen for the study are from the workshop on benchmark problems for airframe noise computations. The large eddy simulation database is used to examine the adequacy of scaling laws that quantify the dependence of trailing edge noise on Mach number, Reynolds number and angle of attack. Simplifying assumptions invoked in engineering approaches towards predicting trailing edge noise are critically evaluated. We gratefully acknowledge financial support from GE Global Research and thank Cascade Technologies Inc. for providing access to their massively-parallel large eddy simulation framework.
NASA Astrophysics Data System (ADS)
Zhu, Xiaowei; Anderson, William
2015-11-01
The inherent spatial heterogeneity exhibited by real urban environments complicates a priori estimation of the roughness height needed to parameterize the inertial layer mean streamwise velocity. A large-eddy simulation study of turbulent flow over 3-D random urban-like topographies is conducted to explore the effects of surface geometry on bulk aerodynamic characterization. In a mean sense, we find that statistical attributes including surface height root mean square and skewness can adequately capture the spatial heterogeneities and randomness of real urban geometries. We find, however, that higher-order statistical moments have a negligible affect on aerodynamic drag (i.e. kurtosis may be omitted). The results enable exploration of applicability of some recently-proposed roughness parameterizations that are relevant to complex, urban-like roughness (including the model proposed by Flack and Schultz, 2010: J. Fluids Eng. 132, 041203-1). We evaluate empirical parameters needed in these models for the present urban-like cases. We find that two empirical parameters (relevant to height rms and skewness) can characterize the bulk aerodynamic roughness of topographies with statistical attributes comparable to dense urban environments. This work was supported by the Army Research Office, Atmospheric Sciences Program (PM: Dr. S. Collier) under Grant # W911NF-13-1-0474. Computational resources were provided by the Texas Advanced Computing Center at the University of Texas.
NASA Astrophysics Data System (ADS)
Chamecki, M.; Pan, Y.; Nepf, H. M.; Follett, E.
2014-12-01
Flexible plants bend in response to fluid motion and this reconfiguration mechanism allows plants to minimize the increase of drag force with increasing velocity, ensuring survival in flow-dominated habitats. The effect of reconfiguration on the flow field can be modeled by introducing a drag coefficient that decreases with increasing velocity. Typically, a power-law decrease of the drag coefficient with increasing velocity is used, and the negative exponent is known as the Vogel number. In practice, the Vogel number is a function of canopy rigidity and flow conditions. In this work we show that accounting for the effect of reconfiguration is required for large-eddy simulation (LES) models to reproduce the skewness of the streamwise and vertical velocity components and the distribution of sweeps and ejections observed in a large cornfield. Additional LES runs are conducted to investigate the structure of turbulence in different reconfiguration regimes, with mean vertical momentum flux constrained by measurements. The change of the Vogel number has negligible effects on LES predictions of the total vertical momentum flux and the components of turbulent kinetic energy, but produces profound changes in the mechanisms of momentum transport. This work demonstrates the necessity to model the effect of reconfiguration in LES studies of canopy flows. It also highlights the impacts of reconfiguration on the structure of turbulence and the dynamics of momentum fluxes, as well as any other process that depends on velocity fluctuations above and within the canopy region.
NASA Astrophysics Data System (ADS)
Wang, Shouping; Stevens, Bjorn
2000-02-01
Large eddy simulation is used to study top-hat parameterizations of second- and third-order scalar statistics in cumulus and stratocumulus cloud-topped boundary layers (CTBLs). Although the top-hat parameterizations based on commonly used conditional sampling methods are a useful approach to modeling the vertical fluxes in the simulated CTBLs, they fail to realistically represent the scalar variances. The reason is that the common sampling methods are based at least in part on the sign of vertical velocity, but not on the sign of the scalars whose variances are represented and that scalars and velocity are not perfectly correlated. Furthermore, the self-correlation nature for a variance means that all the fluctuations contribute to its value, while the top-hat models completely ignore the deviations from the top-hat means and thus considerably degrade the representation of the variance. For the fluxes, however, only the coherent convective elements make the most contribution. Analysis of analytic models and `toy' time series indicates in a more generic setting that the effect of poor correlations between the signal upon which the sampling is based and the signal whose variance is to be represented tends to degrade the ability of top-hat parameterizations to capture the variance. The analysis of toy time series also indicates that variability among individual events within a composite degrades the top-hat representation of the variance more than variability within events. For the vertical velocity-scalar-related third-order moments, the top-hat model gives reasonable estimates for the cumulus CTBL but not for the stratocumulus CTBL. These differences are explained by structural differences (tied to circulation differences in the two CTBLs) in their respective joint probability density functions of vertical velocity and various scalars.
NASA Astrophysics Data System (ADS)
Zhang, Jiangshan; Li, Jingshe; Yan, Yi; Chen, Zhixin; Yang, Shufeng; Zhao, Jingwei; Jiang, Zhengyi
2016-02-01
The advantages of trumpet-shaped ladle shrouds (TLS) have been frequently demonstrated over conventional straight-bore ladle shrouds (CLS) with respect to production efficiency and molten steel quality in continuous casting practices. The present study is to shed some lights on why the TLS are better than the CLS design by examining the fluid dynamics and mass transfer using large eddy simulation. The obtained numerical results were validated with particle imaging velocimetry experiments. Flow velocity, deformation, turbulent energy dissipation, and mixing kinetics of tracer were discussed. The results showed that the entering jet of the CLS flowed straight down into the tundish with a relatively high speed (average at 0.710 to 0.815 m/s) and turbulent kinetic energy. However, the trumpet section of a TLS intensified velocity differences, strain rates, and vortices, and promoted an increase on turbulence dissipation rate in the interior of the ladle shroud. The average speed of the entering jet to the tundish was decreased to 0.270 to 0.410 m/s from the 0.708 m/s of the inlet speed. The entering jet from the TLS swung, twisted and well mixed with surrounding fluid in the tundish, and dissipated its kinetic energy. Consequently, the turbulence of the whole flow field as well as the mean skin friction coefficient of tundish wall and the velocity of free liquid surface were reduced. A tracer experiment was carried out to study mass transfer and flow mixing behavior, and the results demonstrated that the use of the TLS increased the plug volume and decreased the dead zone, thereby enhancing inclusion flotation.
NASA Technical Reports Server (NTRS)
Senocak, I.; Ackerman, A. S.; Kirkpatrick, M. P.; Stevens, D. E.; Mansour, N. N.
2004-01-01
Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.
Large Eddy Simulation of Multiple Turbulent Round Jets
NASA Astrophysics Data System (ADS)
Balajee, G. K.; Panchapakesan, Nagangudy
2015-11-01
Turbulent round jet flow was simulated as a large eddy simulation with OpenFoam software package for a jet Reynolds number of 11000. The intensity of the fluctuating motion in the incoming nozzle flow was adjusted so that the initial shear layer development compares well with available experimental data. The far field development of averages of higher order moments up to fourth order were compared with experiments. The agreement is good indicating that the large eddy motions were being computed satisfactorily by the simulation. Turbulent kinetic energy budget as well as the quality of the LES simulations were also evaluated. These conditions were then used to perform a multiple turbulent round jets simulation with the same initial momentum flux. The far field of the flow was compared with the single jet simulation and experiments to test approach to self similarity. The evolution of the higher order moments in the development region where the multiple jets interact were studied. We will also present FTLE fields computed from the simulation to educe structures and compare it with those educed by other scalar measures. Support of AR&DB CIFAAR, and VIRGO cluster at IIT Madras is gratefully acknowledged.
Large-Eddy Simulation Study of the Effects on Flow of a Heterogeneous Forest at Sub-Tree Resolution
NASA Astrophysics Data System (ADS)
Schlegel, Fabian; Stiller, Jörg; Bienert, Anne; Maas, Hans-Gerd; Queck, Ronald; Bernhofer, Christian
2015-01-01
The effect of three-dimensional plant heterogeneity on flow past a clearing is investigated by means of large-eddy simulation. A detailed representation of the canopy has been acquired by terrestrial laser scanning for a patch of approximately length and width at the field site "Tharandter Wald", near the city of Dresden, Germany. The scanning data are used to produce a highly resolved, three-dimensional plant area distribution representing the actual canopy. Hence, the vegetation maintains a rich horizontal and vertical structure including the three-dimensional clearing. The scanned plant area density is embedded in a larger domain, which is filled with a heterogeneous forest generated by the virtual canopy generator of Bohrer et al. (Tellus B 59:566-576, 2007). Based on forest inventory maps and airborne laser scanning, the characteristics of the actual canopy are preserved. Furthermore, the topography is extracted from a digital terrain model with some modifications to accommodate for periodic boundary conditions. A large-eddy simulation is performed for neutral atmospheric conditions and compared to simulations of a two-dimensional plant area density and an one-year-long field experiment conducted at the corresponding field site. The results reveal a considerable influence of the plant heterogeneity on the mean velocity field as well as on the turbulent quantities. The three-dimensional environment, e.g., the oblique edges combined with horizontal and vertical variations in plant area density and the topography create a sustained vertical and cross-flow velocity. Downstream of the windward forest edge an enhanced gust zone develops, whose intensity and relative position are influenced by the local canopy density and, therefore, is not constant along the edge. These results lead us to the conclusion that the usage of a three-dimensional plant area distribution is essential for capturing the flow features inside the canopy and within the mixing layer above.
Finecasting for renewable energy with large-eddy simulation
NASA Astrophysics Data System (ADS)
Jonker, Harmen; Verzijlbergh, Remco
2016-04-01
We present results of a single, continuous Large-Eddy Simulation of actual weather conditions during the timespan of a full year, made possible through recent computational developments (Schalkwijk et al, MWR, 2015). The simulation is coupled to a regional weather model in order to provide an LES dataset that is representative of the daily weather of the year 2012 around Cabauw, the Netherlands. This location is chosen such that LES results can be compared with both the regional weather model and observations from the Cabauw observational supersite. The run was made possible by porting our Large-Eddy Simulation program to run completely on the GPU (Schalkwijk et al, BAMS, 2012). GPU adaptation allows us to reach much improved time-to-solution ratios (i.e. simulation speedup versus real time). As a result, one can perform runs with a much longer timespan than previously feasible. The dataset resulting from the LES run provides many avenues for further study. First, it can provide a more statistical approach to boundary-layer turbulence than the more common case-studies by simulating a diverse but representative set of situations, as well as the transition between situations. This has advantages in designing and evaluating parameterizations. In addition, we discuss the opportunities of high-resolution forecasts for the renewable energy sector, e.g. wind and solar energy production.
NASA Astrophysics Data System (ADS)
Carper, Matthew A.; Porté-Agel, Fernando
2008-04-01
The ability of subfilter-scale (SFS) models to reproduce the statistical properties of SFS stresses and energy transfers over heterogeneous surface roughness is key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. In this study, several SFS models are evaluated a priori using experimental data acquired downwind of a rough-to-smooth transition in a wind tunnel. The SFS models studied include the eddy-viscosity, similarity, non-linear and a mixed model consisting of a combination of the eddy-viscosity and non-linear models. The dynamic eddy-viscosity model is also evaluated. The experimental data consist of vertical and horizontal planes of high-spatial-resolution velocity fields measured using particle image velocimetry. These velocity fields are spatially filtered and used to calculate SFS stresses and SFS transfer rates of resolved kinetic energy. Coefficients for each SFS model are calculated by matching the measured and modelled SFS energy transfer rates. For the eddy-viscosity model, the Smagorinsky coefficient is also evaluated using a dynamic procedure. The model coefficients are found to be scale dependent when the filter scales are larger than the vertical measurement height and fall into the production subrange of the turbulence where the flow scales are anisotropic. Near the surface, the Smagorinsky coefficient is also found to decrease with distance downwind from the transition, in response to the increase in mean shear as the flow adjusts to the smooth surface. In a priori tests, the ability of each model to reproduce statistical properties of the SFS stress is assessed. While the eddy-viscosity model has low spatial correlation with the measured stress, it predicts mean stresses with the same accuracy as the other models. However, the deficiency of the eddy-viscosity model is apparent in the underestimation of the standard deviation of the SFS stresses and the inability to predict transfers of kinetic energy from
Large-Eddy Simulation of Wind-Plant Aerodynamics: Preprint
Churchfield, M. J.; Lee, S.; Moriarty, P. J.; Martinez, L. A.; Leonardi, S.; Vijayakumar, G.; Brasseur, J. G.
2012-01-01
In this work, we present results of a large-eddy simulation of the 48 multi-megawatt turbines composing the Lillgrund wind plant. Turbulent inflow wind is created by performing an atmospheric boundary layer precursor simulation and turbines are modeled using a rotating, variable-speed actuator line representation. The motivation for this work is that few others have done wind plant large-eddy simulations with a substantial number of turbines, and the methods for carrying out the simulations are varied. We wish to draw upon the strengths of the existing simulations and our growing atmospheric large-eddy simulation capability to create a sound methodology for performing this type of simulation. We have used the OpenFOAM CFD toolbox to create our solver.
Large-eddy simulation of flow past a circular cylinder
NASA Technical Reports Server (NTRS)
Mittal, R.
1995-01-01
Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.
NASA Astrophysics Data System (ADS)
Ooi, Seng-Keat
2005-11-01
Lock-exchange gravity current flows produced by the instantaneous release of a heavy fluid are investigated using 3-D well resolved Large Eddy Simulation simulations at Grashof numbers up to 8*10^9. It is found the 3-D simulations correctly predict a constant front velocity over the initial slumping phase and a front speed decrease proportional to t-1/3 (the time t is measured from the release) over the inviscid phase, in agreement with theory. The evolution of the current in the simulations is found to be similar to that observed experimentally by Hacker et al. (1996). The effect of the dynamic LES model on the solutions is discussed. The energy budget of the current is discussed and the contribution of the turbulent dissipation to the total dissipation is analyzed. The limitations of less expensive 2D simulations are discussed; in particular their failure to correctly predict the spatio-temporal distributions of the bed shear stresses which is important in determining the amount of sediment the gravity current can entrain in the case in advances of a loose bed.
NASA Astrophysics Data System (ADS)
Cheng, Wai Chi; Porté-Agel, Fernando
2014-05-01
Accurate prediction of atmospheric boundary layer (ABL) flow and its interaction with urban surfaces is critical for understanding the transport of momentum and scalars within and above cities. This, in turn, is essential for predicting the local climate and pollutant dispersion patterns in urban areas. Large-eddy simulation (LES) explicitly resolves the large-scale turbulent eddy motions and, therefore, can potentially provide improved understanding and prediction of flows inside and above urban canopies. This study focuses on the validation and the use of a recently-developed LES framework to simulate a turbulent boundary layer flow through idealized urban canopies represented by uniform arrays of cubes. The LES framework is first validated with wind tunnel experimental data. Good agreement between the simulation results and the experimental data are found in the vertical and spanwise profiles of mean velocities and velocity standard deviations at different streamwise locations. Next, the model is used to simulate ABL flows over surface transitions from a flat homogeneous terrain to aligned and staggered arrays of cubes with height h. For both configurations, five different frontal area densities (Λf), equal to 0.028, 0.063, 0.111, 0.174 and 0.250, are considered. Within the arrays, the flow is found to adjust quickly and shows similar structure of the wake of the cubes after the second row. Above the arrays, an internal boundary layer (IBL) is identified. No significant difference in the depth of the IBL among different cases is observed. The drag exerted by the cubes on the flow (Df) and the drag coefficients of the cubes (Cd) are calculated explicitly using the LES results. For the downstream cubes, Df is found to increases with decreasing density for both configurations, and larger values of Cd are found for the cubes of staggered arrays than those of the aligned arrays with the same Λf. At a downstream location where the flow immediately above the cube
NASA Astrophysics Data System (ADS)
Wang, BoBin; Wang, ZhiShi; Cui, GuiXiang; Zhang, ZhaoShun
2014-06-01
In this paper, the dynamic characteristics of building clusters are simulated by large eddy simulation at high Reynolds number for both homogeneous and heterogeneous building clusters. To save the computational cost a channel-like flow model is applied to the urban canopy with free slip condition at the upper boundary. The results show that the domain height is an important parameter for correct evaluation of the dynamic characteristics. The domain height must be greater than 8 h ( h is the average building height) in order to obtain correct roughness height while displacement height and roughness sublayer are less sensitive to the domain height. The Reynolds number effects on the dynamic characteristics and flow patterns are investigated. The turbulence intensity is stronger inside building cluster at high Reynolds number while turbulence intensity is almost unchanged with Reynolds number above the building cluster. Roughness height increases monotonously with Reynolds number by 20% from Re*=103 to Re*=105 but displacement height is almost unchanged. Within the canopy layer of heterogeneous building clusters, flow structures vary between buildings and turbulence is more active at high Reynolds number.
NASA Astrophysics Data System (ADS)
Margairaz, Fabien; Giometto, Marco; Parlange, Marc; Calaf, Marc
2015-11-01
The performance of dealiasing schemes and their computational cost on a pseudo-spectral code are analyzed. Dealiasing is required to limit the error that occurs when two discretized variables are multiplied, polluting the accuracy of the result. In this work three different dealiasing methods are explored: the 2/3 rule, the 3/2 rule, and a high order Fourier smoothing based method. We compare the cost of the traditionally accepted 3/2 rule (Canuto et al., 1988), where an expansion of the computational domain to a larger grid is required, to the cost of the other two techniques that do not require this expansion. This analysis is performed in the framework of Large-Eddy Simulations (LES) of incompressible flows using the constant Smagorinsky sub-grid model with a wall damping function and a wall model based on the log-law. A highly efficient LES code parallelized using a 2D pencil decomposition has been developed. The code employs the traditional pseudo-spectral approach to integrate the incompressible Navier-Stokes equations. Several simulations of a neutral atmospheric boundary layer using different degrees of numerical resolution are considered. Results show a net difference in computational cost between the different techniques without relevant changes in statistics.
Large Eddy Simulations and Turbulence Modeling for Film Cooling
NASA Technical Reports Server (NTRS)
Acharya, Sumanta
1999-01-01
The objective of the research is to perform Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) for film cooling process, and to evaluate and improve advanced forms of the two equation turbulence models for turbine blade surface flow analysis. The DNS/LES were used to resolve the large eddies within the flow field near the coolant jet location. The work involved code development and applications of the codes developed to the film cooling problems. Five different codes were developed and utilized to perform this research. This report presented a summary of the development of the codes and their applications to analyze the turbulence properties at locations near coolant injection holes.
Large-Eddy Simulation of Wind-Plant Aerodynamics
Churchfield, M. J.; Lee, S.; Moriarty, P. J.; Martinez, L. A.; Leonardi, S.; Vijayakumar, G.; Brasseur, J. G.
2012-01-01
In this work, we present results of a large-eddy simulation of the 48 multi-megawatt turbines composing the Lillgrund wind plant. Turbulent inflow wind is created by performing an atmospheric boundary layer precursor simulation, and turbines are modeled using a rotating, variable-speed actuator line representation. The motivation for this work is that few others have done large-eddy simulations of wind plants with a substantial number of turbines, and the methods for carrying out the simulations are varied. We wish to draw upon the strengths of the existing simulations and our growing atmospheric large-eddy simulation capability to create a sound methodology for performing this type of simulation. We used the OpenFOAM CFD toolbox to create our solver. The simulated time-averaged power production of the turbines in the plant agrees well with field observations, except with the sixth turbine and beyond in each wind-aligned. The power produced by each of those turbines is overpredicted by 25-40%. A direct comparison between simulated and field data is difficult because we simulate one wind direction with a speed and turbulence intensity characteristic of Lillgrund, but the field observations were taken over a year of varying conditions. The simulation shows the significant 60-70% decrease in the performance of the turbines behind the front row in this plant that has a spacing of 4.3 rotor diameters in this direction. The overall plant efficiency is well predicted. This work shows the importance of using local grid refinement to simultaneously capture the meter-scale details of the turbine wake and the kilometer-scale turbulent atmospheric structures. Although this work illustrates the power of large-eddy simulation in producing a time-accurate solution, it required about one million processor-hours, showing the significant cost of large-eddy simulation.
Large-eddy simulation using the finite element method
McCallen, R.C.; Gresho, P.M.; Leone, J.M. Jr.; Kollmann, W.
1993-10-01
In a large-eddy simulation (LES) of turbulent flows, the large-scale motion is calculated explicitly (i.e., approximated with semi-empirical relations). Typically, finite difference or spectral numerical schemes are used to generate an LES; the use of finite element methods (FEM) has been far less prominent. In this study, we demonstrate that FEM in combination with LES provides a viable tool for the study of turbulent, separating channel flows, specifically the flow over a two-dimensional backward-facing step. The combination of these methodologies brings together the advantages of each: LES provides a high degree of accuracy with a minimum of empiricism for turbulence modeling and FEM provides a robust way to simulate flow in very complex domains of practical interest. Such a combination should prove very valuable to the engineering community.
Large eddy simulation using the general circulation model ICON
NASA Astrophysics Data System (ADS)
Dipankar, Anurag; Stevens, Bjorn; Heinze, Rieke; Moseley, Christopher; Zängl, Günther; Giorgetta, Marco; Brdar, Slavko
2015-09-01
ICON (ICOsahedral Nonhydrostatic) is a unified modeling system for global numerical weather prediction (NWP) and climate studies. Validation of its dynamical core against a test suite for numerical weather forecasting has been recently published by Zängl et al. (2014). In the present work, an extension of ICON is presented that enables it to perform as a large eddy simulation (LES) model. The details of the implementation of the LES turbulence scheme in ICON are explained and test cases are performed to validate it against two standard LES models. Despite the limitations that ICON inherits from being a unified modeling system, it performs well in capturing the mean flow characteristics and the turbulent statistics of two simulated flow configurations—one being a dry convective boundary layer and the other a cumulus-topped planetary boundary layer.
Large-eddy simulation of turbulent circular jet flows
Jones, S. C.; Sotiropoulos, F.; Sale, M. J.
2002-07-01
This report presents a numerical method for carrying out large-eddy simulations (LES) of turbulent free shear flows and an application of a method to simulate the flow generated by a nozzle discharging into a stagnant reservoir. The objective of the study was to elucidate the complex features of the instantaneous flow field to help interpret the results of recent biological experiments in which live fish were exposed to the jet shear zone. The fish-jet experiments were conducted at the Pacific Northwest National Laboratory (PNNL) under the auspices of the U.S. Department of Energy’s Advanced Hydropower Turbine Systems program. The experiments were designed to establish critical thresholds of shear and turbulence-induced loads to guide the development of innovative, fish-friendly hydropower turbine designs.
Large Eddy Simulation of Vertical Axis Wind Turbine Wakes
NASA Astrophysics Data System (ADS)
Shamsoddin, Sina; Porté-Agel, Fernando
2014-05-01
In this study, large-eddy simulation (LES) is combined with a turbine model to investigate the wake behind a vertical-axis wind turbine (VAWT) in a three dimensional turbulent flow. Two methods are used to model the subgrid-scale (SGS) stresses: (a) the Smagorinsky model, and (b) the modulated gradient model. To parameterize the effects of the VAWT on the flow, two VAWT models are developed: (a) the actuator surface model (ASM), in which the time-averaged turbine-induced forces are distributed on a surface swept by the turbine blades, i.e. the actuator surface, and (b) the actuator line model (ALM), in which the instantaneous blade forces are only spatially distributed on lines representing the blades, i.e. the actuator lines. This is the first time that LES is applied and validated for simulation of VAWT wakes by using either the ASM or the ALM techniques. In both models, blade-element theory is used to calculate the lift and drag forces on the blades. The results are compared with flow measurements in the wake of a model straight-bladed VAWT, carried out in the Institute de Méchanique et Statistique de la Turbulence (IMST) water channel. Different combinations of SGS models with VAWT models are studied and a fairly good overall agreement between simulation results and measurement data is observed. In general, the ALM is found to better capture the unsteady-periodic nature of the wake and shows a better agreement with the experimental data compared with the ASM. The modulated gradient model is also found to be a more reliable SGS stress modeling technique, compared with the Smagorinsky model, and it yields reasonable predictions of the mean flow and turbulence characteristics of a VAWT wake using its theoretically-determined model coefficient. Keywords: Vertical-axis wind turbines (VAWTs); VAWT wake; Large-eddy simulation; Actuator surface model; Actuator line model; Smagorinsky model; Modulated gradient model
NASA Astrophysics Data System (ADS)
Dipankar, A.; Stevens, B. B.; Zängl, G.; Pondkule, M.; Brdar, S.
2014-12-01
The effect of clouds on large scale dynamics is represented in climate models through parameterization of various processes, of which the parameterization of shallow and deep convection are particularly uncertain. The atmospheric boundary layer, which controls the coupling to the surface, and which defines the scale of shallow convection, is typically 1 km in depth. Thus, simulations on a O(100 m) grid largely obviate the need for such parameterizations. By crossing this threshold of O(100m) grid resolution one can begin thinking of large-eddy simulation (LES), wherein the sub-grid scale parameterization have a sounder theoretical foundation. Substantial initiatives have been taken internationally to approach this threshold. For example, Miura et al., 2007 and Mirakawa et al., 2014 approach this threshold by doing global simulations, with (gradually) decreasing grid resolution, to understand the effect of cloud-resolving scales on the general circulation. Our strategy, on the other hand, is to take a big leap forward by fixing the resolution at O(100 m), and gradually increasing the domain size. We believe that breaking this threshold would greatly help in improving the parameterization schemes and reducing the uncertainty in climate predictions. To take this forward, the German Federal Ministry of Education and Research has initiated a project on HD(CP)2 that aims for a limited area LES at resolution O(100 m) using the new unified modeling system ICON (Zängl et al., 2014). In the talk, results from the HD(CP)2 evaluation simulation will be shown that targets high resolution simulation over a small domain around Jülich, Germany. This site is chosen because high resolution HD(CP)2 Observational Prototype Experiment took place in this region from 1.04.2013 to 31.05.2013, in order to critically evaluate the model. Nesting capabilities of ICON is used to gradually increase the resolution from the outermost domain, which is forced from the COSMO-DE data, to the
NASA's Large-Eddy Simulation Research for Jet Noise Applications
NASA Technical Reports Server (NTRS)
DeBonis, James R.
2009-01-01
Research into large-eddy simulation (LES) for application to jet noise is described. The LES efforts include in-house code development and application at NASA Glenn along with NASA Research Announcement sponsored work at Stanford University and Florida State University. Details of the computational methods used and sample results for jet flows are provided.
Large-eddy simulation in complex domains using the finite element method
McCallen, R.C.; Kornblum, B.T.; Kollman, W.
1996-11-12
Finite element methods (FEM) are demonstrated in combination with large-eddy simulations (LES) as a valuable tool for the study of turbulent, separating channel flows, specifically the flow over a backward facing step.
Large-eddy simulation of cavitating nozzle and jet flows
NASA Astrophysics Data System (ADS)
Örley, F.; Trummler, T.; Hickel, S.; Mihatsch, M. S.; Schmidt, S. J.; Adams, N. A.
2015-12-01
We present implicit large-eddy simulations (LES) to study the primary breakup of cavitating liquid jets. The considered configuration, which consists of a rectangular nozzle geometry, adopts the setup of a reference experiment for validation. The setup is a generic reproduction of a scaled-up automotive fuel injector. Modelling of all components (i.e. gas, liquid, and vapor) is based on a barotropic two-fluid two-phase model and employs a homogenous mixture approach. The cavitating liquid model assumes thermodynamic- equilibrium. Compressibility of all phases is considered in order to capture pressure wave dynamics of collapse events. Since development of cavitation significantly affects jet break-up characteristics, we study three different operating points. We identify three main mechanisms which induce primary jet break-up: amplification of turbulent fluctuations, gas entrainment, and collapse events near the liquid-gas interface.
NASA Astrophysics Data System (ADS)
Lew, Phoi-Tack
Improvements in computing speed over the past decade have made Large Eddy Simulations (LES) amenable to the study of jet noise. The study of turbulent hot jets is required jets since all jet engines fitted on aircraft operate at hot exhaust conditions. The primary goal of this research was to further advance the science of jet noise prediction with a specific emphasis on heated jets using 3-D LES. For the 3-D LES methodology, spatial filtering is used as an implicit subgrid scale (SGS) model in place of an explicit SGS model, such as the classical Smagorinsky or Dynamic Smagorinsky models. To study the far-field noise, the porous FfowcsWilliams-Hawkings (FWH) surface integral acoustic formulation is employed. Results obtained for the heated jets in terms of jet development are in good agreement with other LES results and experimental data. The predicted overall sound pressure level (OASPL) values for heated jets exhibited the same trend as experimental data. The levels were over-predicted by approximately 3 dB, which was deemed satisfactory. An investigation of noise sources for heated jets was also performed within the framework of Lighthill's acoustic analogy. It is discovered that when a high-speed is jet heated, significant cancellations occur between shear and entropy noise sources compared to an unheated high speed jet. This could explain why a high speed heated jet is quieter than an unheated jet at the same ambient Mach number. High-order compact finite difference schemes along with high-order filters are used extensively in LES, especially for aeroacoustics problems, since these schemes have very high accuracy and spectral-like resolution as well as low-dispersion and diffusion errors. Due to the implicit nature of compact schemes, one technique of parallelization is based on the data transposition strategy. However, such transposition strategy is near impossible to apply to jets with complex geometries. Hence, an alternative parallelization methodology
Large eddy simulations of a forced semiconfined circular impinging jet
NASA Astrophysics Data System (ADS)
Olsson, M.; Fuchs, L.
1998-02-01
Large eddy simulations (LES) of a forced semiconfined circular impinging jet were carried out. The Reynolds number was 104 and the inflow was forced at a Strouhal number of 0.27. The separation between the jet inlet and the opposing wall was four jet inlet diameters. Four different simulations were made. Two simulations were performed without any explicit sub-grid-scale (SGS) model using 1283 and 963 grid points, respectively. Two simulations were performed with two different SGS-models using 963 grid points; one with a dynamic Smagorinsky based model and one with a stress-similarity model. The simulations were performed to study the mean velocity, the turbulence statistics, the SGS-model effects, the dynamic behavior of the jet with a focus on the near wall region. The existence of separation vortices in the wall jet region was confirmed. These secondary vortices were found to be related to the radially deflected primary vortices generated by the circular shear layer of the jet. It was also shown that the primary vortex structures that reach the wall were helical and not axisymmetric. A quantitative gain was found in the simulations with SGS-models. The stress-similarity model simulation correlated slightly better with the higher resolution simulation than the other coarse grid simulations. The variations in the results predicted by the different simulations were larger for the turbulence statistics than for the mean velocity. However, the variation among the different simulations in terms of the turbulence intensity was less than 10%.
NASA Astrophysics Data System (ADS)
Cheng, Wai-Chi; Porté-Agel, Fernando
2015-05-01
Large-eddy simulations (LES) are performed to simulate the atmospheric boundary-layer (ABL) flow through idealized urban canopies represented by uniform arrays of cubes in order to better understand atmospheric flow over rural-to-urban surface transitions. The LES framework is first validated with wind-tunnel experimental data. Good agreement between the simulation results and the experimental data are found for the vertical and spanwise profiles of the mean velocities and velocity standard deviations at different streamwise locations. Next, the model is used to simulate ABL flows over surface transitions from a flat homogeneous terrain to aligned and staggered arrays of cubes with height . For both configurations, five different frontal area densities , equal to 0.028, 0.063, 0.111, 0.174 and 0.250, are considered. Within the arrays, the flow is found to adjust quickly and shows similar structure to the wake of the cubes after the second row of cubes. An internal boundary layer is identified above the cube arrays and found to have a similar depth in all different cases. At a downstream location where the flow immediately above the cube array is already adjusted to the surface, the spatially-averaged velocity is found to have a logarithmic profile in the vertical. The values of the displacement height are found to be quite insensitive to the canopy layout (aligned vs. staggered) and increase roughly from to as increases from 0.028 to 0.25. Relatively larger values of the aerodynamic roughness length are obtained for the staggered arrays, compared with the aligned cases, and a maximum value of is found at for both configurations. By explicitly calculating the drag exerted by the cubes on the flow and the drag coefficients of the cubes using our LES results, and comparing the results with existing theoretical expressions, we show that the larger values of for the staggered arrays are related to the relatively larger drag coefficients of the cubes for that
Large-eddy simulation of turbulence in steam generators
Bagwell, T.G.; Hassan, Y.A. ); Steininger D.A. )
1989-11-01
A major problem associated with steam generators is excessive tube vibration caused by turbulent-flow buffeting and fluid-elastic excitation. Vibration can lead to tube rupture or wear, necessitating tube plugging and reducing the availability of the steam generator. The fluid/structure interaction phenomenon that causes fluid-elastic tube excitation is unknown at present. The current investigation defines the spectral characteristics of turbulent flow entering the Westinghouse D4 steam generator tube bundles using the large-eddy simulation (LES) technique. Due to the recent availability of supercomputers, LES is being considered as a possible engineering design analysis tool. The information from this study will provide input for defining the temporally fluctuating forces on steam generator tube banks. The GUST code was used to analyze the water box of a Westinghouse model D4 steam generator.
Constrained Large Eddy Simulation of Separated Turbulent Flows
NASA Astrophysics Data System (ADS)
Xia, Zhenhua; Shi, Yipeng; Wang, Jianchun; Xiao, Zuoli; Yang, Yantao; Chen, Shiyi
2011-11-01
Constrained Large-eddy Simulation (CLES) has been recently proposed to simulate turbulent flows with massive separation. Different from traditional large eddy simulation (LES) and hybrid RANS/LES approaches, the CLES simulates the whole flow domain by large eddy simulation while enforcing a RANS Reynolds stress constraint on the subgrid-scale (SGS) stress models in the near-wall region. Algebraic eddy-viscosity models and one-equation Spalart-Allmaras (S-A) model have been used to constrain the Reynolds stress. The CLES approach is validated a posteriori through simulation of flow past a circular cylinder and periodic hill flow at high Reynolds numbers. The simulation results are compared with those from RANS, DES, DDES and other available hybrid RANS/LES methods. It is shown that the capability of the CLES method in predicting separated flows is comparable to that of DES. Detailed discussions are also presented about the effects of the RANS models as constraint in the near-wall layers. Our results demonstrate that the CLES method is a promising alternative towards engineering applications.
Toward the large-eddy simulation of compressible turbulent flows
NASA Technical Reports Server (NTRS)
Erlebacher, G.; Hussaini, M. Y.; Speziale, C. G.; Zang, T. A.
1990-01-01
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.
Large Eddy Simulations using Lattice Boltzmann algorithms. Final report
Serling, J.D.
1993-09-28
This report contains the results of a study performed to implement eddy-viscosity models for Large-Eddy-Simulations (LES) into Lattice Boltzmann (LB) algorithms for simulating fluid flows. This implementation requires modification of the LB method of simulating the incompressible Navier-Stokes equations to allow simulation of the filtered Navier-Stokes equations with some subgrid model for the Reynolds stress term. We demonstrate that the LB method can indeed be used for LES by simply locally adjusting the value of the BGK relaxation time to obtain the desired eddy-viscosity. Thus, many forms of eddy-viscosity models including the standard Smagorinsky model or the Dynamic model may be implemented using LB algorithms. Since underresolved LB simulations often lead to instability, the LES model actually serves to stabilize the method. An alternative method of ensuring stability is presented which requires that entropy increase during the collision step of the LB method. Thus, an alternative collision operator is locally applied if the entropy becomes too low. This stable LB method then acts as an LES scheme that effectively introduces its own eddy viscosity to damp short wavelength oscillations.
Toward large eddy simulation of turbulent flow over an airfoil
NASA Technical Reports Server (NTRS)
Choi, Haecheon
1993-01-01
The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.
Large Eddy Simulation of Flow and Sediment Transport over Dunes
NASA Astrophysics Data System (ADS)
Agegnehu, G.; Smith, H. D.
2012-12-01
Understanding the nature of flow over bedforms has a great importance in fluvial and coastal environments. For example, a bedform is one source of energy dissipation in water waves outside the surf zone in coastal environments. In rivers, the migration of dunes often affects the stability of the river bed and banks. In general, when a fluid flows over a sediment bed, the sediment transport generated by the interaction of the flow field with the bed results in the periodic deformation of the bed in the form of dunes. Dunes generally reach an equilibrium shape, and slowly propagate in the direction of the flow, as sand is lifted in the high shear regions, and redeposited in the separated flow areas. Different numerical approaches have been used in the past to study the flow and sediment transport over bedforms. In most research works, Reynolds Averaged Navier Stokes (RANS) equations are employed to study fluid motions over ripples and dunes. However, evidences suggests that these models can not represent key turbulent quantities in unsteady boundary layers. The use of Large Eddy Simulation (LES) can resolve a much larger range of smaller scales than RANS. Moreover, unsteady simulations using LES give vital turbulent quantities which can help to study fluid motion and sediment transport over dunes. For this steady, we use a three-dimensional, non-hydrostatic model, OpenFOAM. It is a freely available tool which has different solvers to simulate specific problems in engineering and fluid mechanics. Our objective is to examine the flow and sediment transport from numerical stand point for bed geometries that are typical of fixed dunes. At the first step, we performed Large Eddy Simulation of the flow over dune geometries based on the experimental data of Nelson et al. (1993). The instantaneous flow field is investigated with special emphasis on the occurrence of coherent structures. To assess the effect of bed geometries on near bed turbulence, we considered different
Large eddy simulation of a pumped- storage reservoir
NASA Astrophysics Data System (ADS)
Launay, Marina; Leite Ribeiro, Marcelo; Roman, Federico; Armenio, Vincenzo
2016-04-01
The last decades have seen an increasing number of pumped-storage hydropower projects all over the world. Pumped-storage schemes move water between two reservoirs located at different elevations to store energy and to generate electricity following the electricity demand. Thus the reservoirs can be subject to important water level variations occurring at the daily scale. These new cycles leads to changes in the hydraulic behaviour of the reservoirs. Sediment dynamics and sediment budgets are modified, sometimes inducing problems of erosion and deposition within the reservoirs. With the development of computer performances, the use of numerical techniques has become popular for the study of environmental processes. Among numerical techniques, Large Eddy Simulation (LES) has arisen as an alternative tool for problems characterized by complex physics and geometries. This work uses the LES-COAST Code, a LES model under development in the framework of the Seditrans Project, for the simulation of an Upper Alpine Reservoir of a pumped-storage scheme. Simulations consider the filling (pump mode) and emptying (turbine mode) of the reservoir. The hydraulic results give a better understanding of the processes occurring within the reservoir. They are considered for an assessment of the sediment transport processes and of their consequences.
Large eddy simulation of boundary layer flow under cnoidal waves
NASA Astrophysics Data System (ADS)
Li, Yin-Jun; Chen, Jiang-Bo; Zhou, Ji-Fu; Zhang, Qiang
2016-02-01
Water waves in coastal areas are generally nonlinear, exhibiting asymmetric velocity profiles with different amplitudes of crest and trough. The behaviors of the boundary layer under asymmetric waves are of great significance for sediment transport in natural circumstances. While previous studies have mainly focused on linear or symmetric waves, asymmetric wave-induced flows remain unclear, particularly in the flow regime with high Reynolds numbers. Taking cnoidal wave as a typical example of asymmetric waves, we propose to use an infinite immersed plate oscillating cnoidally in its own plane in quiescent water to simulate asymmetric wave boundary layer. A large eddy simulation approach with Smagorinsky subgrid model is adopted to investigate the flow characteristics of the boundary layer. It is verified that the model well reproduces experimental and theoretical results. Then a series of numerical experiments are carried out to study the boundary layer beneath cnoidal waves from laminar to fully developed turbulent regimes at high Reynolds numbers, larger than ever studied before. Results of velocity profile, wall shear stress, friction coefficient, phase lead between velocity and wall shear stress, and the boundary layer thickness are obtained. The dependencies of these boundary layer properties on the asymmetric degree and Reynolds number are discussed in detail.
NASA Astrophysics Data System (ADS)
Guda, Venkata Subba Sai Satish
There have been several advancements in the aerospace industry in areas of design such as aerodynamics, designs, controls and propulsion; all aimed at one common goal i.e. increasing efficiency --range and scope of operation with lesser fuel consumption. Several methods of flow control have been tried. Some were successful, some failed and many were termed as impractical. The low Reynolds number regime of 104 - 105 is a very interesting range. Flow physics in this range are quite different than those of higher Reynolds number range. Mid and high altitude UAV's, MAV's, sailplanes, jet engine fan blades, inboard helicopter rotor blades and wind turbine rotors are some of the aerodynamic applications that fall in this range. The current study deals with using dynamic roughness as a means of flow control over a NACA 0012 airfoil at low Reynolds numbers. Dynamic 3-D surface roughness elements on an airfoil placed near the leading edge aim at increasing the efficiency by suppressing the effects of leading edge separation like leading edge stall by delaying or totally eliminating flow separation. A numerical study of the above method has been carried out by means of a Large Eddy Simulation, a mathematical model for turbulence in Computational Fluid Dynamics, owing to the highly unsteady nature of the flow. A user defined function has been developed for the 3-D dynamic roughness element motion. Results from simulations have been compared to those from experimental PIV data. Large eddy simulations have relatively well captured the leading edge stall. For the clean cases, i.e. with the DR not actuated, the LES was able to reproduce experimental results in a reasonable fashion. However DR simulation results show that it fails to reattach the flow and suppress flow separation compared to experiments. Several novel techniques of grid design and hump creation are introduced through this study.
Large eddy simulation of controlled transition to turbulence
NASA Astrophysics Data System (ADS)
Sayadi, Taraneh; Moin, Parviz
2012-11-01
Large eddy simulation of H- and K-type transitions in a spatially developing zero-pressure-gradient boundary layer at Ma∞ = 0.2 is investigated using several subgrid scale (SGS) models including constant coefficient Smagorinsky and Vreman models and their dynamic extensions, dynamic mixed scale-similarity, dynamic one-equation kinetic energy model, and global coefficient Vreman models. A key objective of this study is to assess the capability of SGS models to predict the location of transition and the skin friction throughout the transition process. The constant coefficient models fail to detect transition, but the dynamic procedure allows for a negligible turbulent viscosity in the early transition region. As a result, the "point" of transition is estimated correctly. However, after secondary instabilities set in and result in the overshoot in the skin friction profile, all models fail to produce sufficient subgrid scale shear stress required for the correct prediction of skin friction and the mean velocity profile. The same underprediction of skin friction persists into the turbulent region. Spatially filtered direct numerical simulation data in the same boundary layers are used to provide guidelines for SGS model development and validation.
Large eddy simulation of flame flashback in a turbulent channel
NASA Astrophysics Data System (ADS)
Hassanaly, Malik; Lietz, Christopher; Raman, Venkat; Kolla, Hemanth; Chen, Jacqueline; Gruber, Andrea; Computational Flow Physics Group Team
2014-11-01
In high-hydrogen content gas turbines, the propagation of a premixed flame along with boundary layers on the combustor walls is a source of failure, whereby the flame could enter the fuel-air premixing region that is not designed to hold high-temperature fluid. In order to develop models for predicting this phenomenon, a large eddy simulation (LES) based study is carried out here. The flow configuration is based on a direct numerical simulation (DNS) of a turbulent channel, where an initial planar flame is allowed to propagate upstream in a non-periodic channel. The LES approach uses a flamelet-based combustion model along with standard models for the unresolved subfilter flux terms. It is found that the LES are very accurate in predicting the structure of the turbulent flame front. However, there was a large discrepancy for the transient evolution of the flame, indicating that the flame-boundary layer interaction modulates flame propagation significantly, and the near-wall flame behavior may be non-flamelet like due to the anisotropic of the flow in this region.
A subfilter-scale stress model for large eddy simulations
NASA Astrophysics Data System (ADS)
Rouhi, Amirreza; Piomelli, Ugo
2013-11-01
In most large eddy simulations, the filter width is related to the grid. This method of specification, however, causes problems in complex flows where local refinement results in grid discontinuities. Following the work of Piomelli and Geurts (Proce. 8th Workshop on DLES, 2010) we propose an eddy-viscosity approach in which the filter width is based on the flow parameters only, with no explicit relationship to the grid size. This model can achieve grid-independent LES solutions, vanishing dynamically in the regions of low turbulence activity and a computational cost less than the dynamic models. The Successive Inverse Polynomial Interpolation (Geurts & Meyers Phys. Fluids 18, 2006) was used to calculate the model parameter. Calculating implicitly the eddy-viscosity at each time-step removes the numerical instabilities found in previous studies, while maintaining the local character of the model. Results of simulations of channel flow at Reτ up to 2,000, and forced homogeneous isotropic turbulence will be presented.
Large Eddy Simulation of a Sooting Jet Diffusion Flame
NASA Astrophysics Data System (ADS)
Blanquart, Guillaume; Pitsch, Heinz
2007-11-01
The understanding of soot particle dynamics in combustion systems is a key issue in the development of low emission engines. Of particular importance are the processes shaping the soot particle size distribution function (PSDF). However, it is not always necessary to represent exactly the full distribution, and often information about its moments only is sufficient. The Direct Quadrature Method of Moments (DQMOM) allows for an efficient and accurate prediction of the moments of the soot PSDF. This method has been validated for laminar premixed and diffusion flames with detailed chemistry and is now implemented in a semi-implicit low Mach-number Navier-Stokes solver. A Large Eddy Simulation (LES) of a piloted sooting jet diffusion flame (Delft flame) is performed to study the dynamics of soot particles in a turbulent environment. The profiles of temperature and major species are compared with the experimental measurements. Soot volume fraction profiles are compared with the recent data of Qamar et al. (2007). Aggregate properties such as the diameter and the fractal shape are studied in the scope of DQMOM.
Large-eddy simulations of contrails in a turbulent atmosphere
NASA Astrophysics Data System (ADS)
Picot, J.; Paoli, R.; Thouron, O.; Cariolle, D.
2014-11-01
In this work, the evolution of contrails in the vortex and dissipation regimes is studied by means of fully three-dimensional large-eddy simulation (LES) coupled to a Lagrangian particle tracking method to treat the ice phase. This is the first paper where fine-scale atmospheric turbulence is generated and sustained by means of a stochastic forcing that mimics the properties of stably stratified turbulent flows as those occurring in the upper troposphere lower stratosphere. The initial flow-field is composed by the turbulent background flow and a wake flow obtained from separate LES of the jet regime. Atmospheric turbulence is the main driver of the wake instability and the structure of the resulting wake is sensitive to the intensity of the perturbations, primarily in the vertical direction. A stronger turbulence accelerates the onset of the instability, which results in shorter contrail decent and more effective mixing in the interior of the plume. However, the self-induced turbulence that is produced in the wake after the vortex break-up dominates over background turbulence at the end of the vortex regime and dominates the mixing with ambient air. This results in global microphysical characteristics such as ice mass and optical depth that are be slightly affected by the intensity of atmospheric turbulence. On the other hand, the background humidity and temperature have a first order effect on the survival of ice crystals and particle size distribution, which is in line with recent and ongoing studies in the literature.
Large-eddy simulation of trans- and supercritical injection
NASA Astrophysics Data System (ADS)
Müller, H.; Niedermeier, C. A.; Jarczyk, M.; Pfitzner, M.; Hickel, S.; Adams, N. A.
2016-07-01
In a joint effort to develop a robust numerical tool for the simulation of injection, mixing, and combustion in liquid rocket engines at high pressure, a real-gas thermodynamics model has been implemented into two computational fluid dynamics (CFD) codes, the density-based INCA and a pressure-based version of OpenFOAM. As a part of the validation process, both codes have been used to perform large-eddy simulations (LES) of trans- and supercritical nitrogen injection. Despite the different code architecture and the different subgrid scale turbulence modeling strategy, both codes yield similar results. The agreement with the available experimental data is good.
Applications of large eddy simulation methods to gyrokinetic turbulence
Bañón Navarro, A. Happel, T.; Teaca, B. [Applied Mathematics Research Centre, Coventry University, Coventry CV1 5FB; Max-Planck für Sonnensystemforschung, Max-Planck-Str. 2, D-37191 Katlenburg-Lindau; Max-Planck Jenko, F. [Max-Planck-Institut für Plasmaphysik, EURATOM Association, D-85748 Garching; Max-Planck Hammett, G. W. [Max-Planck Collaboration: ASDEX Upgrade Team
2014-03-15
The large eddy simulation (LES) approach—solving numerically the large scales of a turbulent system and accounting for the small-scale influence through a model—is applied to nonlinear gyrokinetic systems that are driven by a number of different microinstabilities. Comparisons between modeled, lower resolution, and higher resolution simulations are performed for an experimental measurable quantity, the electron density fluctuation spectrum. Moreover, the validation and applicability of LES is demonstrated through a series of diagnostics based on the free energetics of the system.
NASA Astrophysics Data System (ADS)
Huang, G.; Newchurch, M.; Kuang, S.; Wang, L.; Ouwersloot, H.
2014-12-01
We investigate the diurnal variation of mixed-layer ozone in Huntsville AL, Southeast United States on September, 6, 2013 during the SEAC4RS field campaign. The dynamics and chemistry of the mixed layer are studied with a Large-Eddy Simulation model coupled with a chemical module and Ozone DIAL observations. In this study, we will present calculations of ozone entrainment fluxes using continuous observation by co-located ozone DIAL and Compact Wind Aerosol Lidar (CWAL) at the campus of University of Alabama in Huntsville (UAH). As a part of Tropospheric Ozone Lidar NETwork (TOLNET), UAH ozone DIAL can provide continuous ozone observation in the altitude range from 125 m AGL to 12 km, with 10-min temporal resolution and 150 - 550 m vertical resolution. We also perform an ozone budget study using the Dutch Atmospheric Large-Eddy Simulation (DALES), reasonable approximations of dry deposition, in conjunction with ozone entrainment observations. In this case study, the enhancement of ozone in the mixed layer results from the local emissions of NOx and VOCs. The NOx and VOCs emitted at surface entered into mixed layer by atmospheric turbulence and produced ozone within the whole mixed layer. Simultaneously, non-turbulent air in the residual layer, which is at top of the morning mixed layer, participates in convective mixing through entrainment processes. The clean air in the residual layer decreases the ozone enhancement rate in the mixed layer. After the mixed layer reaches its stable height at 1700m, the large-scale subsidence not only decreases PBL growth but also enhances the entrainment process. The PBL NOx and VOCs mix into the free troposphere through detrainment before producing ozone by photochemical reaction. We have following conclusions from this case study: 1) the relationship between boundary layer height and PBL ozone is complicated. Higher PBL height does not always mean lower PBL ozone. 2) The LES calculation results illuminate the interaction between
Final Report: "Large-Eddy Simulation of Anisotropic MHD Turbulence"
Zikanov, Oleg
2008-06-23
To acquire better understanding of turbulence in flows of liquid metals and other electrically conducting fluids in the presence of steady magnetic fields and to develop an accurate and physically adequate LES (large-eddy simulation) model for such flows. The scientific objectives formulated in the project proposal have been fully completed. Several new directions were initiated and advanced in the course of work. Particular achievements include a detailed study of transformation of turbulence caused by the imposed magnetic field, development of an LES model that accurately reproduces this transformation, and solution of several fundamental questions of the interaction between the magnetic field and fluid flows. Eight papers have been published in respected peer-reviewed journals, with two more papers currently undergoing review, and one in preparation for submission. A post-doctoral researcher and a graduate student have been trained in the areas of MHD, turbulence research, and computational methods. Close collaboration ties have been established with the MHD research centers in Germany and Belgium.
Large-eddy simulation of pulverized coal swirl jet flame
NASA Astrophysics Data System (ADS)
Muto, Masaya; Watanabe, Hiroaki; Kurose, Ryoichi; Komori, Satoru; Balusamy, Saravanan; Hochgreb, Simone
2013-11-01
Coal is an important energy resource for future demand for electricity, as coal reserves are much more abundant than those of other fossil fuels. In pulverized coal fired power plants, it is very important to improve the technology for the control of environmental pollutants such as nitrogen oxide, sulfur oxide and ash particles including unburned carbon. In order to achieve these requirements, understanding the pulverized coal combustion mechanism is necessary. However, the combustion process of the pulverized coal is not well clarified so far since pulverized coal combustion is a complicated phenomenon in which the maximum flame temperature exceeds 1500 degrees Celsius and some substances which can hardly be measured, for example, radical species and highly reactive solid particles are included. Accordingly, development of new combustion furnaces and burners requires high cost and takes a long period. In this study, a large-eddy simulation (LES) is applied to a pulverized coal combustion field and the results will be compared with the experiment. The results show that present LES can capture the general feature of the pulverized coal swirl jet flame.
On the Computation of Sound by Large-Eddy Simulations
NASA Technical Reports Server (NTRS)
Piomelli, Ugo; Streett, Craig L.; Sarkar, Sutanu
1997-01-01
The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.
Large-eddy simulations of unidirectional water flow over dunes
NASA Astrophysics Data System (ADS)
Grigoriadis, D. G. E.; Balaras, E.; Dimas, A. A.
2009-06-01
The unidirectional, subcritical flow over fixed dunes is studied numerically using large-eddy simulation, while the immersed boundary method is implemented to incorporate the bed geometry. Results are presented for a typical dune shape and two Reynolds numbers, Re = 17,500 and Re = 93,500, on the basis of bulk velocity and water depth. The numerical predictions of velocity statistics at the low Reynolds number are in very good agreement with available experimental data. A primary recirculation region develops downstream of the dune crest at both Reynolds numbers, while a secondary region develops at the toe of the dune crest only for the low Reynolds number. Downstream of the reattachment point, on the dune stoss, the turbulence intensity in the developing boundary layer is weaker than in comparable equilibrium boundary layers. Coherent vortical structures are identified using the fluctuating pressure field and the second invariant of the velocity gradient tensor. Vorticity is primarily generated at the dune crest in the form of spanwise "roller" structures. Roller structures dominate the flow dynamics near the crest, and are responsible for perturbing the boundary layer downstream of the reattachment point, which leads to the formation of "horseshoe" structures. Horseshoe structures dominate the near-wall dynamics after the reattachment point, do not rise to the free surface, and are distorted by the shear layer of the next crest. The occasional interaction between roller and horseshoe structures generates tube-like "kolk" structures, which rise to the free surface and persist for a long time before attenuating.
Domain nesting for multi-scale large eddy simulation
NASA Astrophysics Data System (ADS)
Fuka, Vladimir; Xie, Zheng-Tong
2016-04-01
The need to simulate city scale areas (O(10 km)) with high resolution within street canyons in certain areas of interests necessitates different grid resolutions in different part of the simulated area. General purpose computational fluid dynamics codes typically employ unstructured refined grids while mesoscale meteorological models more often employ nesting of computational domains. ELMM is a large eddy simulation model for the atmospheric boundary layer. It employs orthogonal uniform grids and for this reason domain nesting was chosen as the approach for simulations in multiple scales. Domains are implemented as sets of MPI processes which communicate with each other as in a normal non-nested run, but also with processes from another (outer/inner) domain. It should stressed that the duration of solution of time-steps in the outer and in the inner domain must be synchronized, so that the processes do not have to wait for the completion of their boundary conditions. This can achieved by assigning an appropriate number of CPUs to each domain, and to gain high efficiency. When nesting is applied for large eddy simulation, the inner domain receives inflow boundary conditions which lack turbulent motions not represented by the outer grid. ELMM remedies this by optional adding of turbulent fluctuations to the inflow using the efficient method of Xie and Castro (2008). The spatial scale of these fluctuations is in the subgrid-scale of the outer grid and their intensity will be estimated from the subgrid turbulent kinetic energy in the outer grid.
Contrail Formation in Aircraft Wakes Using Large-Eddy Simulations
NASA Technical Reports Server (NTRS)
Paoli, R.; Helie, J.; Poinsot, T. J.; Ghosal, S.
2002-01-01
In this work we analyze the issue of the formation of condensation trails ("contrails") in the near-field of an aircraft wake. The basic configuration consists in an exhaust engine jet interacting with a wing-tip training vortex. The procedure adopted relies on a mixed Eulerian/Lagrangian two-phase flow approach; a simple micro-physics model for ice growth has been used to couple ice and vapor phases. Large eddy simulations have carried out at a realistic flight Reynolds number to evaluate the effects of turbulent mixing and wake vortex dynamics on ice-growth characteristics and vapor thermodynamic properties.
Forced turbulence in large-eddy simulation of compressible magnetohydrodynamic turbulence
Chernyshov, A. A.; Karelsky, K. V.; Petrosyan, A. S.
2010-10-15
We present the large-eddy simulation method for studying forced compressible magnetohydrodynamic turbulence. The proposed method is based on a solution of the filtered basic equations of magnetohydrodynamics by finite-difference methods and on a linear representation of the driving forces in the momentum conservation equation and the magnetic induction equation. These forces supply the production of kinetic and magnetic energies. The emphasis is placed upon the important, and not investigated, question about the ability of the large-eddy simulation approach to reproduce Kolmogorov and Iroshnikov-Kraichnan scale-invariant spectra in compressible magnetohydrodynamic flows.
Time-Domain Filtering for Spatial Large-Eddy Simulation
NASA Technical Reports Server (NTRS)
Pruett, C. David
1997-01-01
An approach to large-eddy simulation (LES) is developed whose subgrid-scale model incorporates filtering in the time domain, in contrast to conventional approaches, which exploit spatial filtering. The method is demonstrated in the simulation of a heated, compressible, axisymmetric jet, and results are compared with those obtained from fully resolved direct numerical simulation. The present approach was, in fact, motivated by the jet-flow problem and the desire to manipulate the flow by localized (point) sources for the purposes of noise suppression. Time-domain filtering appears to be more consistent with the modeling of point sources; moreover, time-domain filtering may resolve some fundamental inconsistencies associated with conventional space-filtered LES approaches.
NASA Astrophysics Data System (ADS)
Unterstrasser, S.
2014-06-01
Large-eddy simulations (LES) with Lagrangian ice microphysics were used to study the early contrail evolution during the vortex phase. Microphysical and geometrical properties of a contrail produced by a large-sized aircraft (type B777/A340) were investigated systematically for a large parameter range. Crystal loss due to adiabatic heating in the downward moving vortices was found to depend strongly on relative humidity and temperature, qualitatively similar to previous 2-D simulation results. Contrail depth is as large as 450 m for the investigated parameter range and was found to be underestimated in a previous 2-D study. Further sensitivity studies show a nonnegligible effect of the initial ice crystal size distribution and the initial ice crystal number on the crystal loss, whereas the contrail structure and ice mass evolution is only barely affected by these variations. Variation of fuel flow has the smallest effect on crystal loss. At high supersaturations, our choice of contrail spatial initialization may underestimate the ice crystal loss. The set of presented sensitivity studies is a first step toward a quantitative description of young contrails in terms of vertical extent and crystal loss. Concluding contrail-to-cirrus simulations demonstrate the relevance of vortex phase processes and its three-dimensional modeling on the later contrail-cirrus properties.
Large Eddy Simulation of Cryogenic Injection Processes at Supercritical Pressure
NASA Technical Reports Server (NTRS)
Oefelein, Joseph C.; Garcia, Roberto (Technical Monitor)
2002-01-01
This paper highlights results from the first of a series of hierarchical simulations aimed at assessing the modeling requirements for application of the large eddy simulation technique to cryogenic injection and combustion processes in liquid rocket engines. The focus is on liquid-oxygen-hydrogen coaxial injectors at a condition where the liquid-oxygen is injected at a subcritical temperature into a supercritical environment. For this situation a diffusion dominated mode of combustion occurs in the presence of exceedingly large thermophysical property gradients. Though continuous, these gradients approach the behavior of a contact discontinuity. Significant real gas effects and transport anomalies coexist locally in colder regions of the flow, with ideal gas and transport characteristics occurring within the flame zone. The current focal point is on the interfacial region between the liquid-oxygen core and the coaxial hydrogen jet where the flame anchors itself.
Implicit large eddy simulation of shock-driven material mixing.
Grinstein, F F; Gowardhan, A A; Ristorcelli, J R
2013-11-28
Under-resolved computer simulations are typically unavoidable in practical turbulent flow applications exhibiting extreme geometrical complexity and a broad range of length and time scales. An important unsettled issue is whether filtered-out and subgrid spatial scales can significantly alter the evolution of resolved larger scales of motion and practical flow integral measures. Predictability issues in implicit large eddy simulation of under-resolved mixing of material scalars driven by under-resolved velocity fields and initial conditions are discussed in the context of shock-driven turbulent mixing. The particular focus is on effects of resolved spectral content and interfacial morphology of initial conditions on transitional and late-time turbulent mixing in the fundamental planar shock-tube configuration. PMID:24146010
Large-Eddy Simulation Code Developed for Propulsion Applications
NASA Technical Reports Server (NTRS)
DeBonis, James R.
2003-01-01
A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of
Large Eddy Simulation of FDA's Idealized Medical Device.
Delorme, Yann T; Anupindi, Kameswararao; Frankel, Steven H
2013-12-01
A hybrid large eddy simulation (LES) and immersed boundary method (IBM) computational approach is used to make quantitative predictions of flow field statistics within the Food and Drug Administration's (FDA) idealized medical device. An in-house code is used, hereafter (W enoHemo(™) ), that combines high-order finite-difference schemes on structured staggered Cartesian grids with an IBM to facilitate flow over or through complex stationary or rotating geometries and employs a subgrid-scale (SGS) turbulence model that more naturally handles transitional flows [2]. Predictions of velocity and wall shear stress statistics are compared with previously published experimental measurements from Hariharan et al. [6] for the four Reynolds numbers considered. PMID:24187599
Large Eddy Simulation in a Channel with Exit Boundary Conditions
NASA Technical Reports Server (NTRS)
Cziesla, T.; Braun, H.; Biswas, G.; Mitra, N. K.
1996-01-01
The influence of the exit boundary conditions (vanishing first derivative of the velocity components and constant pressure) on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbance which is mostly damped by the grid stretching. The difference between the fully developed turbulent channel flow obtained with LES with periodicity condition and the inlet and exit and the LES with fully developed flow at the inlet and the exit boundary condition is less than 10% for equidistant grids and less than 5% for the case grid stretching. The chosen boundary condition is of interest because it may be used in complex flows with backflow at exit.
Improved engine wall models for Large Eddy Simulation (LES)
NASA Astrophysics Data System (ADS)
Plengsaard, Chalearmpol
Improved wall models for Large Eddy Simulation (LES) are presented in this research. The classical Werner-Wengle (WW) wall shear stress model is used along with near-wall sub-grid scale viscosity. A sub-grid scale turbulent kinetic energy is employed in a model for the eddy viscosity. To gain better heat flux results, a modified classical variable-density wall heat transfer model is also used. Because no experimental wall shear stress results are available in engines, the fully turbulent developed flow in a square duct is chosen to validate the new wall models. The model constants in the new wall models are set to 0.01 and 0.8, respectively and are kept constant throughout the investigation. The resulting time- and spatially-averaged velocity and temperature wall functions from the new wall models match well with the law-of-the-wall experimental data at Re = 50,000. In order to study the effect of hot air impinging walls, jet impingement on a flat plate is also tested with the new wall models. The jet Reynolds number is equal to 21,000 and a fixed jet-to-plate spacing of H/D = 2.0. As predicted by the new wall models, the time-averaged skin friction coefficient agrees well with experimental data, while the computed Nusselt number agrees fairly well when r/D > 2.0. Additionally, the model is validated using experimental data from a Caterpillar engine operated with conventional diesel combustion. Sixteen different operating engine conditions are simulated. The majority of the predicted heat flux results from each thermocouple location follow similar trends when compared with experimental data. The magnitude of peak heat fluxes as predicted by the new wall models is in the range of typical measured values in diesel combustion, while most heat flux results from previous LES wall models are over-predicted. The new wall models generate more accurate predictions and agree better with experimental data.
Film cooling from inclined cylindrical holes using large eddy simulations
NASA Astrophysics Data System (ADS)
Peet, Yulia V.
2006-12-01
The goal of the present study is to investigate numerically the physics of the flow, which occurs during the film cooling from inclined cylindrical holes, Film cooling is a technique used in gas turbine industry to reduce heat fluxes to the turbine blade surface. Large Eddy Simulation (LES) is performed modeling a realistic film cooling configuration, which consists of a large stagnation-type reservoir, feeding an array of discrete cooling holes (film holes) flowing into a flat plate turbulent boundary layer. Special computational methodology is developed for this problem, involving coupled simulations using multiple computational codes. A fully compressible LES code is used in the area above the flat plate, while a low Mach number LES code is employed in the plenum and film holes. The motivation for using different codes comes from the essential difference in the nature of the flow in these different regions. Flowfield is analyzed inside the plenum, film hole and a crossflow region. Flow inside the plenum is stagnating, except for the region close to the exit, where it accelerates rapidly to turn into the hole. The sharp radius of turning at the trailing edge of the plenum pipe connection causes the flow to separate from the downstream wall of the film hole. After coolant injection occurs, a complex flowfield is formed consisting of coherent vortical structures responsible for bringing hot crossflow fluid in contact with the walls of either the film hole or the blade, thus reducing cooling protection. Mean velocity and turbulent statistics are compared to experimental measurements, yielding good agreement for the mean flowfield and satisfactory agreement for the turbulence quantities. LES results are used to assess the applicability of basic assumptions of conventional eddy viscosity turbulence models used with Reynolds-averaged (RANS) approach, namely the isotropy of an eddy viscosity and thermal diffusivity. It is shown here that these assumptions do not hold
Large-eddy simulation of unidirectional turbulent flow over dunes
NASA Astrophysics Data System (ADS)
Omidyeganeh, Mohammad
We performed large eddy simulation of the flow over a series of two- and three-dimensional dune geometries at laboratory scale using the Lagrangian dynamic eddy-viscosity subgrid-scale model. First, we studied the flow over a standard 2D transverse dune geometry, then bedform three-dimensionality was imposed. Finally, we investigated the turbulent flow over barchan dunes. The results are validated by comparison with simulations and experiments for the 2D dune case, while the results of the 3D dunes are validated qualitatively against experiments. The flow over transverse dunes separates at the dune crest, generating a shear layer that plays a crucial role in the transport of momentum and energy, as well as the generation of coherent structures. Spanwise vortices are generated in the separated shear; as they are advected, they undergo lateral instabilities and develop into horseshoe-like structures and finally reach the surface. The ejection that occurs between the legs of the vortex creates the upwelling and downdrafting events on the free surface known as "boils". The three-dimensional separation of flow at the crestline alters the distribution of wall pressure, which may cause secondary flow across the stream. The mean flow is characterized by a pair of counter-rotating streamwise vortices, with core radii of the order of the flow depth. Staggering the crestlines alters the secondary motion; two pairs of streamwise vortices appear (a strong one, centred about the lobe, and a weaker one, coming from the previous dune, centred around the saddle). The flow over barchan dunes presents significant differences to that over transverse dunes. The flow near the bed, upstream of the dune, diverges from the centerline plane; the flow close to the centerline plane separates at the crest and reattaches on the bed. Away from the centerline plane and along the horns, flow separation occurs intermittently. The flow in the separation bubble is routed towards the horns and leaves
Large-eddy simulation of sand dune morphodynamics
NASA Astrophysics Data System (ADS)
Khosronejad, Ali; Sotiropoulos, Fotis; St. Anthony Falls Laboratory, University of Minnesota Team
2015-11-01
Sand dunes are natural features that form under complex interaction between turbulent flow and bed morphodynamics. We employ a fully-coupled 3D numerical model (Khosronejad and Sotiropoulos, 2014, Journal of Fluid Mechanics, 753:150-216) to perform high-resolution large-eddy simulations of turbulence and bed morphodynamics in a laboratory scale mobile-bed channel to investigate initiation, evolution and quasi-equilibrium of sand dunes (Venditti and Church, 2005, J. Geophysical Research, 110:F01009). We employ a curvilinear immersed boundary method along with convection-diffusion and bed-morphodynamics modules to simulate the suspended sediment and the bed-load transports respectively. The coupled simulation were carried out on a grid with more than 100 million grid nodes and simulated about 3 hours of physical time of dune evolution. The simulations provide the first complete description of sand dune formation and long-term evolution. The geometric characteristics of the simulated dunes are shown to be in excellent agreement with observed data obtained across a broad range of scales. This work was supported by NSF Grants EAR-0120914 (as part of the National Center for Earth-Surface Dynamics). Computational resources were provided by the University of Minnesota Supercomputing Institute.
Large eddy simulation and its implementation in the COMMIX code.
Sun, J.; Yu, D.-H.
1999-02-15
Large eddy simulation (LES) is a numerical simulation method for turbulent flows and is derived by spatial averaging of the Navier-Stokes equations. In contrast with the Reynolds-averaged Navier-Stokes equations (RANS) method, LES is capable of calculating transient turbulent flows with greater accuracy. Application of LES to differing flows has given very encouraging results, as reported in the literature. In recent years, a dynamic LES model that presented even better results was proposed and applied to several flows. This report reviews the LES method and its implementation in the COMMIX code, which was developed at Argonne National Laboratory. As an example of the application of LES, the flow around a square prism is simulated, and some numerical results are presented. These results include a three-dimensional simulation that uses a code developed by one of the authors at the University of Notre Dame, and a two-dimensional simulation that uses the COMMIX code. The numerical results are compared with experimental data from the literature and are found to be in very good agreement.
Towards Large Eddy Simulation of gas turbine compressors
NASA Astrophysics Data System (ADS)
McMullan, W. A.; Page, G. J.
2012-07-01
With increasing computing power, Large Eddy Simulation could be a useful simulation tool for gas turbine axial compressor design. This paper outlines a series of simulations performed on compressor geometries, ranging from a Controlled Diffusion Cascade stator blade to the periodic sector of a stage in a 3.5 stage axial compressor. The simulation results show that LES may offer advantages over traditional RANS methods when off-design conditions are considered - flow regimes where RANS models often fail to converge. The time-dependent nature of LES permits the resolution of transient flow structures, and can elucidate new mechanisms of vorticity generation on blade surfaces. It is shown that accurate LES is heavily reliant on both the near-wall mesh fidelity and the ability of the imposed inflow condition to recreate the conditions found in the reference experiment. For components embedded in a compressor this requires the generation of turbulence fluctuations at the inlet plane. A recycling method is developed that improves the quality of the flow in a single stage calculation of an axial compressor, and indicates that future developments in both the recycling technique and computing power will bring simulations of axial compressors within reach of industry in the coming years.
Inviscid Wall-Modeled Large Eddy Simulations for Improved Efficiency
NASA Astrophysics Data System (ADS)
Aikens, Kurt; Craft, Kyle; Redman, Andrew
2015-11-01
The accuracy of an inviscid flow assumption for wall-modeled large eddy simulations (LES) is examined because of its ability to reduce simulation costs. This assumption is not generally applicable for wall-bounded flows due to the high velocity gradients found near walls. In wall-modeled LES, however, neither the viscous near-wall region or the viscous length scales in the outer flow are resolved. Therefore, the viscous terms in the Navier-Stokes equations have little impact on the resolved flowfield. Zero pressure gradient flat plate boundary layer results are presented for both viscous and inviscid simulations using a wall model developed previously. The results are very similar and compare favorably to those from another wall model methodology and experimental data. Furthermore, the inviscid assumption reduces simulation costs by about 25% and 39% for supersonic and subsonic flows, respectively. Future research directions are discussed as are preliminary efforts to extend the wall model to include the effects of unresolved wall roughness. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1053575. Computational resources on TACC Stampede were provided under XSEDE allocation ENG150001.
Large Eddy Simulation of Mixing within a Hypervelocity Scramjet Combustor
NASA Astrophysics Data System (ADS)
Petty, David; Wheatley, Vincent; Pantano, Carlos; Smart, Michael
2013-11-01
The turbulent mixing of parallel hypervelocity (U = 3230 m/sec, M = 3.86) air-streams with a sonic stream of gaseous hydrogen is simulated using large eddy simulation. The resultant mixing layers are characterized by a convective Mach number of 1.20. This configuration represents parallel slot injection of hydrogen via an intrusive centerbody within a constant area rectangular combustor. A hybrid shock-capturing/zero numerical dissipation (WENO/TCD) switch method designed for simulations of compressible turbulent flows was utilized. Sub-grid scale turbulence was modeled using the stretched vortex model. Visualizations of the three dimensional turbulent structures generated behind the centerbody will be presented. It has been observed that a span-wise instability of the wake behind the centerbody is initially dominant. Further downstream, the shear-layers coalesce into a mixing wake and develop the expected large-scale coherent span-wise vortices. Ph.D. Candidate, School of Mechanical and Mining Engineering, Centre for Hypersonics.
Unsteady RANS and Large Eddy simulations of multiphase diesel injection
NASA Astrophysics Data System (ADS)
Philipp, Jenna; Green, Melissa; Akih-Kumgeh, Benjamin
2015-11-01
Unsteady Reynolds Averaged Navier-Stokes (URANS) and Large Eddy Simulations (LES) of two-phase flow and evaporation of high pressure diesel injection into a quiescent, high temperature environment is investigated. Unsteady RANS and LES are turbulent flow simulation approaches used to determine complex flow fields. The latter allows for more accurate predictions of complex phenomena such as turbulent mixing and physio-chemical processes associated with diesel combustion. In this work we investigate a high pressure diesel injection using the Euler-Lagrange method for multiphase flows as implemented in the Star-CCM+ CFD code. A dispersed liquid phase is represented by Lagrangian particles while the multi-component gas phase is solved using an Eulerian method. Results obtained from the two approaches are compared with respect to spray penetration depth and air entrainment. They are also compared with experimental data taken from the Sandia Engine Combustion Network for ``Spray A''. Characteristics of primary and secondary atomization are qualitatively evaluated for all simulation modes.
Large Eddy Simulations of Colorless Distributed Combustion Systems
NASA Astrophysics Data System (ADS)
Abdulrahman, Husam F.; Jaberi, Farhad; Gupta, Ashwani
2014-11-01
Development of efficient and low-emission colorless distributed combustion (CDC) systems for gas turbine applications require careful examination of the role of various flow and combustion parameters. Numerical simulations of CDC in a laboratory-scale combustor have been conducted to carefully examine the effects of these parameters on the CDC. The computational model is based on a hybrid modeling approach combining large eddy simulation (LES) with the filtered mass density function (FMDF) equations, solved with high order numerical methods and complex chemical kinetics. The simulated combustor operates based on the principle of high temperature air combustion (HiTAC) and has shown to significantly reduce the NOx, and CO emissions while improving the reaction pattern factor and stability without using any flame stabilizer and with low pressure drop and noise. The focus of the current work is to investigate the mixing of air and hydrocarbon fuels and the non-premixed and premixed reactions within the combustor by the LES/FMDF with the reduced chemical kinetic mechanisms for the same flow conditions and configurations investigated experimentally. The main goal is to develop better CDC with higher mixing and efficiency, ultra-low emission levels and optimum residence time. The computational results establish the consistency and the reliability of LES/FMDF and its Lagrangian-Eulerian numerical methodology.
Large Eddy Simulation of Crashback in Marine Propulsors
NASA Astrophysics Data System (ADS)
Jang, Hyunchul
Crashback is an operating condition to quickly stop a propelled vehicle, where the propeller is rotated in the reverse direction to yield negative thrust. The crashback condition is dominated by the interaction of the free stream flow with the strong reverse flow. This interaction forms a highly unsteady vortex ring, which is a very prominent feature of crashback. Crashback causes highly unsteady loads and flow separation on the blade surface. The unsteady loads can cause propulsor blade damage, and also affect vehicle maneuverability. Crashback is therefore well known as one of the most challenging propeller states to analyze. This dissertation uses Large-Eddy Simulation (LES) to predict the highly unsteady flow field in crashback. A non-dissipative and robust finite volume method developed by Mahesh et al. (2004) for unstructured grids is applied to flow around marine propulsors. The LES equations are written in a rotating frame of reference. The objectives of this dissertation are: (1) to understand the flow physics of crashback in marine propulsors with and without a duct, (2) to develop a finite volume method for highly skewed meshes which usually occur in complex propulsor geometries, and (3) to develop a sliding interface method for simulations of rotor-stator propulsor on parallel platforms. LES is performed for an open propulsor in crashback and validated against experiments performed by Jessup et al. (2004). The LES results show good agreement with experiments. Effective pressures for thrust and side-force are introduced to more clearly understand the physical sources of thrust and side-force. Both thrust and side-force are seen to be mainly generated from the leading edge of the suction side of the propeller. This implies that thrust and side-force have the same source---the highly unsteady leading edge separation. Conditional averaging is performed to obtain quantitative information about the complex flow physics of high- or low-amplitude events. The
Possible modifications to implicit large-eddy simulation
NASA Astrophysics Data System (ADS)
McDonough, J. M.
2009-11-01
Implicit large-eddy simulation (ILES) provides an advantage over more usual LES approaches in that its construction does not involve filtering of the governing equations and, as a consequence, removal of the need to develop sub-grid scale (SGS) models to represent artificial stresses arising from this filtering. At the same time, it is clear that ILES is simply an under-resolved direct numerical simulation with advanced treatments of advection terms to better control numerical stability via dissipation that otherwise would have been provided by a SGS model. As such it cannot be expected to accurately predict interactions of fluid turbulence with other physical phenomena (e.g., heat and mass transfer, chemical kinetics) on subgrid scales---as is also true of usual forms of LES. In this talk we describe a straightforward technique, based on formal multi-scale methods, whereby SGS interactions can be introduced to enhance resolved-scale results computed as in ILES, and we discuss derivation of a class of efficient models based on the ``poor man's Navier--Stokes equation'' (McDonough, Phys.Rev. E 79, 2009; McDonough and Huang, Int.J.Numer. Meth. Fluids 44, 2004). Properties of these models will be presented for a moderate-Re 3-D lid-driven cavity problem.
Large-eddy simulation of density currents on inclined beds
NASA Astrophysics Data System (ADS)
Chawdhary, Saurabh; Khosronejad, Ali; Christodoulou, George; Sotiropoulos, Fotis
2013-11-01
Density currents are stratified flow in presence of density differential and gravity field. We carry out Large-Eddy Simulation (LES) to simulate the flow of a density current formed over sloped bed due to an incoming jet of heavy density salty water for two different cases of bed slope: (a) 5 degrees and (b) 15 degrees. The Reynolds and Richardson numbers based on inlet height and inlet velocity were (a) 1100 and 0.471, and (b) 2000 and 0.0355, respectively. The Schmidt number is set equal to 620, which corresponds to the value for salt-water. The computed results are compared with laboratory experiments in terms of overall shape of the heavy-density plume and its spreading rate and are shown to be in reasonable agreement. The instantaneous LES flow fields are further analyzed to gain novel insights into the rich dynamics of coherent vortical structures in the flow. The half-width of the plume is plotted as a function of downstream length and found to exhibit three different regions on a log scale, in agreement with previous experimental findings. We acknowledge computational support from the Minnesota Supercomputing Institute.
Large eddy simulation of incompressible turbulent channel flow
NASA Technical Reports Server (NTRS)
Moin, P.; Reynolds, W. C.; Ferziger, J. H.
1978-01-01
The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.
NASA Astrophysics Data System (ADS)
Berner, A. H.; Bretherton, C. S.; Wood, R.
2015-05-01
For the first time, a large eddy simulation (LES) coupled to a bulk aerosol scheme is used to simulate an aircraft-sampled ship track. The track was formed by the M/V Sanko Peace on 13 June 1994 in a shallow drizzling boundary layer with high winds but very low background aerosol concentrations (10 cm-3). A Lagrangian framework is used to simulate the evolution of a short segment of track as it is advected away from the ship for 8 h (a downwind distance exceeding 570 km). Using aircraft observations for initialization, good agreement is obtained between the simulated and observed features of the ambient boundary layer outside the track, including the organization of the cloud into mesoscale rolls. After 8 h, a line of aerosol is injected to start the ship track. The simulation successfully reproduces the significant albedo enhancement and suppression of drizzle observed within the track. The aerosol concentration within the track dilutes as it broadens due to turbulent mixing. A sensitivity study shows the broadening rate strongly depends on the alignment between the track and the wind-aligned boundary layer rolls, as satellite images of ship tracks suggest. Entrainment is enhanced within the simulated track, but the observed 100 m elevation of the ship track above the surrounding layer is not simulated, possibly because the LES quickly sharpens the rather weak observed inversion. Liquid water path within the simulated track increases with time even as the ambient liquid water path is decreasing. The albedo increase in the track from liquid water and cloud fraction enhancement (second indirect effect) eventually exceeds that from cloud droplet number increases (first indirect or Twomey effect). In a sensitivity study with a higher initial ambient aerosol concentration, stronger ship track aerosol source, and much weaker drizzle, there is less liquid water inside the track than outside for several hours downwind, consistent with satellite estimates for such
NASA Astrophysics Data System (ADS)
Berner, A. H.; Bretherton, C. S.; Wood, R.
2014-09-01
For the first time, a large-eddy simulation (LES) coupled to a bulk aerosol scheme is used to simulate an aircraft-sampled ship track. The track was formed by the M/V Sanko Peace on 13 June 1994 in a shallow drizzling boundary layer with high winds but very low background aerosol concentrations (10 cm-3). A Lagrangian framework is used to simulate the evolution of a short segment of track as it is advected away from the ship for eight hours (a downwind distance exceeding 570 km). Using aircraft observations for initialization, good agreement is obtained between the simulated and observed features of the ambient boundary layer outside the track, including the organization of cloud into mesoscale rolls. After eight hours, a line of aerosol is injected to start the ship track. The simulation successfully reproduces the significant albedo enhancement and suppression of drizzle observed within the track. The aerosol concentration within the track dilutes as it broadens due to turbulent mixing. A sensitivity study shows the broadening rate strongly depends on the alignment between the track and the wind-aligned boundary layer rolls, as satellite images of ship tracks suggest. Entrainment is enhanced within the simulated track, but the observed 100 m elevation of the ship track above the surrounding layer is not simulated, possibly because the LES quickly sharpens the rather weak observed inversion. Liquid water path within the simulated track increases with time even as the ambient liquid water path is decreasing. The albedo increase in the track from liquid water and cloud fraction enhancement (second indirect effect) eventually exceeds that from cloud droplet number increases (first indirect or Twomey effect). In a sensitivity study with a higher initial ambient aerosol concentration, stronger ship track aerosol source, and much weaker drizzle, there is less liquid water inside the track than outside for several hours downwind, consistent with satellite estimates for
Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows
NASA Technical Reports Server (NTRS)
Givi, Peyman; Madnia, Cyrus K.; Steinberger, Craig J.
1990-01-01
This research is involved with the implementation of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program to extend the present capabilities of this method was initiated for the treatment of chemically reacting flows. In the DNS efforts, the focus is on detailed investigations of the effects of compressibility, heat release, and non-equilibrium kinetics modelings in high speed reacting flows. Emphasis was on the simulations of simple flows, namely homogeneous compressible flows, and temporally developing high speed mixing layers.
NASA Astrophysics Data System (ADS)
Juhui, Chen; Yanjia, Tang; Dan, Li; Pengfei, Xu; Huilin, Lu
2013-07-01
Flow behavior of gas and particles is predicted by the large eddy simulation of gas-second order moment of solid model (LES-SOM model) in the simulation of flow behavior in CFB. This study shows that the simulated solid volume fractions along height using a two-dimensional model are in agreement with experiments. The velocity, volume fraction and second-order moments of particles are computed. The second-order moments of clusters are calculated. The solid volume fraction, velocity and second order moments are compared at the three different model constants.
Minimum-dissipation models for large-eddy simulation
NASA Astrophysics Data System (ADS)
Rozema, Wybe; Bae, Hyun J.; Moin, Parviz; Verstappen, Roel
2015-08-01
Minimum-dissipation eddy-viscosity models are a class of sub-filter models for large-eddy simulation that give the minimum eddy dissipation required to dissipate the energy of sub-filter scales. A previously derived minimum-dissipation model is the QR model. This model is based on the invariants of the resolved rate-of-strain tensor and has many desirable properties. It appropriately switches off for laminar and transitional flows, has low computational complexity, and is consistent with the exact sub-filter tensor on isotropic grids. However, the QR model proposed in the literature gives insufficient eddy dissipation. It is demonstrated that this can be corrected by increasing the model constant. The corrected QR model gives good results in simulations of decaying grid turbulence on an isotropic grid. On anisotropic grids the QR model is not consistent with the exact sub-filter tensor and requires an approximation of the filter width. It is demonstrated that the results of the QR model on anisotropic grids are primarily determined by the used filter width approximation, and that no approximation gives satisfactory results in simulations of both a temporal mixing layer and turbulent channel flow. A new minimum-dissipation model for anisotropic grids is proposed. This anisotropic minimum-dissipation (AMD) model generalizes the desirable practical and theoretical properties of the QR model to anisotropic grids and does not require an approximation of the filter width. The AMD model is successfully applied in simulations of decaying grid turbulence on an isotropic grid and in simulations of a temporal mixing layer and turbulent channel flow on anisotropic grids.
Large Eddy Simulation and the Filtered Probability Density Function Method
NASA Astrophysics Data System (ADS)
Jones, W. P.; Navarro-Martinez, S.
2009-12-01
Recently there is has been increased interest in modelling combustion processes with high-levels of extinction and re-ignition. Such system often lie beyond the scope of conventional single scalar-based models. Large Eddy Simulation (LES) has shown a large potential for describing turbulent reactive systems, though combustion occurs at the smallest unresolved scales of the flow and must be modelled. In the sub-grid Probability Density Function (pdf) method approximations are devised to close the evolution equation for the joint-pdf which is then solved directly. The paper describes such an approach and concerns, in particular, the Eulerian stochastic field method of solving the pdf equation. The paper examines the capabilities of the LES-pdf method in capturing auto-ignition and extinction events in different partially premixed configurations with different fuels (hydrogen, methane and n-heptane). The results show that the LES-pdf formulation can capture different regimes without any parameter adjustments, independent of Reynolds numbers and fuel type.
Large eddy simulation modelling of combustion for propulsion applications.
Fureby, C
2009-07-28
Predictive modelling of turbulent combustion is important for the development of air-breathing engines, internal combustion engines, furnaces and for power generation. Significant advances in modelling non-reactive turbulent flows are now possible with the development of large eddy simulation (LES), in which the large energetic scales of the flow are resolved on the grid while modelling the effects of the small scales. Here, we discuss the use of combustion LES in predictive modelling of propulsion applications such as gas turbine, ramjet and scramjet engines. The LES models used are described in some detail and are validated against laboratory data-of which results from two cases are presented. These validated LES models are then applied to an annular multi-burner gas turbine combustor and a simplified scramjet combustor, for which some additional experimental data are available. For these cases, good agreement with the available reference data is obtained, and the LES predictions are used to elucidate the flow physics in such devices to further enhance our knowledge of these propulsion systems. Particular attention is focused on the influence of the combustion chemistry, turbulence-chemistry interaction, self-ignition, flame holding burner-to-burner interactions and combustion oscillations. PMID:19531515
A family of dynamic models for large-eddy simulation
NASA Technical Reports Server (NTRS)
Carati, D.; Jansen, K.; Lund, T.
1995-01-01
Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.
Large eddy simulation predictions of absolutely unstable round hot jet
NASA Astrophysics Data System (ADS)
Boguslawski, A.; Tyliszczak, A.; Wawrzak, K.
2016-02-01
The paper presents a novel view on the absolute instability phenomenon in heated variable density round jets. As known from literature the global instability mechanism in low density jets is released when the density ratio is lower than a certain critical value. The existence of the global modes was confirmed by an experimental evidence in both hot and air-helium jets. However, some differences in both globally unstable flows were observed concerning, among others, a level of the critical density ratio. The research is performed using the Large Eddy Simulation (LES) method with a high-order numerical code. An analysis of the LES results revealed that the inlet conditions for the velocity and density distributions at the nozzle exit influence significantly the critical density ratio and the global mode frequency. Two inlet velocity profiles were analyzed, i.e., the hyperbolic tangent and the Blasius profiles. It was shown that using the Blasius velocity profile and the uniform density distribution led to a significantly better agreement with the universal scaling law for global mode frequency.
NASA Astrophysics Data System (ADS)
Stoll, J. R.; Shingleton, N. D.; Bosveld, F.
2010-12-01
Accurately reproducing the dynamic two-way interaction between the land surface and the atmosphere in the stable boundary layer (SBL) requires detailed treatment of the governing physical processes. Increasingly, large-eddy simulation (LES) is used for this purpose. In many studies, the dominant treatment of surface boundary conditions is to specify a known state or flux. This results in one-way or weak two-way coupling between the land surface and the boundary layer. The impact of how this coupling is modeled on atmospheric boundary layer (ABL) dynamics is still not fully understood, especially under transitional and weakly turbulent conditions. Here, LES that is fully coupled to a land-surface model (LSM) is used to investigate the nocturnal and the transitional periods of the diurnal cycle. The LSM explicitly solves for the transport of heat and water in a one-dimensional column of the upper soil and is coupled to the atmosphere through a surface energy budget. The fully coupled LES-LSM is used to simulate the third GEWEX (Global Energy and Water Cycle Experiment) ABL (GABLS3) LES intercomparison case. Turbulent boundary layer profiles and surface fluxes are compared to field data and results from simulations using three different levels of physical description as lower boundary conditions. These include simulations with prescribed temperature and moisture state, with a LSM that uses a bare-soil approximation and a LSM that include a skin layer. Overall, simulations with all three types of boundary conditions compare fairly well with the general trends observed in the field data for surface fluxes and boundary layer turbulence statistical profiles during the intercomparison time period (night to early morning) with a few differences. The LES-LSM model under-predicts the latent heat flux during the night and over-predicts the ground heat and moisture fluxes. The addition of a skin layer improves flux predictions during the night and early morning. Surface fluxes
Large eddy simulation of a lifted turbulent jet flame
Ferraris, S.A.; Wen, J.X.
2007-09-15
The flame index concept for large eddy simulation developed by Domingo et al. [P. Domingo, L. Vervisch, K. Bray, Combust. Theory Modell. 6 (2002) 529-551] is used to capture the partially premixed structure at the leading point and the dual combustion regimes further downstream on a turbulent lifted flame, which is composed of premixed and nonpremixed flame elements each separately described under a flamelet assumption. Predictions for the lifted methane/air jet flame experimentally tested by Mansour [M.S. Mansour, Combust. Flame 133 (2003) 263-274] are made. The simulation covers a wide domain from the jet exit to the far flow field. Good agreement with the data for the lift-off height and the mean mixture fraction has been achieved. The model has also captured the double flames, showing a configuration similar to that of the experiment which involves a rich premixed branch at the jet center and a diffusion branch in the outer region which meet at the so-called triple point at the flame base. This basic structure is contorted by eddies coming from the jet exit but remains stable at the lift-off height. No lean premixed branches are observed in the simulation or and experiment. Further analysis on the stabilization mechanism was conducted. A distinction between the leading point (the most upstream point of the flame) and the stabilization point was made. The later was identified as the position with the maximum premixed heat release. This is in line with the stabilization mechanism proposed by Upatnieks et al. [A. Upatnieks, J. Driscoll, C. Rasmussen, S. Ceccio, Combust. Flame 138 (2004) 259-272]. (author)
NASA Astrophysics Data System (ADS)
Cox, Christopher; Liang, Chunlei
2011-11-01
In this investigation, we implement a high-order three-dimensional spectral difference (SD) method to solve the compressible Navier-Stokes equations on an unstructured moving deformable grid. Presently, the SD method is used to perform simulations of compressible flow past an oscillating circular cylinder. Oscillations parallel and normal to the free stream are considered at a fixed Reynolds number of 4000, oscillation frequency of 1 Hz , and oscillation amplitude of 20% cylinder diameter. We extend this study to large eddy simulations with the integration of a Smagorinsky-type subgrid-scale (SGS) model. Computational results will be compared to experimental data. The effectiveness of the large eddy simulation in capturing the vortex dynamics in the wake is analyzed. This work is funded by the Mechanical & Aerospace Engineering Department at George Washington University
NASA Technical Reports Server (NTRS)
Zhang, Zhibo; Ackerman, Andrew S.; Feingold, Graham; Platnick, Steven; Pincus, Robert; Xue, Huiwen
2012-01-01
This study investigates effects of drizzle and cloud horizontal inhomogeneity on cloud effective radius (re) retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS). In order to identify the relative importance of various factors, we developed a MODIS cloud property retrieval simulator based on the combination of large-eddy simulations (LES) and radiative transfer computations. The case studies based on synthetic LES cloud fields indicate that at high spatial resolution (100 m) 3-D radiative transfer effects, such as illumination and shadowing, can induce significant differences between retrievals ofre based on reflectance at 2.1 m (re,2.1) and 3.7 m (re,3.7). It is also found that 3-D effects tend to have stronger impact onre,2.1 than re,3.7, leading to positive difference between the two (re,3.72.1) from illumination and negative re,3.72.1from shadowing. The cancellation of opposing 3-D effects leads to overall reasonable agreement betweenre,2.1 and re,3.7 at high spatial resolution as far as domain averages are concerned. At resolutions similar to MODIS, however, re,2.1 is systematically larger than re,3.7when averaged over the LES domain, with the difference exhibiting a threshold-like dependence on bothre,2.1and an index of the sub-pixel variability in reflectance (H), consistent with MODIS observations. In the LES cases studied, drizzle does not strongly impact reretrievals at either wavelength. It is also found that opposing 3-D radiative transfer effects partly cancel each other when cloud reflectance is aggregated from high spatial resolution to MODIS resolution, resulting in a weaker net impact of 3-D radiative effects onre retrievals. The large difference at MODIS resolution between re,3.7 and re,2.1 for highly inhomogeneous pixels with H 0.4 can be largely attributed to what we refer to as the plane-parallelrebias, which is attributable to the impact of sub-pixel level horizontal variability of cloud optical thickness onre retrievals
Thermobaric cabbeling over Maud Rise: Theory and large eddy simulation
NASA Astrophysics Data System (ADS)
Harcourt, Ramsey R.
2005-10-01
A Large Eddy Simulation (LES) of the wintertime upper ocean below seasonal Antarctic ice cover over Maud Rise was carried out using observed time-dependent surface forcing from 1994 Antarctic Zone Flux Experiment (ANZFLUX) observations. Surface ice formation increases the density of the cold, fresher Surface Mixed Layer (SML), that overlies warmer, saltier Weddell Deep Water (WDW). This reduces the stability of the thermocline until it reaches a critical point for instabilities arising from the nonlinear equation of state (NES) for seawater density ρ. This simulation was intended to model the thermobaric detrainment of SML fluid, a NES instability predicted to result from the dependence of seawater density on the product θP of temperature and pressure. Instead, model results demonstrate a different instability arising from the combination of thermobaricity with cabbeling, the NES effect due primarily to the dependence of ρ on θ2. This combined thermobaric cabbeling instability drives turbulent convection in a deep interior mixed layer (IML) that may grow hundreds of meters thick below the thermocline, largely decoupled from SML dynamics. In the LES, thermobaric cabbeling and IML convection shoals the SML through entrainment from below until ice motion increases in the observationally-based model forcing. Increased upper ocean model heat flux due to higher ice speed melts surface ice, increasing thermocline stratification and eventually bringing the simulated instability to a halt. In an auxiliary simulation the lull preceding strong ice motion in field observations is artificially extended by temporarily holding model surface forcing constant until the SML shoals entirely, bringing the modified WDW of the IML, 2 °C above freezing, directly to the surface. Subsequently, reverting to the observed surface forcing and its attendant strong ice motion melts the ice cover entirely, demonstrating a possible mechanism for open ocean Antarctic polynya formation. The
Structural subgrid-scale modeling for large-eddy simulation: A review
NASA Astrophysics Data System (ADS)
Lu, Hao; Rutland, Christopher J.
2016-03-01
Accurately modeling nonlinear interactions in turbulence is one of the key challenges for large-eddy simulation (LES) of turbulence. In this article, we review recent studies on structural subgrid scale modeling, focusing on evaluating how well these models predict the effects of small scales. The article discusses a priori and a posteriori test results. Other nonlinear models are briefly discussed, and future prospects are noted.
Structural subgrid-scale modeling for large-eddy simulation: A review
NASA Astrophysics Data System (ADS)
Lu, Hao; Rutland, Christopher J.
2016-08-01
Accurately modeling nonlinear interactions in turbulence is one of the key challenges for large-eddy simulation (LES) of turbulence. In this article, we review recent studies on structural subgrid scale modeling, focusing on evaluating how well these models predict the effects of small scales. The article discusses a priori and a posteriori test results. Other nonlinear models are briefly discussed, and future prospects are noted.
Prediction and control of turbulent aero-optical distortion using large eddy simulation
NASA Astrophysics Data System (ADS)
Childs, Robert E.
1993-06-01
The problem of aero-optical distortion caused by turbulence in high speed mixing layers was studied using large eddy simulation (LES) as the model of turbulence. The accuracy of LES is established for global features of the mixing layer, such as mean growth rate and statistics of turbulent velocity fluctuations. LES was then used to assess two concepts for suppressing density fluctuations and aero-optical distortion, lateral convergence and streamline curvature, and one of these was found to be reasonably effective.
Large eddy simulation of mixing between hot and cold sodium flows - comparison with experiments
Simoneau, J.P.; Noe, H.; Menant, B.
1995-09-01
The large eddy simulation is becoming a potential powerful tool for the calculation of turbulent flows. In nuclear liquid metal cooled fast reactors, the knowledge of the turbulence characteristics is of great interest for the prediction and the analysis of thermal stripping phenomena. The objective of this paper is to give a contribution in the evaluation of the large eddy simulation technique is an individual case. The problem chosen is the case of the mixing between hot and cold sodium flows. The computations are compared with available sodium tests. This study shows acceptable qualitative results but the simple model used is not able to predict the turbulence characteristics. More complex models including larger domains around the fluctuating zone and fluctuating boundary conditions could be necessary. Validation works are continuing.
NASA Technical Reports Server (NTRS)
Bardino, J.; Ferziger, J. H.; Reynolds, W. C.
1983-01-01
The physical bases of large eddy simulation and subgrid modeling are studied. A subgrid scale similarity model is developed that can account for system rotation. Large eddy simulations of homogeneous shear flows with system rotation were carried out. Apparently contradictory experimental results were explained. The main effect of rotation is to increase the transverse length scales in the rotation direction, and thereby decrease the rates of dissipation. Experimental results are shown to be affected by conditions at the turbulence producing grid, which make the initial states a function of the rotation rate. A two equation model is proposed that accounts for effects of rotation and shows good agreement with experimental results. In addition, a Reynolds stress model is developed that represents the turbulence structure of homogeneous shear flows very well and can account also for the effects of system rotation.
Numerical simulation of a closed rotor-stator system using Large Eddy Simulation
NASA Astrophysics Data System (ADS)
Amouyal, Solal Abraham Teva
A large eddy simulation of an enclosed annular rotor stator cavity is presented. The geometry is characterized by a large aspect ratio G = (b-a)/h = 18.32 and a small radius ratio a/b = 0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate o under consideration is equivalent to the rotational Reynolds number Re = o b2 /nu= 9.5x104 , where nu is the kinematic viscosity. The main objective of this study is to correctly simulate the rotor stator cavity using a low order numerical scheme on unstructured grids. The numerical simulations were run on the software AVBP developed by the Centre Europeen de Recherche et de Formation Avancee en Calcul Scientific. The results were compared to the experimental results obtained by Sebastien Poncet of Universit e Aix-Marseille. Two large eddy simulations techniques were used: the Smagorinsky and Wall-adapting local eddy-viscosity models. The simulations were run on three set of grids, each with a different cell resolution-14, 35 and 50- along the thickness of the system. Results from each mesh show a good qualitative agreement of the mean velocity field with Poncet's experimental results. It was found that the Samgorinsky model to be more appropriate for this configuration.
Calibration of a new very large eddy simulation (VLES) methodology for turbulent flow simulation
NASA Astrophysics Data System (ADS)
Han, XingSi; Ye, TaoHong; Chen, YiLiang
2012-10-01
Following the idea of Speziale's Very Large Eddy Simulation (VLES) method, a new unified hybrid simulation approach was proposed which can change seamlessly from RANS (Reynolds-Averaged Navier-Stokes) to LES (Large Eddy Simulation) method depending on the numerical resolution. The model constants were calibrated in accordance with other hybrid methods. Besides being able to approach the two limits of RANS and LES, the new model also provides a proper VLES mode between the two limits, and thus can be used for a wide range of mesh resolutions. Also RANS simulation can be recovered near the wall which is similar to the Detached Eddy Simulation (DES) concept. This new methodology was implemented into Wilcox's k-ω model and applications were conducted for fully developed turbulent channel flow at Re τ = 395 and turbulent flow past a square cylinder at Re = 22000. Results were compared with LES predictions and other studies. The new method is found to be quite efficient in resolving large flow structures, and can predict satisfactory results on relative coarse mesh.
Large Eddy Simulation of stratified flows over structures
NASA Astrophysics Data System (ADS)
Fuka, V.; Brechler, J.
2013-04-01
We tested the ability of the LES model CLMM (Charles University Large-Eddy Microscale Model) to model the stratified flow around three dimensional hills. We compared the quantities, as the height of the dividing streamline, recirculation zone length or length of the lee waves with experiments by Hunt and Snyder[3] and numerical computations by Ding, Calhoun and Street[5]. The results mostly agreed with the references, but some important differences are present.
NASA Astrophysics Data System (ADS)
Watanabe, Tomoaki; Sakai, Yasuhiko; Nagata, Koji; Ito, Yasumasa
2016-04-01
Spatially developing planar jets with passive scalar transports are simulated for various Reynolds (Re = 2200, 7000, and 22 000) and Schmidt numbers (Sc = 1, 4, 16, 64, and 128) by the implicit large eddy simulation (ILES) using low-pass filtering as an implicit subgrid-scale model. The budgets of resolved turbulent kinetic energy k and scalar variance < {φ }\\prime 2> are explicitly evaluated from the ILES data except for the dissipation terms, which are obtained from the balance in the transport equations. The budgets of k and < {φ }\\prime 2> in the ILES agree well with the DNS and experiments for both high and low Re cases. The streamwise decay of the mean turbulent kinetic energy dissipation rate obeys the power low obtained by the scaling argument. The mechanical-to-scalar timescale ratio C ϕ is evaluated in the self-similar region. For the high Re case, C ϕ is close to the isotropic value (C ϕ = 2) near the jet centerline. However, when Re is not large, C ϕ is smaller than 2 and depends on the Schmidt number. The T/NT interface is also investigated by using the scalar isosurface. The velocity and scalar fields near the interface depend on the interface orientation for all Re. The velocity toward the interface is observed near the interface facing in the streamwise, cross-streamwise, and spanwise directions in the planar jet in the resolved velocity field.
NASA Astrophysics Data System (ADS)
Heinze, Rieke; Mironov, Dmitrii; Raasch, Siegfried
2016-03-01
A detailed analysis of the pressure-scrambling terms (i.e., the pressure-strain and pressure gradient-scalar covariances) in the Reynolds-stress and scalar-flux budgets for cloud-topped boundary layers (CTBLs) is performed using high-resolution large-eddy simulation (LES). Two CTBLs are simulated — one with trade wind shallow cumuli, and the other with nocturnal marine stratocumuli. The pressure-scrambling terms are decomposed into contributions due to turbulence-turbulence interactions, mean velocity shear, buoyancy, and Coriolis effects. Commonly used models of these contributions, including a simple linear model most often used in geophysical applications and a more sophisticated two-component-limit (TCL) nonlinear model, are tested against the LES data. The decomposition of the pressure-scrambling terms shows that the turbulence-turbulence and buoyancy contributions are most significant for cloud-topped boundary layers. The Coriolis contribution is negligible. The shear contribution is generally of minor importance inside the cloudy layers, but it is the leading-order contribution near the surface. A comparison of models of the pressure-scrambling terms with the LES data suggests that the more complex TCL model is superior to the simple linear model only for a few contributions. The linear model is able to reproduce the principal features of the pressure-scrambling terms reasonably well. It can be applied in the second-order turbulence modeling of cloud-topped boundary layer flows, provided some uncertainties are tolerated.
Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows
NASA Technical Reports Server (NTRS)
Givi, Peyman; Madnia, C. K.; Steinberger, C. J.; Tsai, A.
1991-01-01
This research is involved with the implementations of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program was initiated to extend the present capabilities of this method for the treatment of chemically reacting flows, whereas in the DNS efforts, focus was on detailed investigations of the effects of compressibility, heat release, and nonequilibrium kinetics modeling in high speed reacting flows. The efforts to date were primarily focussed on simulations of simple flows, namely, homogeneous compressible flows and temporally developing hign speed mixing layers. A summary of the accomplishments is provided.
Large eddy simulation of vortex shedding and pressure fluctuation in aerostatic bearings
NASA Astrophysics Data System (ADS)
Zhu, Jincheng; Chen, Han; Chen, Xuedong
2013-07-01
In aerostatic bearings, high speed air flow may induce small vibration, which has been harmful to the improvement of moving and positioning accuracy of aerostatically supported devices in ultra-precision applications. In this paper, the transient flow field in the aerostatic bearing is numerically investigated using the large eddy simulation method. Turbulent structures are studied and vortex shedding phenomenon is discovered in the bearing recess. Our computational results demonstrate that vortex shedding causes pressure fluctuation in the bearing clearance. Relationship between pressure fluctuation and bearing vibration is established based on our simulation results and experimentally measured vibration strength.
Hybrid Reynolds-Averaged/Large-Eddy Simulations of a Coaxial Supersonic Free-Jet Experiment
NASA Technical Reports Server (NTRS)
Baurle, Robert A.; Edwards, Jack R.
2010-01-01
Reynolds-averaged and hybrid Reynolds-averaged/large-eddy simulations have been applied to a supersonic coaxial jet flow experiment. The experiment was designed to study compressible mixing flow phenomenon under conditions that are representative of those encountered in scramjet combustors. The experiment utilized either helium or argon as the inner jet nozzle fluid, and the outer jet nozzle fluid consisted of laboratory air. The inner and outer nozzles were designed and operated to produce nearly pressure-matched Mach 1.8 flow conditions at the jet exit. The purpose of the computational effort was to assess the state-of-the-art for each modeling approach, and to use the hybrid Reynolds-averaged/large-eddy simulations to gather insight into the deficiencies of the Reynolds-averaged closure models. The Reynolds-averaged simulations displayed a strong sensitivity to choice of turbulent Schmidt number. The initial value chosen for this parameter resulted in an over-prediction of the mixing layer spreading rate for the helium case, but the opposite trend was observed when argon was used as the injectant. A larger turbulent Schmidt number greatly improved the comparison of the results with measurements for the helium simulations, but variations in the Schmidt number did not improve the argon comparisons. The hybrid Reynolds-averaged/large-eddy simulations also over-predicted the mixing layer spreading rate for the helium case, while under-predicting the rate of mixing when argon was used as the injectant. The primary reason conjectured for the discrepancy between the hybrid simulation results and the measurements centered around issues related to the transition from a Reynolds-averaged state to one with resolved turbulent content. Improvements to the inflow conditions were suggested as a remedy to this dilemma. Second-order turbulence statistics were also compared to their modeled Reynolds-averaged counterparts to evaluate the effectiveness of common turbulence closure
Coelho, P.J.
2009-05-15
An analysis of the relevance of turbulence-radiation interaction in the numerical simulation of turbulent reactive flows is presented. A semi-causal stochastic model was used to generate a time-series of turbulent scalar fluctuations along optical paths of Sandia flame D, a widely studied piloted turbulent jet nonpremixed flame. The radiative transfer equation was integrated along these paths for every realization using a grid resolution typical of a direct numerical simulation. The correlated k-distribution method was employed to compute the radiative properties of the medium. The results were used to determine the ensemble average, as well as the extreme values, of quantities that indicate the importance of the turbulence-radiation interaction. Several approximate methods are then proposed to solve the filtered radiative transfer equation in the framework of large eddy simulations. The proposed methods are applicable along with combustion models that either assume the filtered probability density function of a conserved scalar or solve a transport equation for a joint scalar or joint scalar/velocity filtered density function. It is concluded that the errors resulting from neglecting the turbulence-radiation interaction in large eddy simulations are much lower than those found in Reynolds-averaged Navier-Stokes calculations. The optically thin fluctuation approximation may be extended to large eddy simulations yielding predictions in excellent agreement with the reference solution. If the turbulence-radiation interaction is accounted for using this approximation, the average relative error of the filtered total radiation intensity is generally below 0.3% for the studied flame. (author)
Large eddy simulation of Rayleigh-Taylor instability using the arbitrary Lagrangian-Eulerian method
Darlington, R
1999-12-01
This research addresses the application of a large eddy simulation (LES) to Arbitrary Lagrangian Eulerian (ALE) simulations of Rayleigh-Taylor instability. First, ALE simulations of simplified Rayleigh-Taylor instability are studied. The advantages of ALE over Eulerian simulations are shown. Next, the behavior of the LES is examined in a more complicated ALE simulation of Rayleigh-Taylor instability. The effects of eddy viscosity and stochastic backscatter are examined. The LES is also coupled with ALE to increase grid resolution in areas where it is needed. Finally, the methods studied above are applied to two sets of experimental simulations. In these simulations, ALE allows the mesh to follow expanding experimental targets, while LES can be used to mimic the effect of unresolved instability modes.
Large-Eddy Simulation of Flow Around a NACA 4412 Airfoil Using Unstructured Grids
NASA Technical Reports Server (NTRS)
Jansen, Kenneth
1996-01-01
Large-eddy simulation (LES) has matured to the point where application to complex flows is described. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential.
Reynolds-constrained large-eddy simulation of compressible flow over a compression ramp
NASA Astrophysics Data System (ADS)
Xiao, Zuoli; Chen, Liang
2015-11-01
A novel large-eddy simulation (LES) method is introduced for numerical simulation of wall-bounded compressible turbulent flows. The subgrid-scale (SGS) model in this method is designed to be composed of two parts depending on the distance to the nearest wall. In the near-wall region, both the mean SGS stress and heat flux are constrained by external Reynolds stress and heat flux to ensure the total target quantities, while the fluctuating SGS stress and heat flux are closed in a traditional fashion but using residual model parameterizations. In the far-wall region, the conventional SGS model is directly employed with necessary smoothing operation in the neighborhood of the constrained-unconstrained interface, which might be different for the stress and heat flux depending on the flow configuration. Compressible flow over a compression ramp is numerically studied using the new LES technique. The results are compared with the available experimental and direct numerical simulation (DNS) data, and those from traditional LES and detached-eddy simulation (DES). It turns out that the Reynolds-constrained large-eddy simulation (RCLES) method can predict the size of the separation bubble, mean flow profile, and friction force, etc. more accurately than traditional LES and DES techniques. Moreover, the RCLES method proves to be much less sensitive to the grid resolution than traditional LES method, and makes pure LES of flows of engineering interest feasible with moderate grids.
Large Eddy Simulations of Severe Convection Induced Turbulence
NASA Technical Reports Server (NTRS)
Ahmad, Nash'at; Proctor, Fred
2011-01-01
Convective storms can pose a serious risk to aviation operations since they are often accompanied by turbulence, heavy rain, hail, icing, lightning, strong winds, and poor visibility. They can cause major delays in air traffic due to the re-routing of flights, and by disrupting operations at the airports in the vicinity of the storm system. In this study, the Terminal Area Simulation System is used to simulate five different convective events ranging from a mesoscale convective complex to isolated storms. The occurrence of convection induced turbulence is analyzed from these simulations. The validation of model results with the radar data and other observations is reported and an aircraft-centric turbulence hazard metric calculated for each case is discussed. The turbulence analysis showed that large pockets of significant turbulence hazard can be found in regions of low radar reflectivity. Moderate and severe turbulence was often found in building cumulus turrets and overshooting tops.
NASA Astrophysics Data System (ADS)
Nejadmalayeri, Alireza; Vezolainen, Alexei; Vasilyev, Oleg V.
2013-11-01
In view of the ongoing longtime pursuit of numerical approaches that can capture important flow physics of high Reynolds number flows with fewest degrees of freedom, two important wavelet-based multi-resolution schemes are thoroughly examined, namely, the Coherent Vortex Simulation (CVS) and the Stochastic Coherent Adaptive Large Eddy Simulation (SCALES) with constant and spatially/temporarily variable thresholding. Reynolds number scaling of active spatial modes for CVS and SCALES of linearly forced homogeneous turbulence at high Reynolds numbers is investigated in dynamic study for the first time. This dynamic computational complexity study demonstrates that wavelet-based methods can capture flow-physics while using substantially fewer degrees of freedom than both direct numerical simulation and marginally resolved LES with the same level of fidelity or turbulence resolution, defined as ratio of subgrid scale and the total dissipations. The study provides four important observations: (1) the linear Reynolds number scaling of energy containing structures at a fixed level of kinetic energy, (2) small, close to unity, fractal dimension for constant-threshold CVS and SCALES simulations, (3) constant, close to two, fractal dimension for constant-dissipation SCALES that is insensitive to the level of fidelity, and (4) faster than quadratic decay of the compression ratio as a function of turbulence resolution. The very promising slope for Reynolds number scaling of CVS and SCALES demonstrates the potential of the wavelet-based methodologies for hierarchical multiscale space/time adaptive variable fidelity simulations of high Reynolds number turbulent flows.
NASA Astrophysics Data System (ADS)
Awasthi, Ankit; Anderson, William
2015-11-01
Large-scale motions in the logarithmic region of turbulent boundary layers amplitude modulate the viscous sublayer (Marusic et al., 2010: Science; Mathis et al., 2009: J. Fluid Mech.). This finding has promising implications for large-eddy simulation of wall-bounded turbulence at high Reynolds number (wherein the turbulence integral length exhibits linear proportionality with wall-normal elevation). Existing amplitude modulation studies have addressed smooth wall flows, though high Reynolds number rough wall flows are ubiquitous. Under such conditions, roughness-scale vortices ablate the viscous sublayer and result in the roughness sublayer. The roughness sublayer depth scales with aggregate element height, k, and is typically 2k ~ 3k. Above this, Townsend's Hypothesis dictates that the logarithmic layer is unaffected by the roughness sublayer. Here, we present large-eddy simulation results of turbulent channel flow over rough walls. We follow the decoupling procedure of Mathis et al., 2009: J. Fluid Mech., and present evidence that outer-layer dynamics amplitude modulate the roughness sublayer. Below the roughness element height, we report enormous sensitivity to element proximity. Above the elements, but within the roughness sublayer, topography dependence rapidly declines. This work was supported by the Air Force Office of Scientific Research, Turbulence and Transition Program (PM: Dr. R. Ponnoppan) under Grant # FA9550-14-1-0101. Computational resources were provided by the Texas Adv. Comp. Center at the Univ. of Texas.
An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder
NASA Technical Reports Server (NTRS)
Mankbadi, M. R.; Georgiadis, N. J.
2014-01-01
Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.
Large-eddy simulations of axisymmetric excitation effects on a row of impinging jets
NASA Astrophysics Data System (ADS)
Rizk, Magdi H.; Menon, Suresh
1988-07-01
Numerical simulations of a row of impinging jets are performed. Both the impinging jets and the fountains caused by the collision of the wall jets are modeled in the simulation. The problem considered contains the essential features of twin jets impinging on the ground, simulating the hovering configuration of a vertical takeoff and landing (VTOL) aircraft. The flow is assumed to be governed by the time-dependent, incompressible Navier-Stokes equations. The large-eddy simulation approach is followed. The present study focuses on the motion and dynamics of large-scale structures that have been experimentally observed in jet flows. The behavior of the jets and the fountain caused by the introduction of axisymmetric disturbances at the jet exits are investigated.
Large-eddy simulation of a turbulent flow over a heavy vehicle with drag reduction devices
NASA Astrophysics Data System (ADS)
Lee, Sangseung; Kim, Myeongkyun; You, Donghyun
2015-11-01
Aerodynamic drag contributes to a considerable amount of energy loss of heavy vehicles. To reduce the energy loss, drag reduction devices such as side skirts and boat tails, are often installed to the side and the rear of a heavy vehicle. In the present study, turbulent flow around a heavy vehicle with realistic geometric details is simulated using large-eddy simulation (LES), which is capable of providing unsteady flow physics responsible for aerodynamic in sufficient detail. Flow over a heavy vehicle with and without a boat tail and side skirts as drag reduction devices is simulated. The simulation results are validated against accompanying in-house experimental measurements. Effects of a boat tail and side skirts on drag reduction are discussed in detail. Supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) Grant NTIS 1615007940.
NASA Technical Reports Server (NTRS)
Givi, Peyman; Jaberi, Farhad A.
2001-01-01
The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.
Three regularization models as large-eddy simulations
NASA Astrophysics Data System (ADS)
Graham, Jonathan; Holm, Darryl; Mininni, Pablo; Pouquet, Annick
2006-11-01
We test three regularizations, the α-model, Leray-α, and Clark-α, as sub-grid models for LES by comparison with a 1024^3 direction numerical simulation (DNS), Rλ 800, with a Taylor-Green forcing. Both the α-model and Clark-α are able to reproduce the large-scale anisotropy of the flow as well as the time scale of developing turbulence. Leray-α fails in both these regards. We study intermittency corrections through pdfs and the anomalous scaling of the velocity increment structure functions. Leray-α is somewhat less intermittent than the DNS and produces an energy spectrum that is too shallow in the inertial range, while Clark-α produces a broad k-5/3 spectrum and stronger intermittency corrections. Finally, the agreement of the DNS and α-model spectra, in disparity with results for lower Reynolds number simulations, is worse than in the Clark-α model. We conjecture that this enhanced intermittency in the α model is related to the steeper than k-5/3 spectrum now reported for the very highest Reynolds number simulations and atmospheric observations.
Large Eddy Simulations of Kelvin-Helmholtz Instabilities in Stratified Ocean Flows
NASA Astrophysics Data System (ADS)
Brown, Dana; Goodman, Louis; Raessi, Mehdi
2012-11-01
Numerical simulations of turbulence in the ocean environment are used to supplement and enhance understanding of observational data. Here, using the NGA framework (Dejardins et al., JCP 2008), direct numerical simulations (DNS) and large eddy simulations (LES) of Kelvin-Helmholtz instabilities are employed to study turbulence in presence of density stratification. Kelvin-Helmholtz instabilities have been shown to be a common source of turbulence in the ocean. Past DNS studies of Kelvin Helmholtz instabilities have compared favorably with observational data, but were limited to moderate Reynolds numbers. Here, LES is used to solve the filtered incompressible NS equations at a higher Reynolds number, Re = 10,000. The effect of increased Reynolds number on the turbulence behavior is examined in terms of velocity spectra and energy budgets.
Formation mechanism of dust devil-like vortices in a large eddy simulation
NASA Astrophysics Data System (ADS)
Ito, J.; Niino, H.; Nakanishi, M.
2012-12-01
Dust devils are small-scale vertical vortices that often occur over deserts in fine weather conditions in which a convective mixed layer develops. Why such strong vortices are generated remains an issue in the dynamics of the atmospheric boundary layer, and several hypotheses for the origin of strong vertical vorticity in dust devils have been proposed. However, no quantitative study on the source of vertical vorticity of dust devils has been made. In this study, a large eddy simulation model with grid spacing of 5m is used to simulate dust devil-like vortices (DDVs) embedded in a convective mixed layer and a quantitative analysis on their source of its vertical vorticity is made. In order to investigate the origin of vertical vorticity in the simulated DDV, the circulation, which is a conserved quantity in the absence of turbulent transport and baroclinic production of horizontal vorticity, is examined, where the circulation is calculated as a surface integral of vorticity vector on a material surface. The deformation of the material surface as it flows into the DDV shows gives quantitative information about how stretching and tilting of vorticity contribute to the formation of the DDV. Material surface is initially placed horizontally in the core of the simulated DDV. It is divided into about 20000 triangular patches and vertices of the patches are tracked backward for 128 seconds. Our analysis shows that the material surface converges, while approximately conserving circulation, toward the DDV from a wide horizontal plane. A standard deviation of circulations over horizontal circles of several hundred meters in radius near the surface shows that presence of circulations is an inherent property of the convective mixed layer and its magnitude is reasonably scaled by the product of the depth of convective mixed layer and the convective velocity. As a result of horizontal convergence of the circulation, strength of formed DDVs can be scaled with the convective
Large-eddy simulation of supercritical fluid flow and combustion
NASA Astrophysics Data System (ADS)
Huo, Hongfa
The present study focuses on the modeling and simulation of injection, mixing, and combustion of real fluids at supercritical conditions. The objectives of the study are: (1) to establish a unified theoretical framework that can be used to study the turbulent combustion of real fluids; (2) to implement the theoretical framework and conduct numerical studies with the aim of improving the understanding of the flow and combustion dynamics at conditions representative of contemporary liquid-propellant rocket engine operation; (3) to identify the key design parameters and the flow variables which dictate the dynamics characteristics of swirl- and shear- coaxial injectors. The theoretical and numerical framework is validated by simulating the Sandia Flame D. The calculated axial and radial profiles of velocity, temperature, and mass fractions of major species are in reasonably good agreement with the experimental measurements. The conditionally averaged mass fraction profiles agree very well with the experimental results at different axial locations. The validated model is first employed to examine the flow dynamics of liquid oxygen in a pressure swirl injector at supercritical conditions. Emphasis is placed on analyzing the effects of external excitations on the dynamic response of the injector. The high-frequency fluctuations do not significantly affect the flow field as they are dissipated shortly after being introduced into the flow. However, the lower-frequency fluctuations are amplified by the flow. As a result, the film thickness and the spreading angle at the nozzle exit fluctuate strongly for low-frequency external excitations. The combustion of gaseous oxygen/gaseous hydrogen in a high-pressure combustion chamber for a shear coaxial injector is simulated to assess the accuracy and the credibility of the computer program when applied to a sub-scale model of a combustor. The predicted heat flux profile is compared with the experimental and numerical studies. The
The Jefferson Project: Large-eddy simulations of a watershed
NASA Astrophysics Data System (ADS)
Watson, C.; Cipriani, J.; Praino, A. P.; Treinish, L. A.; Tewari, M.; Kolar, H.
2015-12-01
The Jefferson Project is a new endeavor at Lake George, NY by IBM Research, Rensselaer Polytechnic Institute (RPI) and The Fund for Lake George. Lake George is an oligotrophic lake - one of low nutrients - and a 30-year study recently published by RPI's Darrin Fresh Water Institute highlighted the renowned water quality is declining from the injection of salt (from runoff), algae, and invasive species. In response, the Jefferson Project is developing a system to provide extensive data on relevant physical, chemical and biological parameters that drive ecosystem function. The system will be capable of real-time observations and interactive modeling of the atmosphere, watershed hydrology, lake circulation and food web dynamics. In this presentation, we describe the development of the operational forecast system used to simulate the atmosphere in the model stack, Deep ThunderTM (a configuration of the ARW-WRF model). The model performs 48-hr forecasts twice daily in a nested configuration, and in this study we present results from ongoing tests where the innermost domains are dx = 333-m and 111-m. We discuss the model's ability to simulate boundary layer processes, lake surface conditions (an input into the lake model), and precipitation (an input into the hydrology model) during different weather regimes, and the challenges of data assimilation and validation at this scale. We also explore the potential for additional nests over select regions of the watershed to better capture turbulent boundary layer motions.
Large-eddy simulation of turbulence in the free atmosphere and behind aircraft
NASA Astrophysics Data System (ADS)
Schumann, U.; Dörnbrack, A.; Dürbeck, T.; Gerz, T.
1997-02-01
The method of large-eddy simulation has been used for a wide variety of atmospheric flow problems. This paper gives an overview on recent applications of this method to turbulence in the free atmosphere under stably stratified conditions. In particular, flows in the wake of aircraft are studied in light of the potential impact of aircraft exhausts on the chemical and climatological state of the atmosphere. It is shown that different profiles of heat and moisture in the initial conditions of a jet representing engine exhaust gases may cause larger water saturation and hence earlier contrail formation than assumed up to now. The instability of trailing vortices in the wake of an aircraft is simulated up to the fully turbulent regime. The vertical diffusivity of aircraft exhaust is large in the vortex regime and much smaller than horizontal diffusivities in the later diffusion regime. The three-dimensional formation of a critical layer and breaking of gravity waves is simulated.
Large Eddy Simulation of complex sidearms subject to solar radiation and surface cooling.
Dittko, Karl A; Kirkpatrick, Michael P; Armfield, Steven W
2013-09-15
Large Eddy Simulation (LES) is used to model two lake sidearms subject to heating from solar radiation and cooling from a surface flux. The sidearms are part of Lake Audrey, NJ, USA and Lake Alexandrina, SA, Australia. The simulation domains are created using bathymetry data and the boundary is modelled with an Immersed Boundary Method. We investigate the cooling and heating phases with separate quasi-steady state simulations. Differential heating occurs in the cavity due to the changing depth. The resulting temperature gradients drive lateral flows. These flows are the dominant transport process in the absence of wind. Study in this area is important in water quality management as the lateral circulation can carry particles and various pollutants, transporting them to and mixing them with the main lake body. PMID:23863384
Large Eddy Simulation and PIV Visualization of a Vertical Hydrogen Jet
NASA Astrophysics Data System (ADS)
Pedro, G.; Peneau, F.; Wu, T. C.; Oshkai, P.; Djilali, N.
2006-11-01
Increasing concerns about green house gas emissions and deteriorating local air quality will necessitate substantial emission reductions, particularly from road vehicles. Canada has made important contributions in paving the way for the use of hydrogen in the transportation sector, which could lead to a substantial reduction of urban pollution and CO2 emissions. However production and storage issues, as well as the absence of specific standards for hydrogen are regarded as obstacles to the introduction of hydrogen in the energy market. A hydrogen jet exiting into quiescent air in both the supersonic and subsonic regimes is simulated using large eddy simulation with a Smagorinski subgrid model. The subsonic results are compared with experimental results obtained by Panchapakesan et al. and So et al. Using a high speed PIV system, a subsonic air-in-air jet is studied and the time averaged flow-field is compared to the one obtained in the simulation.
Modelling artificial sea salt emission in large eddy simulations.
Maalick, Z; Korhonen, H; Kokkola, H; Kühn, T; Romakkaniemi, S
2014-12-28
We study the dispersion of sea salt particles from artificially injected sea spray at a cloud-resolving scale. Understanding of how different aerosol processes affect particle dispersion is crucial when designing emission sources for marine cloud brightening. Compared with previous studies, we include for the first time an explicit treatment of aerosol water, which takes into account condensation, evaporation and their effect on ambient temperature. This enables us to capture the negative buoyancy caused by water evaporation from aerosols. Additionally, we use a higher model resolution to capture aerosol loss through coagulation near the source point. We find that, with a seawater flux of 15 kg s(-1), the cooling due to evaporation can be as much as 1.4 K, causing a delay in particle dispersion of 10-20 min. This delay enhances particle scavenging by a factor of 1.14 compared with simulations without aerosol water. We further show that both cooling and particle dispersion depend on the model resolution, with a maximum particle scavenging efficiency of 20% within 5 h after emission at maximum resolution of 50 m. Based on these results, we suggest further regional high-resolution studies which model several injection periods over several weeks. PMID:25404679
Modelling artificial sea salt emission in large eddy simulations
Maalick, Z.; Korhonen, H.; Kokkola, H.; Kühn, T.; Romakkaniemi, S.
2014-01-01
We study the dispersion of sea salt particles from artificially injected sea spray at a cloud-resolving scale. Understanding of how different aerosol processes affect particle dispersion is crucial when designing emission sources for marine cloud brightening. Compared with previous studies, we include for the first time an explicit treatment of aerosol water, which takes into account condensation, evaporation and their effect on ambient temperature. This enables us to capture the negative buoyancy caused by water evaporation from aerosols. Additionally, we use a higher model resolution to capture aerosol loss through coagulation near the source point. We find that, with a seawater flux of 15 kg s−1, the cooling due to evaporation can be as much as 1.4 K, causing a delay in particle dispersion of 10–20 min. This delay enhances particle scavenging by a factor of 1.14 compared with simulations without aerosol water. We further show that both cooling and particle dispersion depend on the model resolution, with a maximum particle scavenging efficiency of 20% within 5 h after emission at maximum resolution of 50 m. Based on these results, we suggest further regional high-resolution studies which model several injection periods over several weeks. PMID:25404679
NASA Astrophysics Data System (ADS)
Lee, Joon Sang
The compressible filtered Navier-Stokes equations were solved using a second order accurate finite volume method with low Mach number preconditioning. A dynamic subgrid-scale stress model accounted for the subgrid-scale turbulence. The study focused on the effects of buoyancy and rotation on the structure of turbulence and transport processes including heat transfer. Several different physical arrangements were studied as outlined below. The effects of buoyancy were first studied in a vertical channel using large eddy simulation (LES). The walls were maintained at constant temperatures, one heated and the other cooled. Results showed that aiding and opposing buoyancy forces emerge near the heated and cooled walls, respectively. In the aiding flow, the turbulent intensities and heat transfer were suppressed at large values of Grashof number. In the opposing flow, however, turbulence was enhanced with increased velocity fluctuations. Another buoyancy study considered turbulent flow in a vertically oriented annulus. Isoflux wall boundary conditions with low and high heating were imposed on the inner wall while the outer wall was adiabatic. The results showed that the strong heating and buoyancy force caused distortions of the flow structure resulting in reduction of turbulent intensities, shear stress, and turbulent heat flux, particularly near the heated wall. Flow in an annular pipe with and without an outer wall rotation about its axis was first investigated at moderate Reynolds numbers. When the outer pipe wall was rotated, a significant reduction of turbulent kinetic energy was realized near the rotating wall. Secondly, a large eddy simulation has been performed to investigate the effect of swirl on the heat and momentum transfer in an annular pipe flow with a rotating inner wall. The simulations indicated that the Nusselt number and the wall friction coefficient increased with increasing rotation speed of the wall. It was also observed that the axial velocity
Technical note: Large-eddy simulation of cloudy boundary layer with the Advanced Research WRF model
NASA Astrophysics Data System (ADS)
Yamaguchi, Takanobu; Feingold, Graham
2012-03-01
A thorough evaluation of the large-eddy simulation (LES) mode of the Advanced Research WRF model is performed with use of three cloudy boundary layer cases developed as LES intercomparison cases by the GEWEX Cloud System Study. Our evaluation reveals two problems that must be recognized and carefully addressed before proceeding with production runs. These are (i) sensitivity of results to the prescribed number of acoustic time steps per physical time step; and (ii) the assumption of saturation adjustment in the initial cloudy state. A temporary, but effective method of how to cope with these issues is suggested. With the proper treatment, the simulation results are comparable to the ensemble mean of the other LES models, and sometimes closer to the observational estimate than the ensemble mean. In order to ease the burden for configuration and post-processing, two new packages are developed and implemented. A detailed description of each package is presented. These packages are freely available to the public.
NASA Technical Reports Server (NTRS)
Jaberi, Farhad A.; Givi, Peyman
2003-01-01
The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.
Large-eddy simulation and multiscale modelling of a Richtmyer Meshkov instability with reshock
NASA Astrophysics Data System (ADS)
Hill, D. J.; Pantano, C.; Pullin, D. I.
2006-06-01
Large-eddy simulations of the Richtmyer Meshkov instability with reshock are pre- sented and the results are compared with experiments. Several configurations of shocks initially travelling from light (air) to heavy (sulfur hexafluoride, SF6) have been simulated to match previous experiments and good agreement is found in the growth rates of the turbulent mixing zone (TMZ). The stretched-vortex subgrid model used in this study allows for subgrid continuation modelling, where statistics of the unresolved scales of the flow are estimated. In particular, this multiscale modelling allows the anisotropy of the flow to be extended to the dissipation scale, eta, and estimates to be formed for the subgrid probability density function of the mixture fraction of air/SF6 based on the subgrid variance, including the effect of Schmidt number.
Large eddy simulation of flows after a bluff body: Coherent structures and mixing properties
NASA Astrophysics Data System (ADS)
Zhang, Pei; Han, Chao; Chen, Yiliang
2013-10-01
This paper performs large eddy simulations (LES) to investigate coherent structures in the flows after the Sydney bluff-body burner, a circular bluff body with an orifice at its center. The simulations are validated by comparison to existing experimental data. The Q function method is used to visualize the instantaneous vortex structures. Three kinds of structures are found, a cylindrical shell structure in the outer shear layer, a ring structure and some hairpin-like structures in the inner shear layer. An eduction scheme is employed to investigate the coherent structures in this flow. Some large streaks constituted by counter-rotating vortices are found in the outer shear layer and some well-organized strong structures are found in the inner shear layer. Finally, the influences of coherent structures on scalar mixing are studied and it is shown that scalar in the recirculation region is transported outward by coherent structures.
Further Aspects of Large Eddy Simulation Model Statistics and Inconsistencies with Field Data.
NASA Astrophysics Data System (ADS)
Agee, Ernest; Gluhovsky, Alexander
1999-08-01
This study has employed the concept of bandpass filtering of aircraft field data that are obtained within convective planetary boundary layers, for the purpose of comparing the characteristic turbulence statistics with the results of large eddy simulation (LES) models. Field data from Project LESS (10 January 1984) conducted over Lake Michigan in wintertime cold air outbreaks have been used to demonstrate the validity of the concept presented. These results show excellent agreement in variance statistics for LESS filtered data and LES model simulations; however, agreement in the skewness statistics is unsatisfactory. This is attributed to either the authors' inability to design a proper bandpass filter or the lack of an adequate dataset for computing third-moment statistics.
Large-eddy simulations of turbulent flow for grid-to-rod fretting in nuclear reactors
Bakosi, J.; Christon, M. A.; Lowrie, R. B.; Pritchett-Sheats, L. A.; Nourgaliev, R. R.
2013-07-12
The grid-to-rod fretting (GTRF) problem in pressurized water reactors is a flow-induced vibration problem that results in wear and failure of the fuel rods in nuclear assemblies. In order to understand the fluid dynamics of GTRF and to build an archival database of turbulence statistics for various configurations, implicit large-eddy simulations of time-dependent single-phase turbulent flow have been performed in 3 × 3 and 5 × 5 rod bundles with a single grid spacer. To assess the computational mesh and resolution requirements, a method for quantitative assessment of unstructured meshes with no-slip walls is described. The calculations have been carried out using Hydra-TH, a thermal-hydraulics code developed at Los Alamos for the Consortium for Advanced Simulation of Light water reactors, a United States Department of Energy Innovation Hub. Hydra-TH uses a second-order implicit incremental projection method to solve the singlephase incompressible Navier-Stokes equations. The simulations explicitly resolve the large scale motions of the turbulent flow field using first principles and rely on a monotonicity-preserving numerical technique to represent the unresolved scales. Each series of simulations for the 3 × 3 and 5 × 5 rod-bundle geometries is an analysis of the flow field statistics combined with a mesh-refinement study and validation with available experimental data. Our primary focus is the time history and statistics of the forces loading the fuel rods. These hydrodynamic forces are believed to be the key player resulting in rod vibration and GTRF wear, one of the leading causes for leaking nuclear fuel which costs power utilities millions of dollars in preventive measures. As a result, we demonstrate that implicit large-eddy simulation of rod-bundle flows is a viable way to calculate the excitation forces for the GTRF problem.
Large-eddy simulations of turbulent flow for grid-to-rod fretting in nuclear reactors
Bakosi, J.; Christon, M. A.; Lowrie, R. B.; Pritchett-Sheats, L. A.; Nourgaliev, R. R.
2013-07-12
The grid-to-rod fretting (GTRF) problem in pressurized water reactors is a flow-induced vibration problem that results in wear and failure of the fuel rods in nuclear assemblies. In order to understand the fluid dynamics of GTRF and to build an archival database of turbulence statistics for various configurations, implicit large-eddy simulations of time-dependent single-phase turbulent flow have been performed in 3 × 3 and 5 × 5 rod bundles with a single grid spacer. To assess the computational mesh and resolution requirements, a method for quantitative assessment of unstructured meshes with no-slip walls is described. The calculations have been carriedmore » out using Hydra-TH, a thermal-hydraulics code developed at Los Alamos for the Consortium for Advanced Simulation of Light water reactors, a United States Department of Energy Innovation Hub. Hydra-TH uses a second-order implicit incremental projection method to solve the singlephase incompressible Navier-Stokes equations. The simulations explicitly resolve the large scale motions of the turbulent flow field using first principles and rely on a monotonicity-preserving numerical technique to represent the unresolved scales. Each series of simulations for the 3 × 3 and 5 × 5 rod-bundle geometries is an analysis of the flow field statistics combined with a mesh-refinement study and validation with available experimental data. Our primary focus is the time history and statistics of the forces loading the fuel rods. These hydrodynamic forces are believed to be the key player resulting in rod vibration and GTRF wear, one of the leading causes for leaking nuclear fuel which costs power utilities millions of dollars in preventive measures. As a result, we demonstrate that implicit large-eddy simulation of rod-bundle flows is a viable way to calculate the excitation forces for the GTRF problem.« less
Large-eddy simulations of flows in a ramjet combustor
NASA Astrophysics Data System (ADS)
Jou, Wen-Huei; Menon, Suresh
The oscillatory cold flow in a ramjet combustor configuration is presently addressed by a numerical simulation method which gives attention to the interaction between the flowfield's vorticity and acoustic components, when the reduced frequency of the flow, based on the speed of sound, is of the order of unity. The numerical model has indicated that the combustor's interior must be isolated from the external region region by a choked nozzle. The numerical simulations thus obtained are able to exclude the effects of artificially imposed outflow-boundary conditions. The unsteady flow fields near the shear layer separation point in the nozzle region are investigated.
Scale-Similar Models for Large-Eddy Simulations
NASA Technical Reports Server (NTRS)
Sarghini, F.
1999-01-01
Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.
Wind Energy-Related Atmospheric Boundary Layer Large-Eddy Simulation Using OpenFOAM: Preprint
Churchfield, M.J.; Vijayakumar, G.; Brasseur, J.G.; Moriarty, P.J.
2010-08-01
This paper develops and evaluates the performance of a large-eddy simulation (LES) solver in computing the atmospheric boundary layer (ABL) over flat terrain under a variety of stability conditions, ranging from shear driven (neutral stratification) to moderately convective (unstable stratification).
Nesting Large-Eddy Simulations Within Mesoscale Simulations for Wind Energy Applications
NASA Astrophysics Data System (ADS)
Lundquist, J. K.; Mirocha, J. D.; Chow, F. K.; Kosovic, B.; Lundquist, K. A.
2008-12-01
With increasing demand for more accurate atmospheric simulations for wind turbine micrositing, for operational wind power forecasting, and for more reliable turbine design, simulations of atmospheric flow with resolution of tens of meters or higher are required. These time-dependent large-eddy simulations (LES) account for complex terrain and resolve individual atmospheric eddies on length scales smaller than turbine blades. These small-domain high-resolution simulations are possible with a range of commercial and open- source software, including the Weather Research and Forecasting (WRF) model. In addition to "local" sources of turbulence within an LES domain, changing weather conditions outside the domain can also affect flow, suggesting that a mesoscale model provide boundary conditions to the large-eddy simulations. Nesting a large-eddy simulation within a mesoscale model requires nuanced representations of turbulence. Our group has improved the Weather and Research Forecating model's (WRF) LES capability by implementing the Nonlinear Backscatter and Anisotropy (NBA) subfilter stress model following Kosoviæ (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005). We have also implemented an immersed boundary method (IBM) in WRF to accommodate complex terrain. These new models improve WRF's LES capabilities over complex terrain and in stable atmospheric conditions. We demonstrate approaches to nesting LES within a mesoscale simulation for farms of wind turbines in hilly regions. Results are sensitive to the nesting method, indicating that care must be taken to provide appropriate boundary conditions, and to allow adequate spin-up of turbulence in the LES domain. This work is performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
Large-eddy simulation of a turbulent mixing layer
NASA Technical Reports Server (NTRS)
Mansour, N. N.; Ferziger, J. H.; Reynolds, W. C.
1978-01-01
The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.
Minimum-dissipation models for large-eddy simulation
NASA Astrophysics Data System (ADS)
Bae, Hyunji Jane; Rozema, Wybe; Moin, Parviz; Verstappen, Roel
2015-11-01
Minimum-dissipation eddy-viscosity models are a class of subgrid scale models for LES that give the minimum eddy dissipation required to dissipate the energy of subgrid scales. The QR minimum-dissipation model [Verstappen, J. Sci. Comp., 2011] gives good results in simulations of decaying grid turbulence carried out on an isotropic grid. In particular, due to the minimum dissipation property of the model, the predicted energy spectra are in very good agreement with the DNS results up to the cut-off wave number unlike other methods. However, its results on anisotropic grids are often unsatisfactory because the model does not properly incorporate the grid anisotropy. We propose the anisotropic minimum-dissipation (AMD) model [Rozema et al., submitted for publication, 2015], a minimum-dissipation model that generalizes the QR model to anisotropic grids. The AMD model is more cost effective than the dynamic Smagorinsky model, appropriately switches off in laminar and transitional flow on anisotropic grids, and its subgrid scale model is consistent with the theoretic subgrid tensor. Experiments show that the AMD model is as accurate as the dynamic Smagorinsky model and Vreman model in simulations of isotropic turbulence, temporal mixing layer, and turbulent channel flow. H. J. Bae acknowledges support from SGF. W. Rozema and R. Verstappen acknowledge sponsoring by NWO for the use of supercomputing facilities and the financial support to attend the CTR SP 2014.
Nesting large-eddy simulations within mesoscale simulations for wind energy applications
Lundquist, J K; Mirocha, J D; Chow, F K; Kosovic, B; Lundquist, K A
2008-09-08
With increasing demand for more accurate atmospheric simulations for wind turbine micrositing, for operational wind power forecasting, and for more reliable turbine design, simulations of atmospheric flow with resolution of tens of meters or higher are required. These time-dependent large-eddy simulations (LES), which resolve individual atmospheric eddies on length scales smaller than turbine blades and account for complex terrain, are possible with a range of commercial and open-source software, including the Weather Research and Forecasting (WRF) model. In addition to 'local' sources of turbulence within an LES domain, changing weather conditions outside the domain can also affect flow, suggesting that a mesoscale model provide boundary conditions to the large-eddy simulations. Nesting a large-eddy simulation within a mesoscale model requires nuanced representations of turbulence. Our group has improved the Weather and Research Forecasting model's (WRF) LES capability by implementing the Nonlinear Backscatter and Anisotropy (NBA) subfilter stress model following Kosovic (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005). We have also implemented an immersed boundary method (IBM) in WRF to accommodate complex terrain. These new models improve WRF's LES capabilities over complex terrain and in stable atmospheric conditions. We demonstrate approaches to nesting LES within a mesoscale simulation for farms of wind turbines in hilly regions. Results are sensitive to the nesting method, indicating that care must be taken to provide appropriate boundary conditions, and to allow adequate spin-up of turbulence in the LES domain.
NASA Astrophysics Data System (ADS)
Li, X.
2014-12-01
Thermal stratification of the atmospheric surface layer has strong impact on the land-atmosphere exchange of turbulent, heat, and pollutant fluxes. Few studies have been carried out for the interaction of the weakly to moderately stable stratified atmosphere and the urban canopy. This study performs a large-eddy simulation of a modeled street canyon within a weakly to moderately stable atmosphere boundary layer. To better resolve the smaller eddy size resulted from the stable stratification, a higher spatial and temporal resolution is used. The detailed flow structure and turbulence inside the street canyon are analyzed. The relationship of pollutant dispersion and Richardson number of the atmosphere is investigated. Differences between these characteristics and those under neutral and unstable atmosphere boundary layer are emphasized.
Large-Eddy Simulation of Chemically Reactive Pollutant Transport from a Point Source in Urban Area
NASA Astrophysics Data System (ADS)
Du, Tangzheng; Liu, Chun-Ho
2013-04-01
Most air pollutants are chemically reactive so using inert scalar as the tracer in pollutant dispersion modelling would often overlook their impact on urban inhabitants. In this study, large-eddy simulation (LES) is used to examine the plume dispersion of chemically reactive pollutants in a hypothetical atmospheric boundary layer (ABL) in neutral stratification. The irreversible chemistry mechanism of ozone (O3) titration is integrated into the LES model. Nitric oxide (NO) is emitted from an elevated point source in a rectangular spatial domain doped with O3. The LES results are compared well with the wind tunnel results available in literature. Afterwards, the LES model is applied to idealized two-dimensional (2D) street canyons of unity aspect ratio to study the behaviours of chemically reactive plume over idealized urban roughness. The relation among various time scales of reaction/turbulence and dimensionless number are analysed.
Wall-modeling for large-eddy simulation of high Reynolds number supersonic flows
NASA Astrophysics Data System (ADS)
Kawai, Soshi; Larsson, Johan; Lele, Sanjiva
2010-11-01
We present an idea of approximate wall-boundary-condition approach with dynamic procedure for large-eddy simulation of Mach 3 supersonic turbulent boundary layer at various Reynolds numbers (Reδ=2 x10^4, 10^5 and 10^6) on a flat plate. This wall-model is the extension of previous work by Wang and Moin [Phys. Fluid, 14, 2043 (2002)] for incompressible flows to compressible flows. We note that the present study is both the first extension of the dynamic concept to compressible flows and also the first test at high Reynolds number flows. The present study also revisits the issue of numerical errors near wall-region on outer-layer coarse LES mesh. The numerical results are compared with wall-resolved LES data (at low Reynolds number case) and available experimental data (at high Reynolds number case).
Large Eddy Simulation in the Computation of Jet Noise
NASA Technical Reports Server (NTRS)
Mankbadi, R. R.; Goldstein, M. E.; Povinelli, L. A.; Hayder, M. E.; Turkel, E.
1999-01-01
Noise can be predicted by solving Full (time-dependent) Compressible Navier-Stokes Equation (FCNSE) with computational domain. The fluctuating near field of the jet produces propagating pressure waves that produce far-field sound. The fluctuating flow field as a function of time is needed in order to calculate sound from first principles. Noise can be predicted by solving the full, time-dependent, compressible Navier-Stokes equations with the computational domain extended to far field - but this is not feasible as indicated above. At high Reynolds number of technological interest turbulence has large range of scales. Direct numerical simulations (DNS) can not capture the small scales of turbulence. The large scales are more efficient than the small scales in radiating sound. The emphasize is thus on calculating sound radiated by large scales.
Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation
NASA Technical Reports Server (NTRS)
Kim, J.; Moin, P.
1979-01-01
The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').
High Speed Jet Noise Prediction Using Large Eddy Simulation
NASA Technical Reports Server (NTRS)
Lele, Sanjiva K.
2002-01-01
Current methods for predicting the noise of high speed jets are largely empirical. These empirical methods are based on the jet noise data gathered by varying primarily the jet flow speed, and jet temperature for a fixed nozzle geometry. Efforts have been made to correlate the noise data of co-annular (multi-stream) jets and for the changes associated with the forward flight within these empirical correlations. But ultimately these emipirical methods fail to provide suitable guidance in the selection of new, low-noise nozzle designs. This motivates the development of a new class of prediction methods which are based on computational simulations, in an attempt to remove the empiricism of the present day noise predictions.
Large-eddy simulations of a propelled submarine model
NASA Astrophysics Data System (ADS)
Posa, Antonio; Balaras, Elias
2015-11-01
The influence of the propeller on the wake as well as the evolution of the turbulent boundary layers over an appended notional submarine geometry (DARPA SUBOFF) is reported. The present approach utilizes a wall-resolved LES, coupled with an immersed boundary formulation, to simulate the flow model scale Reynolds numbers (Re = 1 . 2 e + 06 , based on the free-stream velocity and the length of the body). Cylindrical coordinates are adopted, and the computational grid is composed of 3.5 billion nodes. Our approach has been validated on the appended submarine body in towed conditions (without propeller), by comparisons to wind tunnel experiments in the literature. The comparison with the towed configuration shows profound modifications in the boundary layer over the stern surface, due to flow acceleration, with higher values of turbulent kinetic energy in the inner layer and lower values in the outer layer. This behavior was found tied to a different topology of the coherent structures between propelled and towed cases. The wake is also highly affected, and the momentum deficit displays a non-monotonic evolution downstream. An axial peak of turbulent kinetic energy replaces the bimodal distribution of the stresses in the wake, observed in the towed configuration. Supported by ONR Grant N000141110455, monitored by Dr. Ki-Han Kim.
NASA Technical Reports Server (NTRS)
Baurle, R. A.
2015-01-01
Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit
Large-Eddy Simulation of plume dispersion within various actual urban areas
NASA Astrophysics Data System (ADS)
Nakayama, H.; Jurcakova, K.; Nagai, H.
2013-02-01
Plume dispersion of hazardous materials within urban area resulting from accidental or intentional releases is of great concern to public health. Many researchers have developed local-scale atmospheric dispersion models using building-resolving computational fluid dynamics. However, an important issue is encountered when determining a reasonable domain size of the computational model in order to capture concentration distribution patterns influenced by urban surface geometries. In this study, we carried out Large-Eddy Simulations (LES) of plume dispersion within various urban areas with a wide range of obstacle density and building height variability. The difference of centerline mean and r.m.s. concentration distributions among various complex urban surface geometries becomes small for downwind distances from the point source greater than 1.0 km. From these results, it can be concluded that a length of a computational model should be at least 1.0 km from a point source.
NASA Technical Reports Server (NTRS)
Spyropoulos, Evangelos T.; Holmes, Bayard S.
1997-01-01
The dynamic subgrid-scale model is employed in large-eddy simulations of flow over a cylinder at a Reynolds number, based on the diameter of the cylinder, of 90,000. The Centric SPECTRUM(trademark) finite element solver is used for the analysis. The far field sound pressure is calculated from Lighthill-Curle's equation using the computed fluctuating pressure at the surface of the cylinder. The sound pressure level at a location 35 diameters away from the cylinder and at an angle of 90 deg with respect to the wake's downstream axis was found to have a peak value of approximately 110 db. Slightly smaller peak values were predicted at the 60 deg and 120 deg locations. A grid refinement study suggests that the dynamic model demands mesh refinement beyond that used here.
Large Eddy Simulations of Turbulent Reacting Flows in an Opposed-Piston Two Stroke Engine
NASA Astrophysics Data System (ADS)
Srivastava, Shalabh; Schock, Harold; Jaberi, Farhad
2013-11-01
The two-phase filtered mass density function (FMDF) subgrid-scale model has been used for large eddy simulation (LES) of turbulent spray combustion in a generic single cylinder, opposed-piston, two-stroke engine configuration. The LES/FMDF is implemented via an efficient, hybrid numerical method in which the filtered compressible Navier-Stokes equations are solved with a high-order, multi-block, compact differencing scheme, and the spray and FMDF are implemented with stochastic Lagrangian methods. The reliability and consistency of the numerical methods are established for the engine configuration by comparing the Eulerian and Lagrangian components of the LES/FMDF. The effects of various operating conditions like boost pressure, heat transfer model, fuel spray temperature, nozzle diameter, injection pressure, and injector configuration on the flow field, heat loss and the evolution of spray and combustion are studied.
Large-eddy simulations of adverse pressure gradient turbulent boundary layers
NASA Astrophysics Data System (ADS)
Bobke, Alexandra; Vinuesa, Ricardo; Örlü, Ramis; Schlatter, Philipp
2016-04-01
Adverse pressure-gradient (APG) turbulent boundary layers (TBL) are studied by performing well-resolved large-eddy simulations. The pressure gradient is imposed by defining the free-stream velocity distribution with the description of a power law. Different inflow conditions, box sizes and upper boundary conditions are tested in order to determine the final set-up. The statistics of turbulent boundary layers with two different power-law coefficients and thus magnitudes of adverse pressure gradients are then compared to zero pressure-gradient (ZPG) data. The effect of the APG on TBLs is manifested in the mean flow through a much more prominent wake region and in the Reynolds stresses through the existence of an outer peak. The pre-multiplied energy budgets show that more energy is transported from the near-wall region to farther away from the wall.
Using Large Eddy Simulation for understanding vented gas explosions in the presence of obstacles.
Di Sarli, Valeria; Di Benedetto, Almerinda; Russo, Gennaro
2009-09-30
In this work, a validated Large Eddy Simulation model of unsteady premixed flame propagation is used to study the phenomenology underlying vented gas explosions in the presence of obstacles. Computations are run of deflagrating flames in a small-scale combustion chamber closed at the bottom end and open at the opposite face. A single obstacle is centred inside the chamber. Methane-air mixtures of various compositions (ranging from lean to stoichiometric and rich), and obstacles with different area blockage ratios (30, 50 and 70%) and shapes (circular, rectangular and square cross-section in the flow direction) are investigated. All cases are initialized from stagnation. The competition between combustion rate and venting rate allows explaining both number and intensity of the overpressure peaks observed. PMID:19409700
Large-eddy simulation of plume dispersion under various thermally stratified boundary layers
NASA Astrophysics Data System (ADS)
Nakayama, H.; Takemi, T.; Nagai, H.
2014-07-01
Contaminant gas dispersion in atmospheric boundary layer is of great concern to public health. For the accurate prediction of the dispersion problem, the present study numerically investigates the behavior of plume dispersion by taking into account the atmospheric stability which is classified into three types; neutral, stable, and convective boundary layers. We first proposed an efficient method to generate spatially-developing, thermally-stratified boundary layers and examined the usefulness of our approach by comparing to wind tunnel experimental data for various thermal boundary layers. The spreads of plume in the spanwise direction are quantitatively underestimated especially at large downwind distances from the point source, owing to the underestimation of turbulence intensities for the spanwise component; however, the dependence of the spanwise spreads to atmospheric stability is well represented in a qualitative sense. It was shown that the large-eddy simulation (LES) model provides physically reasonable results.
Large Eddy simulation of parallel blade-vortex interaction
NASA Astrophysics Data System (ADS)
Felten, Frederic; Lund, Thomas
2002-11-01
Helicopter Blade-Vortex Interaction (BVI) generally occurs under certain conditions of powered descent or during extreme maneuvering. The vibration and acoustic problems associated with the interaction of rotor tip vortices and the following blades is a major aerodynamic concern for the helicopter community. Numerous experimental and computational studies have been done over the last two decades in order to gain a better understanding of the physical mechanisms involved in BVI. The most severe interaction, in terms of generated noise, happens when the vortex filament is parallel to the blade, thus affecting a great portion of it. The majority of the previous numerical studies of parallel BVI fall within a potential flow framework. Some Navier-Stokes approaches using dissipative numerical methods and RANS-type turbulence models have also been attempted, but with limited success. The current investigation makes use of an incompressible, non-dissipative, kinetic energy conserving collocated mesh scheme in conjunction with a dynamic subgrid-scale model. The concentrated tip vortex is not attenuated as it is convected downstream and over a NACA-0012 airfoil. The lift, drag, moment and pressure coefficients induced by the passage of the vortex are monitored in time and compared with experimental data.
Large eddy simulations of coal jet flame ignition using the direct quadrature method of moments
NASA Astrophysics Data System (ADS)
Pedel, Julien
The Direct Quadrature Method of Moments (DQMOM) was implemented in the Large Eddy Simulation (LES) tool ARCHES to model coal particles. LES coupled with DQMOM was first applied to nonreacting particle-laden turbulent jets. Simulation results were compared to experimental data and accurately modeled a wide range of particle behaviors, such as particle jet waviness, spreading, break up, particle clustering and segregation, in different configurations. Simulations also accurately predicted the mean axial velocity along the centerline for both the gas phase and the solid phase, thus demonstrating the validity of the approach to model particles in turbulent flows. LES was then applied to the prediction of pulverized coal flame ignition. The stability of an oxy-coal flame as a function of changing primary gas composition (CO2 and O2) was first investigated. Flame stability was measured using optical measurements of the flame standoff distance in a 40 kW pilot facility. Large Eddy Simulations (LES) of the facility provided valuable insight into the experimentally observed data and the importance of factors such as heterogeneous reactions, radiation or wall temperature. The effects of three parameters on the flame stand-off distance were studied and simulation predictions were compared to experimental data using the data collaboration method. An additional validation study of the ARCHES LES tool was then performed on an air-fired pulverized coal jet flame ignited by a preheated gas flow. The simulation results were compared qualitatively and quantitatively to experimental observations for different inlet stoichiometric ratios. LES simulations were able to capture the various combustion regimes observed during flame ignition and to accurately model the flame stand-off distance sensitivity to the stoichiometric ratio. Gas temperature and coal burnout predictions were also examined and showed good agreement with experimental data. Overall, this research shows that high
Large-eddy simulation of free-surface decaying turbulence with dynamic subgrid-scale models
NASA Astrophysics Data System (ADS)
Salvetti, M. V.; Zang, Y.; Street, R. L.; Banerjee, S.
1997-08-01
This paper describes large-eddy simulations of decaying turbulence in an open channel, using different dynamic subgrade-scale models, viz. the dynamic model of Germano et al. [Phys. Fluids A 3, 1790 (1991)] (DSM), the dynamic mixed model in Zang et al. [Phys. Fluids A 5, 3186 (1993)] (DMM), and the dynamic two-parameter model of Salvetti and Banerjee [Phys. Fluids 7, 2831 (1995)] (DTM). These models are incorporated in a finite-volume solver of the Navier-Stokes equations. A direct numerical simulation of this flow conducted by Pan and Banerjee [Phys. Fluids 7, 1649 (1995)] showed that near the free surface turbulence has a quasi-two-dimensional behavior. Moreover, the quasi-two-dimensional region increases in thickness with the decay time, although the structure remains three-dimensional in the central regions of the flow. The results of the large-eddy simulations show that both the DMM and the DTM are able to reproduce the features of the decay process observed in the direct simulation and to handle the anisotropic nature of the flow. Nevertheless, the addition of the second model coefficient in the DTM improves the agreement with the direct simulation. When the DSM is used, significant discrepancies are observed between the large-eddy and the direct simulations during the decay process at the free surface.
Large-Eddy Simulations of Plasma Flow Control on a GOE735 Wind Turbine Airfoil
NASA Astrophysics Data System (ADS)
Czulak, Alexander; Franck, Jennifer
2015-11-01
Active flow control using plasma actuation was studied for the GOE735 airfoil and compared to non-actuated baseline cases using numerical simulations. This investigation considers two-dimensional simulations at a Reynolds number of 1,000 using direct numerical simulation (DNS) as well as three-dimensional simulations at a Reynolds number of 50,000 and 100,000 using large-eddy simulation (LES). Plasma actuation is applied in terms of a source term within the boundary layer close to the airfoil surface. Angles of attack of 0°, 5° and 15° were considered, and control is shown to be effective at increasing the lift coefficient, decreasing the drag coefficient and reducing the root mean squared deviation of both lift and drag. An analysis of the flow physics reveals that the actuated cases delay the point of separation, reduce the wake width and diminish the size and strength of the shed vortices. For this particular airfoil, there are significant differences in Reynolds number in terms of the baseline flow, control effectiveness and performance factors such as lift and drag.
Large-eddy numerical simulation of an array of three-dimensional impinging jets
NASA Astrophysics Data System (ADS)
Rizk, M. H.; Menon, S.
1987-04-01
Numerical simulations of a row of impinging jets are performed. Both the impinging jets and the fountains caused by the collision of the wall jets are modeled in the simulations. The problem considered contains the essential features of twin jets impinging on the ground, simulating the hovering configuration of a VTOL aircraft. The flow is assumed to be governed by the time-dependent, incompressible Navier-Stokes equations. The large-eddy simulation approach is followed in which all scales resolvable by the grid resolution are computed explicitly, while the small-scale turbulence structures, which are nearly universal in character, are modeled by an eddy viscosity formulation that simulates the energy cascade into the small scales. The Navier-Stokes equations are solved using a staggered computational mesh. Central finite differencing is used to discretize all terms except the convective terms, which are discretized using the QUICK scheme. The Adams-Bashforth scheme is used to advance the solution in time. The pressure Poisson equation is used in place of the continuity equation. Efficient direct solutions are obtained for the pressure field, which allows the continuity equation to be satisfied at each time step. This study focuses on the motion and dynamics of large-scale structures that have been experimentally observed in jet flows. The behavior of the jets and the fountain due to introducing axisymmetric, azimuthal and random disturbances at the jet exits is investigated.
Aero-Thermal Prediction in High Pressure Turbine Cascade using Large Eddy Simulation
NASA Astrophysics Data System (ADS)
Bhaskaran, Rathakrishnan; Lele, Sanjiva
2008-11-01
The aero-thermal performance of an uncooled, smooth high pressure (HP) turbine cascade in the presence of free-stream turbulence is studied using a high-order overset mesh Large Eddy Simulation (LES) procedure. A HP vane cascade designed at the von Karman Institute (VKI) for fluid dynamics, Belgium, is used as the model geometry. Simulations matching experimental conditions, except for the Reynolds number which is about half of the experimental value, have been carried out. Significant enhancement in the blade heat-transfer is seen in the presence of inflow turbulence. Eddies from the free-stream turbulence get stretched around the blade, creating long streaky structures in the blade boundary layer. These structures quickly break down on the suction side, while they persist on the pressure side. The blade heat transfer signature from the simulations does not show transition of the boundary layer at the Reynolds number of the simulation. This is consistent with the trend seen in the experiments where transition is delayed by lowering the Reynolds number. New simulations matching the experimental Reynolds number are currently under way.
Applications of large-eddy simulation: Synthesis of neutral boundary layer models
Ohmstede, W.D.
1987-12-01
The object of this report is to describe progress made towards the application of large-eddy simulation (LES), in particular, to the study of the neutral boundary layer (NBL). The broad purpose of the study is to provide support to the LES project currently underway at LLNL. The specific purpose of this study is to lay the groundwork for the simulation of the SBL through the establishment and implementation of model criteria for the simulation of the NBL. The idealistic NBL is never observed in the atmosphere and therefore has little practical significance. However, it is of considerable theoretical interest for several reasons. The report discusses the concept of Rossby-number similarity theory as it applies to the NBL. A particular implementation of the concept is described. Then, the results from prior simulations of the NBL are summarized. Model design criteria for two versions of the Brost LES (BLES) model are discussed. The general guidelines for the development of Version 1 of the Brost model (BV1) were to implement the model with a minimum of modifications which would alter the design criteria as established by Brost. Two major modifications of BLES incorporated into BV1 pertain to the initialization/parameterization of the model and the generalization of the boundary conditions at the air/earth interface. 18 refs., 4 figs.
NASA Astrophysics Data System (ADS)
Le, Trung; Khosronejad, Ali; Bartelt, Nicole; Woldeamlak, Solomon; Peterson, Bonnie; Dewall, Petronella; Sotiropoulos, Fotis; Saint Anthony Falls Laboratory, University of Minnesota Team; Minnesota Department of Transportation Team
2015-11-01
We study the dynamics of a river confluence on Mississippi River branch in the city of Minneapolis, Minnesota, United States. Field measurements by Acoustic Doppler Current Profiler using on-board GPS tracking were carried out for five campaigns in the summer of 2014 and 2015 to collect both river bed elevation data and flow fields. Large Eddy Simulation is carried out to simulate the flow field with the total of 100 million grid points for the domain length of 3.2 km. The simulation results agree well with field measurements at measured cross-sections. The results show the existence of wake mode on the mixing interface of two branches near the upstream junction corner. The mutual interaction between the shear layers emanating from the river banks leading to the formation of large scale energetic structures that leads to ``switching'' side of the flow coherent structures. Our result here is a feasibility study for the use of eddy-resolving simulations in predicting complex flow dynamics in medium-size natural rivers. This work is funded by Minnesota Dept. Transportation and Minnesota Institute of Supercomputing.
Mean-state acceleration of cloud-resolving models and large eddy simulations
Jones, C. R.; Bretherton, C. S.; Pritchard, M. S.
2015-10-29
In this study, large eddy simulations and cloud-resolving models (CRMs) are routinely used to simulate boundary layer and deep convective cloud processes, aid in the development of moist physical parameterization for global models, study cloud-climate feedbacks and cloud-aerosol interaction, and as the heart of superparameterized climate models. These models are computationally demanding, placing practical constraints on their use in these applications, especially for long, climate-relevant simulations. In many situations, the horizontal-mean atmospheric structure evolves slowly compared to the turnover time of the most energetic turbulent eddies. We develop a simple scheme to reduce this time scale separation to accelerate themore » evolution of the mean state. Using this approach we are able to accelerate the model evolution by a factor of 2–16 or more in idealized stratocumulus, shallow and deep cumulus convection without substantial loss of accuracy in simulating mean cloud statistics and their sensitivity to climate change perturbations. As a culminating test, we apply this technique to accelerate the embedded CRMs in the Superparameterized Community Atmosphere Model by a factor of 2, thereby showing that the method is robust and stable to realistic perturbations across spatial and temporal scales typical in a GCM.« less
Mean-state acceleration of cloud-resolving models and large eddy simulations
Jones, C. R.; Bretherton, C. S.; Pritchard, M. S.
2015-10-29
In this study, large eddy simulations and cloud-resolving models (CRMs) are routinely used to simulate boundary layer and deep convective cloud processes, aid in the development of moist physical parameterization for global models, study cloud-climate feedbacks and cloud-aerosol interaction, and as the heart of superparameterized climate models. These models are computationally demanding, placing practical constraints on their use in these applications, especially for long, climate-relevant simulations. In many situations, the horizontal-mean atmospheric structure evolves slowly compared to the turnover time of the most energetic turbulent eddies. We develop a simple scheme to reduce this time scale separation to accelerate the evolution of the mean state. Using this approach we are able to accelerate the model evolution by a factor of 2–16 or more in idealized stratocumulus, shallow and deep cumulus convection without substantial loss of accuracy in simulating mean cloud statistics and their sensitivity to climate change perturbations. As a culminating test, we apply this technique to accelerate the embedded CRMs in the Superparameterized Community Atmosphere Model by a factor of 2, thereby showing that the method is robust and stable to realistic perturbations across spatial and temporal scales typical in a GCM.
Analysis on Turbulent Flows using Large-eddy Simulation on the Seaside Complex Terrain
NASA Astrophysics Data System (ADS)
Kamio, T.; Iida, M.; Arakawa, C.
2014-12-01
The purpose of this study is the Large-eddy Simulation (LES) of the turbulent wind on the complex terrain, and the first results of the simulation are described. The authors tried to apply the LES code, which was developed as an atmospheric simulator in Japan Agency for the Marine-Earth Science and Technology (JAMSTEC), to the wind prediction for the wind energy. On the wind simulation, the highest problem would be the boundary conditions, and the case in this paper was simplified one. The case study in this paper is the west wind on a complex terrain site, which is the wind from sea for the site. The steady flow was employed for the inlet condition, because the wind on the sea is the low turbulent wind, and almost all the turbulence would be generated by the roughness of the ground surface. The wall function was employed as the surface condition on the ground surface. The computational domain size was about 8 × 3 × 2.5 km3, and the minimum cell size was about 10 × 10 × 3 m3. The computational results, the vertical profile of the averaged wind speed and the turbulence intensity, agreed with the measurement by the meteorological masts. Moreover, the authors tried the analysis of the turbulence characteristics. The power spectrum density model, and the cross spectrum analyses gave the knowledge of the turbulent characteristics on the complex terrain and the hints for the domain and grid of the numerical analysis.
A modified Darcy's law . Large eddy simulation of turbulent flows through a fractal model city
NASA Astrophysics Data System (ADS)
Gisinger, Sonja; Dörnbrack, Andreas; Schröttle, Josef
2015-08-01
An approach to describe the turbulent flow through a complex geometry (e.g., urban area) by means of an analogy to flows through porous media is presented. Therefore, a modification of the original Darcy's law is proposed, and its application is tested in a prototype problem with an idealized complex geometry using large eddy simulations. The numerical results indicate the validity of the modified Darcy's law for the chosen setup.
Direct-Numerical and Large-Eddy Simulations of a Non-Equilibrium Turbulent Kolmogorov Flow
NASA Technical Reports Server (NTRS)
Woodruff, S. L.; Shebalin, J. V.; Hussaini, M. Y.
1999-01-01
A non-equilibrium form of turbulent Kolmogorov flow is set up by making an instantaneous change in the amplitude of the spatially-periodic forcing. It is found that the response of the flow to this instantaneous change becomes more dramatic as the wavenumber of the forcing is increased, and, at the same time, that the faithfulness with which the large-eddy-simulation results agree with the direct-numerical results decreases.
Large-eddy simulation of bubble-driven plume in stably stratified flow.
NASA Astrophysics Data System (ADS)
Yang, Di; Chen, Bicheng; Socolofsky, Scott; Chamecki, Marcelo; Meneveau, Charles
2015-11-01
The interaction between a bubble-driven plume and stratified water column plays a vital role in many environmental and engineering applications. As the bubbles are released from a localized source, they induce a positive buoyancy flux that generates an upward plume. As the plume rises, it entrains ambient water, and when the plume rises to a higher elevation where the stratification-induced negative buoyancy is sufficient, a considerable fraction of the entrained fluid detrains, or peels, to form a downward outer plume and a lateral intrusion layer. In the case of multiphase plumes, the intrusion layer may also trap weakly buoyant particles (e.g., oil droplets in the case of a subsea accidental blowout). In this study, the complex plume dynamics is studied using large-eddy simulation (LES), with the flow field simulated by hybrid pseudospectral/finite-difference scheme, and the bubble and dye concentration fields simulated by finite-volume scheme. The spatial and temporal characteristics of the buoyant plume are studied, with a focus on the effects of different bubble buoyancy levels. The LES data provide useful mean plume statistics for evaluating the accuracy of 1-D engineering models for entrainment and peeling fluxes. Based on the insights learned from the LES, a new continuous peeling model is developed and tested. Study supported by the Gulf of Mexico Research Initiative (GoMRI).
Large-eddy simulation of flow over the Great Plains under stable atmospheric conditions
NASA Astrophysics Data System (ADS)
Zhou, B.; Chow, F. K.
2010-12-01
The Great Plains in the central part of North America hosts enormous wind resources. One of the key meteorological features over the Great Plains is the frequent occurrence of nocturnal low-level jets under stably-stratified conditions. The flow speed up due to the formation of the low-level jets represents great wind power potential. In this study, large-eddy simulations (LES) will be performed over the site where the Cooperative Atmospheric-surface Exchange Study (CASES-99) field experiment took place. Atmospheric boundary layer (ABL) simulations driven by both strongly and weakly forced synoptic flows under stable atmospheric conditions will be investigated. While continuous turbulence is expected under strongly forced conditions, the weakly forced scenario is likely intermittent in nature, with occasional elevated turbulent bursts. The focus of this study includes vertical wind shear profiles, as well as turbulent statistics under stable conditions over the relatively flat, yet complex terrain. We will use an explicit filtering and reconstruction turbulence modeling LES approach. This approach has been proven advantageous in our previous work in terms of turbulence representation and agreement with similarity theory in neutral and stable atmospheric boundary layer flow over flat terrain. The dynamic reconstruction turbulence closure is capable of handling strong atmospheric stability, and predicting intermittent turbulence burst events in previous idealized simulations. This LES study ill provide detailed flow features under stable conditions over the Great Plains that can be valuable to the wind energy industry.
Large-Eddy Simulations of Fuel-Air Mixing and Combustion in an Internal Combustion Engine
NASA Astrophysics Data System (ADS)
Sone, Kazuo; Menon, Suresh
2000-11-01
Past studies of internal combustion (IC) engine steady-state flow field have employed the well-known KIVA code for steady-state predictions. However, it is also well known that this code is incapable of accurately capturing the impact of unsteady fuel-air mixing on the combustion process. Here, the latest KIVA-3V code has been modified to carry out large-eddy simulations (LES). In particular, the RANS k-e model has been replaced by a subgrid kinetic energy model and a fourth-order ENO scheme has been implemented to increase the accuracy of the discretization of the advection term. Finally, a subgrid model to simulate the small-scale turbulent mixing, combustion and heat release is implemented for reacting flows. Simulations using the new version of KIVA3V denoted here as KIVALES of temporal mixing layers and flows past rearward facing step demonstrate the improved accuracy of the LES model. Accuracy of the prediction is demonstrated by comparing with DNS, LES and experimental results obtained in the past. Finally, the new code is employed to simulate fuel-air mixing and combustion in a typical IC engine. Comparison with predictions using the conventional KIVA is used to demonstrate the ability of the new code.
Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder
Morgan, B. E.; Greenough, J. A.
2015-04-08
Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the time at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.
Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder
Morgan, B. E.; Greenough, J. A.
2015-04-08
Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the timemore » at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.« less
Large-eddy simulation of sea breeze at an idealized peninsular site
NASA Astrophysics Data System (ADS)
Rizza, Umberto; Miglietta, Mario Marcello; Anabor, Vagner; Degrazia, Gervasio A.; Maldaner, Silvana
2015-08-01
A high-resolution large-eddy simulation (LES) has been performed to simulate a sea-breeze circulation over an idealized peninsular domain. The simulation is forced with the surface latent/sensible heat fluxes and the large-scale horizontal pressure gradient that are obtained from a mesoscale simulation. This methodology allows the investigation of the physical phenomena that are peculiar for a sea-breeze circulation and that generally require spatial resolution approximately equal to 100 m or less. Here, small-scale dynamical effects associated to these phenomena, i.e. the interaction between the sea-breeze front with the convective turbulence generated over-land, the formation of the zero-velocity layer, and the development of the Kelvin-Helmholtz billows, are investigated. Results from the present numerical study have revealed the formation of a zero-velocity layer that is initially near the ground then it rises to define a well-marked sea-breeze depth. Scaling analysis applied to the LES output fields reveals that during the phase of inland penetration the scaling laws for sea-breeze strength and depth have both a proportionality coefficient equal to 0.15.
Large Eddy/Reynolds-Averaged Navier-Stokes Simulations of CUBRC Base Heating Experiments
NASA Technical Reports Server (NTRS)
Salazar, Giovanni; Edwards, Jack R.; Amar, Adam J.
2012-01-01
ven with great advances in computational techniques and computing power during recent decades, the modeling of unsteady separated flows, such as those encountered in the wake of a re-entry vehicle, continues to be one of the most challenging problems in CFD. Of most interest to the aerothermodynamics community is accurately predicting transient heating loads on the base of a blunt body, which would result in reduced uncertainties and safety margins when designing a re-entry vehicle. However, the prediction of heat transfer can vary widely depending on the turbulence model employed. Therefore, selecting a turbulence model which realistically captures as much of the flow physics as possible will result in improved results. Reynolds Averaged Navier Stokes (RANS) models have become increasingly popular due to their good performance with attached flows, and the relatively quick turnaround time to obtain results. However, RANS methods cannot accurately simulate unsteady separated wake flows, and running direct numerical simulation (DNS) on such complex flows is currently too computationally expensive. Large Eddy Simulation (LES) techniques allow for the computation of the large eddies, which contain most of the Reynolds stress, while modeling the smaller (subgrid) eddies. This results in models which are more computationally expensive than RANS methods, but not as prohibitive as DNS. By complimenting an LES approach with a RANS model, a hybrid LES/RANS method resolves the larger turbulent scales away from surfaces with LES, and switches to a RANS model inside boundary layers. As pointed out by Bertin et al., this type of hybrid approach has shown a lot of promise for predicting turbulent flows, but work is needed to verify that these models work well in hypersonic flows. The very limited amounts of flight and experimental data available presents an additional challenge for researchers. Recently, a joint study by NASA and CUBRC has focused on collecting heat transfer data
Investigation of Turbulent Tip Leakage Vortex in an Axial Water Jet Pump with Large Eddy Simulation
NASA Technical Reports Server (NTRS)
Hah, Chunill; Katz, Joseph
2012-01-01
Detailed steady and unsteady numerical studies were performed to investigate tip clearance flow in an axial water jet pump. The primary objective is to understand physics of unsteady tip clearance flow, unsteady tip leakage vortex, and cavitation inception in an axial water jet pump. Steady pressure field and resulting steady tip leakage vortex from a steady flow analysis do not seem to explain measured cavitation inception correctly. The measured flow field near the tip is unsteady and measured cavitation inception is highly transient. Flow visualization with cavitation bubbles shows that the leakage vortex is oscillating significantly and many intermittent vortex ropes are present between the suction side of the blade and the tip leakage core vortex. Although the flow field is highly transient, the overall flow structure is stable and a characteristic frequency seems to exist. To capture relevant flow physics as much as possible, a Reynolds-averaged Navier-Stokes (RANS) calculation and a Large Eddy Simulation (LES) were applied for the current investigation. The present study reveals that several vortices from the tip leakage vortex system cross the tip gap of the adjacent blade periodically. Sudden changes in local pressure field inside tip gap due to these vortices create vortex ropes. The instantaneous pressure filed inside the tip gap is drastically different from that of the steady flow simulation. Unsteady flow simulation which can calculate unsteady vortex motion is necessary to calculate cavitation inception accurately even at design flow condition in such a water jet pump.
Dynamic dose assessment by Large Eddy Simulation of the near-range atmospheric dispersion.
Vervecken, Lieven; Camps, Johan; Meyers, Johan
2015-03-01
In order to improve the simulation of the near-range atmospheric dispersion of radionuclides, computational fluid dynamics is becoming increasingly popular. In the current study, Large-Eddy Simulation is used to examine the time-evolution of the turbulent dispersion of radioactive gases in the atmospheric boundary layer, and it is coupled to a gamma dose rate model that is based on the point-kernel method with buildup factors. In this way, the variability of radiological dose rate from cloud shine due to instantaneous turbulent mixing processes can be evaluated. The steady release in an open field of (41)Ar and (133)Xe for 4 different release heights is studied, thus covering radionuclides that decay with a high-energy gamma and a low-energy gamma, respectively. Based on these simulations, the variability of dose rates at ground level for different averaging times in the dose measurements is analyzed. It is observed that turbulent variability in the wind field can lead to dose estimates that are underestimated by up to a factor of four when conventional long-term measurements are used to estimate the dose from short-term exposures. PMID:25634888
Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows
NASA Technical Reports Server (NTRS)
Givi, P.; Madnia, C. K.; Steinberger, C. J.; Frankel, S. H.; Vidoni, T. J.
1991-01-01
The main objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. In the efforts related to LES, we were concerned with developing reliable subgrid closures for modeling of the fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we focused our attention to further investigation of the effects of exothermicity in compressible turbulent flows. In our previous work, in the first year of this research, we have considered only 'simple' flows. Currently, we are in the process of extending our analyses for the purpose of modeling more practical flows of current interest at LaRC. A summary of our accomplishments during the third six months of the research is presented.
Large-Eddy Simulation of Flow Through an Array of Cubes with Local Grid Refinement
NASA Astrophysics Data System (ADS)
Goodfriend, Elijah; Katopodes Chow, Fotini; Vanella, Marcos; Balaras, Elias
2016-05-01
High resolution simulations of the transport of urban contaminants are important for disaster response and city planning. Large-eddy simulation (LES) and mesh refinement can each be used to decrease the computational cost of modelling, but combining these techniques can result in additional errors at grid-refinement interfaces. Here, we study the effect of the turbulence closure on the accuracy of LES results, for grids with mesh refinement, in a test case of flow through a periodic array of cubes. It is found that a mixed-model turbulence closure, using both an eddy viscosity and a scale similarity component, reduces energy accumulation at grid-refinement interfaces when used with explicit filtering of the advection term. The mixed model must be used with explicit filtering to control high wavenumber errors generated by the non-linear scale-similarity model. The results demonstrate that the turbulence closure mitigates errors associated with using LES on block-structured grids for urban-flow simulations.
Large-eddy simulation for the prediction of supersonic rectangular jet noise
NASA Astrophysics Data System (ADS)
Nichols, Joseph W.; Ham, Frank E.; Lele, Sanjiva K.; Bridges, James E.
2011-11-01
We investigate the noise from isothermal and heated under-expanded supersonic turbulent jets issuing from a rectangular nozzle of aspect ratio 4:1 using high-fidelity unstructured large-eddy simulation (LES) and acoustic projection based on the Ffowcs-Williams Hawkings (FWH) equations. The nozzle/flow interaction is directly included by simulating the flow in and around the nozzle in addition to the jet plume downstream. A grid resolution study is performed and results are shown for unstructured meshes containing up to 300 million control volumes, generated by a massively parallel code scaled to as many as 65,536 processors. Validated against laboratory measurements using a nozzle of precisely the same geometry, we find that mesh isotropy is a key factor in determining the quality of the far-field aeroacoustic predictions. The full flow fields produced by the simulation, in conjunction with particle image velocimetry (PIV) data measured from experiment, allow for a detailed examination of the interaction of large-scale coherent flow features and the resultant far-field noise, and its subsequent modification in the presence of heating. Supported by NASA grant NNX07AC94A and PSAAP, with computational resources from a DoD HPCMP CAP-2 project.
Large-eddy Simulation of Boundary Layer Flow over Desert Sand Dune Structures
NASA Astrophysics Data System (ADS)
Uhlrich, S.; Anderson, W.; Passalacqua, P.; Mohrig, D. C.; Kocurek, G.
2012-12-01
Complex spatiotemporal coupling exists between desert sand dune topography and surface layer physics of the atmospheric boundary layer (ABL). Although the interactions of individual desert sand dunes have been extensively studied, with categorical interaction mechanisms identified, the aero-mechanical coupling associated with these dune interactions remains an open problem. Large-eddy simulation (LES) is used to simulate turbulent boundary layer flow over dune structures from White Sands, NM. The dunes are resolved with an immersed boundary method (IBM). The flow-forcing (imposed pressure gradient) is varied to simulate the three common prevailing wind conditions at White Sands (southwest, southeast, and northwest, with southwest being the most common). In the present research, comparison between flow statistics (dune wall pressure distribution retrieved from the IBM) and time-difference dune elevation data are used to characterize the mechanisms responsible for erosion (stoss side) and deposition (lee side) of sand. Additionally, statistical details of time series of aerodynamic forcing at different locations on the dune face are evaluated, which may be used to deepen understanding of erosion and deposition events observed in the time-difference lidar data.
Wall modeled large eddy simulation of supersonic flow physics over compression-expansion ramp
NASA Astrophysics Data System (ADS)
Goshtasbi Rad, Ebrahim; Mousavi, Seyed Mahmood
2015-12-01
In the present work, wall modeled large-eddy simulation (WMLES) in the Fluent software is used to investigate the flow physics of a three-dimensional shock-turbulent boundary layer interaction, as an important phenomenon in aerospace science, on a compression-expansion ramp with the angle of 25°. Fine flow structures are obtained via Laplacian of density that called shadowgraph, in which shock wave structures are visible distinctly. The results are compared with the experimental data of Zheltovodov et al., 1990 [33], in the same condition regarding geometry, boundary conditions, etc. as those used by them. Results show that not only there are a good agreement with experimental trends concerning wall pressure, friction coefficient distribution and mean velocity profiles, but also in comparison with those presented by Grilli et al., 2013 [24]. LES simulation, used in this study, presents more accurate results with fewer computational costs. Afterwards, we investigated the influence of discontinuity in wall temperature, varying stagnation pressure and Reynolds number on physics of flow in order to control the shock behavior. Our simulations shows that, discontinuity in wall temperature, varying free stream stagnation pressure and Reynolds number (the free stream Mach number remained essentially constant) influences the starting point of shock, shock strength, separation length and the collision angle of separated and reattachment shock waves.
Large-eddy simulation of the generation and propagation of internal solitary waves
NASA Astrophysics Data System (ADS)
Zhu, Hai; Wang, LingLing; Tang, HongWu
2014-06-01
A modified large-eddy simulation model, the dynamic coherent eddy model (DCEM) is employed to simulate the generation and propagation of internal solitary waves (ISWs) of both depression and elevation type, with wave amplitudes ranging from small, medium to large scales. The simulation results agree well with the existing experimental data. The generation process of ISWs is successfully captured by the DCEM method. Shear instabilities and diapycnal mixing in the initial wave generation phase are observed. The dissipation rate is not equal at different locations of an ISW. ISW-induced velocity field is analyzed in the present study. The structure of the bottom boundary layer (BBL) of internal wave packets is found to be different from that of a single ISW. A reverse boundary jet instead of a separation bubble exists behind the leading internal wave while separation bubbles appear in other parts of the wave-induced velocity field. The boundary jet flow resulting from the adverse pressure gradients has distinctive dynamics compared with free shear jets.
High-resolution large-eddy simulation of turbulent mixing of a river plume
NASA Astrophysics Data System (ADS)
Yu, X.; Hsu, T. J.; Shi, F.; Kirby, J. T., Jr.
2014-12-01
A non-hydrostatic sigma-coordinate numerical model (NHWave) is applied to study the structure of a river plume, and the vertical mixing due to shear instabilities. A 3D large-eddy simulation approach is used with the aim to resolve the flow turbulence in the stratified ambient fluid at high Reynolds number. The domain is of depth 10m, length 500m and width 25m, and initially quiescent containing saltwater of salinity 26 psu. Fresh water plume is sent from the left boundary with a range of internal Froude number. Simulation resulting using Standard Smagorinsky closure demonstrates that the model is able to predict shear instabilities although it could not resolve the secondary instability at high Reynolds number. The characteristic length scale of the shear instabilities is around 10 m, which is consistent with field observation of Connecticut River plume using a 4-channel broadband echo sounder (Geyer et al. 2010, Geophy. Res. Lett., 37, L22607). The mixing efficiency and dissipation rate are obtained from the numerical simulation results, and these results are used to investigate and evaluate the Richardson-number-dependent parameterization of the mixing process. The model can also provide the information on fine structures of surface elevation variations, which enables us to correlate the surface signature with the turbulent billow underneath. The model therefore may be useful to help interpret surface signatures observed using various remote sensing techniques. Supported by Office of Naval Research.
On the large-eddy simulation of transitional wall-bounded flows
NASA Technical Reports Server (NTRS)
Piomelli, Ugo; Zang, Thomas A.; Speziale, Charles G.; Hussaini, M. Y.
1989-01-01
The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed.
Large-eddy simulation of vortex streets and dispersion behind high-rise buildings
NASA Astrophysics Data System (ADS)
Han, Beom-Soon; Park, Seung-Bu; Baik, Jong-Jin
2015-11-01
Understanding flow and dispersion in densely built-up urban areas is one of the important problems in the field of urban fluid mechanics. Nowadays, sophisticated numerical models and high-resolution urban morphology data enable us to study detailed flow structures in real urban areas. Simulations with high-resolution urban morphology data show very complex flow structures in several studies. Here, we examine turbulent flow patterns and associated pollutant dispersion near and, particularly, behind high-rise buildings using the parallelized large-eddy simulation model (PALM) and high-resolution urban morphology data. The study area selected is a highly built-up area of Seoul, South Korea. It is shown that turbulent wakes are produced behind high-rise buildings and vortex streets appear in the places where turbulent wakes occur. The vortex street seems to be related to strong updrafts and ejections that appear downwind of high-rise buildings. The vortex street is found to affect pollutant dispersion. Various factors that influence the evolution and structure of vortex streets will be presented and discussed along with involved dispersion mechanisms.
Large eddy simulation of wind-induced interunit dispersion around multistory buildings.
Ai, Z T; Mak, C M
2016-04-01
Previous studies regarding interunit dispersion used Reynolds-averaged Navier-Stokes (RANS) models and thus obtained only mean dispersion routes and re-entry ratios. Given that the envelope flow around a building is highly fluctuating, mean values could be insufficient to describe interunit dispersion. This study investigates the wind-induced interunit dispersion around multistory buildings using the large eddy simulation (LES) method. This is the first time interunit dispersion has been investigated transiently using a LES model. The quality of the selected LES model is seriously assured through both experimental validation and sensitivity analyses. Two aspects are paid special attention: (i) comparison of dispersion routes with those provided by previous RANS simulations and (ii) comparison of timescales with those of natural ventilation and the survival times of pathogens. The LES results reveal larger dispersion scopes than the RANS results. Such larger scopes could be caused by the fluctuating and stochastic nature of envelope flows, which, however, is canceled out by the inherent Reynolds-averaged treatment of RANS models. The timescales of interunit dispersion are comparable with those of natural ventilation. They are much shorter than the survival time of most pathogens under ordinary physical environments, indicating that interunit dispersion is a valid route for disease transmission. PMID:25787963
ENDLESS: An extended nonperiodic domain large-eddy simulation approach for scalar plumes
NASA Astrophysics Data System (ADS)
Chen, Bicheng; Yang, Di; Meneveau, Charles; Chamecki, Marcelo
2016-05-01
Large-eddy simulation (LES) has proven to be a valuable tool for high-fidelity modeling of environmental and geophysical turbulent flows. An important application of LES is to study the transport of effluents (e.g. oils from a subsea blowout) in the ocean mixed layer (OML). Oil plumes being transported in the OML experience the action of shear-generated turbulence, Langmuir circulations, Ekman transport and submesoscale quasi-geostrophic eddies. To resolve such turbulent processes, grid sizes of a few meters are desirable while horizontal domain sizes of LES are typically restricted from hundreds of meters to a few kilometers, for LES to remain practically affordable. Yet transported oil plumes evolve to large scales extending to tens or even hundreds of kilometers. In this study, the Extended Nonperiodic Domain LES for Scalar transport (ENDLESS) is proposed as a multi-scale approach to tackle this challenge while being computationally affordable. The basic idea is to simulate the shear turbulence and Langmuir circulations on a small horizontal domain with periodic boundary conditions while the resulting transport velocity field is replicated periodically following adaptively the large-scale plume as it evolves spatially towards much larger scales. This approach also permits the superposition of larger-scale quasi two-dimensional flow motions on the oil advection, allowing for coupling with regional circulation models. A validation case and two sample applications to oil plume evolution in the OML are presented in order to demonstrate key features and computational speedup associated with the ENDLESS method.
NASA Astrophysics Data System (ADS)
Schlegel, Fabian; Stiller, Jörg; Bienert, Anne; Maas, Hans-Gerd; Queck, Ronald; Bernhofer, Christian
2012-02-01
The effect of sub-tree forest heterogeneity in the flow past a clearing is investigated by means of large-eddy simulation (LES). For this purpose, a detailed representation of the canopy has been acquired by terrestrial laser scanning for a patch of approximately 190 m length in the field site "Tharandter Wald", near the city of Dresden, Germany. The scanning data are used to produce a high resolution plant area distribution (PAD) that is averaged over approximately one tree height (30 m) along the transverse direction, in order to simplify the LES study. Despite the smoothing involved with this procedure, the resulting two-dimensional PAD maintains a rich vertical and horizontal structure. For the LES study, the PAD is embedded in a larger domain covered with an idealized, horizontally homogeneous canopy. Simulations are performed for neutral conditions and compared to a LES with homogeneous PAD and recent field measurements. The results reveal a considerable influence of small-scale plant distribution on the mean velocity field as well as on turbulence data. Particularly near the edges of the clearing, where canopy structure is highly variable, usage of a realistic PAD appears to be crucial for capturing the local flow structure. Inside the forest, local variations in plant density induce a complex pattern of upward and downward motions, which remain visible in the mean flow and make it difficult to identify the "adjustment zone" behind the windward edge of the clearing.
Maurer, K. D.; Bohrer, G.; Kenny, W. T.; Ivanov, V. Y.
2015-04-30
Surface roughness parameters, namely the roughness length and displacement height, are an integral input used to model surface fluxes. However, most models assume these parameters to be a fixed property of plant functional type and disregard the governing structural heterogeneity and dynamics. In this study, we use large-eddy simulations to explore, in silico, the effects of canopy-structure characteristics on surface roughness parameters. We performed a virtual experiment to test the sensitivity of resolved surface roughness to four axes of canopy structure: (1) leaf area index, (2) the vertical profile of leaf density, (3) canopy height, and (4) canopy gap fraction.more » We found roughness parameters to be highly variable, but uncovered positive relationships between displacement height and maximum canopy height, aerodynamic canopy height and maximum canopy height and leaf area index, and eddy-penetration depth and gap fraction. We also found negative relationships between aerodynamic canopy height and gap fraction, as well as between eddy-penetration depth and maximum canopy height and leaf area index. We generalized our model results into a virtual "biometric" parameterization that relates roughness length and displacement height to canopy height, leaf area index, and gap fraction. Using a decade of wind and canopy-structure observations in a site in Michigan, we tested the effectiveness of our model-driven biometric parameterization approach in predicting the friction velocity over heterogeneous and disturbed canopies. We compared the accuracy of these predictions with the friction-velocity predictions obtained from the common simple approximation related to canopy height, the values calculated with large-eddy simulations of the explicit canopy structure as measured by airborne and ground-based lidar, two other parameterization approaches that utilize varying canopy-structure inputs, and the annual and decadal means of the surface roughness parameters at
Investigation of particle-laden flow in a straight duct using large eddy simulation
Fairweather, M.; Yao, J.
2007-07-01
A particle-laden turbulent flow in a square duct is predicted using large eddy simulation (LES). The simulation is performed for a Reynolds number of 35,500, and correctly predicts the existence of secondary flows and their effects on the mean flow. The results are also in good qualitative agreement with experimental data obtained at different Reynolds numbers. One-way coupling is assumed between solid particles and the fluid, and a particle equation of motion, including Stokes drag, lift, buoyancy and gravity force terms, solved using a Lagrangian particle tracking technique. Three sizes of particle (1, 50 and 100 {mu}m) are considered, and results demonstrate that size has a significant effect on particle dispersion and deposition in the duct flow. As particle size increases, therefore, they tend to settle on the floor of the duct, with less dispersion in the fluid phase. The study demonstrates the usefulness of LES for nuclear waste processing applications since secondary flows occur in many practically-relevant flows, and since it is desirable that the two-phase waste mixture is kept as homogeneous as possible to prevent, or at least discourage, the settling out of solid particles to form a bed which can promote pipe blockages. (authors)
Lantz, Jonas; Gårdhagen, Roland; Karlsson, Matts
2012-10-01
In this study, large-eddy simulation (LES) is employed to calculate the disturbed flow field and the wall shear stress (WSS) in a subject specific human aorta. Velocity and geometry measurements using magnetic resonance imaging (MRI) are taken as input to the model to provide accurate boundary conditions and to assure the physiological relevance. In total, 50 consecutive cardiac cycles were simulated from which a phase average was computed to get a statistically reliable result. A decomposition similar to Reynolds decomposition is introduced, where the WSS signal is divided into a pulsating part (due to the mass flow rate) and a fluctuating part (originating from the disturbed flow). Oscillatory shear index (OSI) is plotted against time-averaged WSS in a novel way, and locations on the aortic wall where elevated values existed could easily be found. In general, high and oscillating WSS values were found in the vicinity of the branches in the aortic arch, while low and oscillating WSS were present in the inner curvature of the descending aorta. The decomposition of WSS into a pulsating and a fluctuating part increases the understanding of how WSS affects the aortic wall, which enables both qualitative and quantitative comparisons. PMID:22209366
Reduction of Cloud Water in Ship Tracks: Observations and Large-Eddy Simulations
NASA Technical Reports Server (NTRS)
Ackerman, A. S.; Stevens, D. E.; Toon, O. B.; Coakley, J. A., Jr.; Gore, Warren J. (Technical Monitor)
2001-01-01
Ship tracks represent a natural laboratory to study the effects of aerosols on clouds. A number of observations and simulations have shown that increased droplet concentrations in ship tracks increase their total cross-sectional area, thereby enhancing cloud albedo and providing a negative radiative forcing at the surface and the top of the atmosphere. In some cases, cloud water has been found to be enhanced in ship tracks, which has been attributed to suppression of drizzle and implies an enhanced susceptibility of cloud albedo to droplet concentrations. However, more recently compiled observations indicate that cloud water is instead reduced in daytime ship tracks on average. Such a response is consistent with cloud-burning due to solar absorption by soot (the semi-direct radiative forcing of aerosols), recently suggested to be suppressing trade cumulus cloud coverage over the Indian Ocean. We will summarize observational evidence and present large-eddy simulations that consider these competing mechanisms in the effects of aerosols on cloud albedo.
From large-eddy simulation to multi-UAVs sampling of shallow cumulus clouds
NASA Astrophysics Data System (ADS)
Lamraoui, Fayçal; Roberts, Greg; Burnet, Frédéric
2016-04-01
In-situ sampling of clouds that can provide simultaneous measurements at satisfying spatio-temporal resolutions to capture 3D small scale physical processes continues to present challenges. This project (SKYSCANNER) aims at bringing together cloud sampling strategies using a swarm of unmanned aerial vehicles (UAVs) based on Large-eddy simulation (LES). The multi-UAV-based field campaigns with a personalized sampling strategy for individual clouds and cloud fields will significantly improve the understanding of the unresolved cloud physical processes. An extensive set of LES experiments for case studies from ARM-SGP site have been performed using MesoNH model at high resolutions down to 10 m. The carried out simulations led to establishing a macroscopic model that quantifies the interrelationship between micro- and macrophysical properties of shallow convective clouds. Both the geometry and evolution of individual clouds are critical to multi-UAV cloud sampling and path planning. The preliminary findings of the current project reveal several linear relationships that associate many cloud geometric parameters to cloud related meteorological variables. In addition, the horizontal wind speed indicates a proportional impact on cloud number concentration as well as triggering and prolonging the occurrence of cumulus clouds. In the framework of the joint collaboration that involves a Multidisciplinary Team (including institutes specializing in aviation, robotics and atmospheric science), this model will be a reference point for multi-UAVs sampling strategies and path planning.
Large-eddy simulation of airflow and heat transfer in a general ward of hospital
NASA Astrophysics Data System (ADS)
Hasan, Md. Farhad; Himika, Taasnim Ahmed; Molla, Md. Mamun
2016-07-01
In this paper, a very popular alternative computational technique, the Lattice Boltzmann Method (LBM) has been used for Large-Eddy Simulation (LES) of airflow and heat transfer in general ward of hospital. Different Reynolds numbers have been used to study the airflow pattern. In LES, Smagorinsky turbulence model has been considered and a discussion has been conducted in brief. A code validation has been performed comparing the present results with benchmark results for lid-driven cavity problem and the results are found to agree very well. LBM is demonstrated through simulation in forced convection inside hospital ward with six beds with a partition in the middle, which acted like a wall. Changes in average rate of heat transfer in terms of average Nusselt numbers have also been recorded in tabular format and necessary comparison has been showed. It was found that partition narrowed the path for airflow and once the air overcame this barrier, it got free space and turbulence appeared. For higher turbulence, the average rate of heat transfer increased and patients near the turbulence zone released maximum heat and felt more comfortable.
Investigations of Subsonic Compressible Boundary Layer Flows using Hybrid Large Eddy Simulations
NASA Astrophysics Data System (ADS)
Taylor, Sara Jo
The objective of this thesis is to investigate the spatially developing turbulent compressible boundary layer on a flat plate using the Spalart-Allmaras Detached Eddy Simulation (SA-DES) model [22] and the Nichols-Nelson hybrid Reynolds-Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) model [13] which have been implemented into the Wind-US 3.0 computational fluid dynamics code [30]; both of the hybrid approaches involve RANS modeling in the near-wall region and LES treatment in the outer region. Generation of unsteady turbulent inflow data is achieved via the prescribed energy spectrum method. The studies illustrated dependence on Reynolds number based on momentum thickness, Reθ, ranging from 3018 to 19430, and dependence on Mach number,
Large Eddy Simulation of Dilute Sediment Suspension in an Open Channel Flow
NASA Astrophysics Data System (ADS)
Agegnehu, Getnet; Smith, Heather D.
2012-11-01
Flow and suspended sediment transport in fully developed turbulent open channel flow has been investigated using Large Eddy Simulation. We used a three-dimensional, non-hydrostatic model, OpenFOAM for this study. Pre-evaluation of three existing turbulence closure schemes is performed by comparing the mean flow and turbulent quantities with the direct numerical simulation results of Moser et al. (1999). It is found that the Dynamic Mixed Smagorinsky model underestimates the wall shear stress compared to the Dynamic Smagorinsky and one equation Eddy Viscosity schemes. Moreover, the Dynamic Smagorinsky scheme gives relatively better results in both the mean and turbulent quantities. The advection-diffusion equation is solved for suspended sediment transport and the effect of sediment roughness is included in the momentum equation based on the rough wall formulation proposed by Cebeci and Bradshaw (1977). A pick up function based on van Rijn (1984) is used to determine the sediment erosion. The settling process is taken into account with a settling velocity appearing in the concentration equation. Sediment and flow quantities are validated by comparing with the experimental data of Lyn (1988). The coupled hydrodynamics results are in good agreement with the experimental data.
NASA Astrophysics Data System (ADS)
Liu, Zhongqiu; Li, Linmin; Li, Baokuan; Jiang, Maofa
2014-07-01
The current study developed a coupled computational model to simulate the transient fluid flow, solidification, and particle transport processes in a slab continuous-casting mold. Transient flow of molten steel in the mold is calculated using the large eddy simulation. An enthalpy-porosity approach is used for the analysis of solidification processes. The transport of bubble and non-metallic inclusion inside the liquid pool is calculated using the Lagrangian approach based on the transient flow field. A criterion of particle entrapment in the solidified shell is developed using the user-defined functions of FLUENT software (ANSYS, Inc., Canonsburg, PA). The predicted results of this model are compared with the measurements of the ultrasonic testing of the rolled steel plates and the water model experiments. The transient asymmetrical flow pattern inside the liquid pool exhibits quite satisfactory agreement with the corresponding measurements. The predicted complex instantaneous velocity field is composed of various small recirculation zones and multiple vortices. The transport of particles inside the liquid pool and the entrapment of particles in the solidified shell are not symmetric. The Magnus force can reduce the entrapment ratio of particles in the solidified shell, especially for smaller particles, but the effect is not obvious. The Marangoni force can play an important role in controlling the motion of particles, which increases the entrapment ratio of particles in the solidified shell obviously.
NASA Astrophysics Data System (ADS)
Gicquel, L. Y. M.; Staffelbach, G.; Sanjose, M.; Boileau, M.
2009-12-01
Being able to ignite or reignite a gas turbine engine in a cold and rarefied atmosphere is a critical issue for many aeronautical gas turbine manufacturers. From a fundamental point of view, the ignition of the first burner and the flame propagation from one burner to another are two phenomena that are usually not studied. The present work presents on-going and past Large Eddy Simulations (LES) on this specific subject and as investigated at CERFACS (European Centre for Research and Advanced Training in Scientific Computation) located in Toulouse, France. Validation steps and potential difficulties are underlined to ensure reliability of LES for such problems. Preliminary LES results on simple burners are then presented, followed by simulations of a complete ignition sequence in an annular helicopter chamber. For all cases and when possible, two-phase or purely gaseous LES have been applied to the experimentally simplified or the full geometries. For the latter, massively parallel computing (700 processors on a Cray XT3 machine) was essential to perform the computation. Results show that liquid fuel injection has a strong influence on the ignition times and the rate at which the flame progresses from burner to burner. The propagation speed characteristic of these phenomena is much higher than the turbulent flame speed. Based on an in-depth analysis of the computational data, the difference in speed is mainly identified as being due to thermal expansion and the flame speed is strongly modified by the main burner aerodynamics issued by the swirled injection.
NASA Astrophysics Data System (ADS)
Singh, Ramnik; Thomas, Brian G.; Vanka, Surya P.
2014-06-01
Transient flow during nominally steady conditions is responsible for many intermittent defects during the continuous casting of steel. The double-ruler electromagnetic field configuration, or "FC-Mold EMBr," is popular in commercial slab casting as it provides independent control of the applied static field near the jet and free surface regions of the mold. In the current study, transient flow in a typical commercial caster is simulated in the absence and in the presence of a double-ruler magnetic field, with rulers of equal strengths. Large eddy simulations with the in-house code CU-FLOW resolve the important transient behavior, using grids of over five million cells with a fast parallel solver. In the absence of a magnetic field, a double-roll pattern is observed, with transient unbalanced behavior, high surface velocities (~0.5 m/s), surface vortex formation, and very large surface-level fluctuations (~±12 mm). Applying the magnetic field suppresses the unbalanced behavior, producing a more complex mold flow pattern, but with much lower surface velocities (~0.1 m/s), and a flat surface level with small level fluctuations (<±1 mm). Nail board measurements taken at this commercial caster, in the absence of the field, matched reasonably well with the calculated results, both quantitatively and qualitatively.
Large eddy simulation of forced ignition of an annular bluff-body burner
Subramanian, V.; Domingo, P.; Vervisch, L.
2010-03-15
The optimization of the ignition process is a crucial issue in the design of many combustion systems. Large eddy simulation (LES) of a conical shaped bluff-body turbulent nonpremixed burner has been performed to study the impact of spark location on ignition success. This burner was experimentally investigated by Ahmed et al. [Combust. Flame 151 (2007) 366-385]. The present work focuses on the case without swirl, for which detailed measurements are available. First, cold-flow measurements of velocities and mixture fractions are compared with their LES counterparts, to assess the prediction capabilities of simulations in terms of flow and turbulent mixing. Time histories of velocities and mixture fractions are recorded at selected spots, to probe the resolved probability density function (pdf) of flow variables, in an attempt to reproduce, from the knowledge of LES-resolved instantaneous flow conditions, the experimentally observed reasons for success or failure of spark ignition. A flammability map is also constructed from the resolved mixture fraction pdf and compared with its experimental counterpart. LES of forced ignition is then performed using flamelet fully detailed tabulated chemistry combined with presumed pdfs. Various scenarios of flame kernel development are analyzed and correlated with typical flow conditions observed in this burner. The correlations between, velocities and mixture fraction values at the sparking time and the success or failure of ignition, are then further discussed and analyzed. (author)
Localized dynamic kinetic-energy-based models for stochastic coherent adaptive large eddy simulation
NASA Astrophysics Data System (ADS)
De Stefano, Giuliano; Vasilyev, Oleg V.; Goldstein, Daniel E.
2008-04-01
Stochastic coherent adaptive large eddy simulation (SCALES) is an extension of the large eddy simulation approach in which a wavelet filter-based dynamic grid adaptation strategy is employed to solve for the most "energetic" coherent structures in a turbulent field while modeling the effect of the less energetic background flow. In order to take full advantage of the ability of the method in simulating complex flows, the use of localized subgrid-scale models is required. In this paper, new local dynamic one-equation subgrid-scale models based on both eddy-viscosity and non-eddy-viscosity assumptions are proposed for SCALES. The models involve the definition of an additional field variable that represents the kinetic energy associated with the unresolved motions. This way, the energy transfer between resolved and residual flow structures is explicitly taken into account by the modeling procedure without an equilibrium assumption, as in the classical Smagorinsky approach. The wavelet-filtered incompressible Navier-Stokes equations for the velocity field, along with the additional evolution equation for the subgrid-scale kinetic energy variable, are numerically solved by means of the dynamically adaptive wavelet collocation solver. The proposed models are tested for freely decaying homogeneous turbulence at Reλ=72. It is shown that the SCALES results, obtained with less than 0.5% of the total nonadaptive computational nodes, closely match reference data from direct numerical simulation. In contrast to classical large eddy simulation, where the energetic small scales are poorly simulated, the agreement holds not only in terms of global statistical quantities but also in terms of spectral distribution of energy and, more importantly, enstrophy all the way down to the dissipative scales.
Large-Eddy Simulation of the Evolving Stable Boundary Layer Over Flat Terrain
Townsend, R
2002-01-02
The stable boundary layer (SBL) in the atmosphere is of considerable interest because it is often the worse case scenario for air pollution studies and health effect assessments associated with the accidental release of toxic material. Traditional modeling approaches used in such studies do not simulate the non-steady character of the velocity field, and hence often overpredict concentrations while underpredicting spatial coverage of potentially harmful concentrations of airborne material. The challenge for LES is to be able to resolve the rather small energy-containing eddies of the SBL while still maintaining an adequate domain size. This requires that the subgrid-scale (SGS) parameterization of turbulence incorporate an adequate representation of turbulent energy transfer. Recent studies have shown that both upscale and downscale energy transfer can occur simultaneously, but that overall the net transfer is downscale. Including the upscale transfer of turbulent energy (energy backscatter) is particularly important near the ground and under stably-stratified conditions. The goal of this research is to improve the ability to realistically simulate the SBL. The large-eddy simulation (LES) approach with its subgrid-scale (SGS) turbulence model does a better job of capturing the temporally and spatially varying features of the SBL than do Reynolds-averaging models. The scientific objectives of this research are: (1) to characterize features of the evolving SBL structure for a range of meteorological conditions (wind speed and surface cooling), (2) to simulate realistically the transfer of energy between resolved and subgrid scales, and (3) to apply results to improve simulation of dispersion in the SBL.
Large Eddy Simulation and Field Experiments of Pollen Transport in the Atmospheric Boundary Layer
NASA Astrophysics Data System (ADS)
Chamecki, M.; Meneveau, C.; Parlange, M. B.; van Hout, R.
2006-12-01
Dispersion of airborne pollen by the wind has been a subject of interest for botanists and allergists for a long time. More recently, the development of genetically modified crops and questions about cross-pollination and subsequent contamination of natural plant populations has brought even more interest to this field. A critical question is how far from the source field pollen grains will be advected. Clearly the answer depends on the aerodynamic properties of the pollen, geometrical properties of the field, topography, local vegetation, wind conditions, atmospheric stability, etc. As a consequence, field experiments are well suited to provide some information on pollen transport mechanisms but are limited to specific field and weather conditions. Numerical simulations do not have this drawback and can be a useful tool to study pollen dispersal in a variety of configurations. It is well known that the dispersion of particles in turbulent fields is strongly affected by the large scale coherent structures. Large Eddy Simulation (LES) is a technique that allows us to study the typical distances reached by pollen grains and, at the same time, resolve the larger coherent structures present in the atmospheric boundary layer. The main objective of this work is to simulate the dispersal of pollen grains in the atmospheric surface layer using LES. Pollen concentrations are simulated by an advection-diffusion equation including gravitational settling. Of extreme importance is the specification of the bottom boundary conditions characterizing the pollen source over the canopy and the deposition process everywhere else. In both cases we make use of the theoretical profile for suspended particles derived by Kind (1992). Field experiments were performed to study the applicability of the theoretical profile to pollen grains and the results are encouraging. Airborne concentrations as well as ground deposition from the simulations are compared to experimental data to validate the
NASA Astrophysics Data System (ADS)
Anupindi, Kameswararao; Delorme, Yann; Shetty, Dinesh A.; Frankel, Steven H.
2013-12-01
Computational fluid dynamics (CFD) simulations are becoming a reliable tool to understand hemodynamics, disease progression in pathological blood vessels and to predict medical device performance. Immersed boundary method (IBM) emerged as an attractive methodology because of its ability to efficiently handle complex moving and rotating geometries on structured grids. However, its application to study blood flow in complex, branching, patient-specific anatomies is scarce. This is because of the dominance of grid nodes in the exterior of the fluid domain over the useful grid nodes in the interior, rendering an inevitable memory and computational overhead. In order to alleviate this problem, we propose a novel multiblock based IBM that preserves the simplicity and effectiveness of the IBM on structured Cartesian meshes and enables handling of complex, anatomical geometries at a reduced memory overhead by minimizing the grid nodes in the exterior of the fluid domain. As pathological and medical device hemodynamics often involve complex, unsteady transitional or turbulent flow fields, a scale resolving turbulence model such as large eddy simulation (LES) is used in the present work. The proposed solver (here after referred as WenoHemo), is developed by enhancing an existing in-house high-order incompressible flow solver that was previously validated for its numerics and several LES models by Shetty et al. (2010) [33]. In the present work, WenoHemo is systematically validated for additional numerics introduced, such as IBM and the multiblock approach, by simulating laminar flow over a sphere and laminar flow over a backward facing step respectively. Then, we validate the entire solver methodology by simulating laminar and transitional flow in abdominal aortic aneurysm (AAA). Finally, we perform blood flow simulations in the challenging clinically relevant thoracic aortic aneurysm (TAA), to gain insights into the type of fluid flow patterns that exist in pathological
Investigation of Rossby-number similarity in the neutral boundary layer using large-eddy simulation
Ohmstede, W.D.; Cederwall, R.T.; Meyers, R.E.
1988-01-01
One special case of particular interest, especially to theoreticians, is the steady-state, horizontally homogeneous, autobarotropic (PLB), hereafter referred to as the neutral boundary layer (NBL). The NBL is in fact a 'rare' atmospheric phenomenon, generally associated with high-wind situations. Nevertheless, there is a disproportionate interest in this problem because Rossby-number similarity theory provides a sound approach for addressing this issue. Rossby-number similarity theory has rather wide acceptance, but because of the rarity of the 'true' NBL state, there remains an inadequate experimental database for quantifying constants associated with the Rossby-number similarity concept. Although it remains a controversial issue, it has been proposed that large-eddy simulation (LES) is an alternative to physical experimentation for obtaining basic atmospherc 'data'. The objective of the study reported here is to investigate Rossby-number similarity in the NBL using LES. Previous studies have not addressed Rossby-number similarity explicitly, although they made use of it in the interpretation of their results. The intent is to calculate several sets of NBL solutions that are ambiguous relative to the their respective Rossby numbers and compare the results for similarity, or the lack of it. 14 refs., 1 fig.
Large-eddy simulation of cavitating nozzle flow and primary jet break-up
Örley, F. Trummler, T.; Mihatsch, M. S.; Schmidt, S. J.; Adams, N. A.; Hickel, S.
2015-08-15
We employ a barotropic two-phase/two-fluid model to study the primary break-up of cavitating liquid jets emanating from a rectangular nozzle, which resembles a high aspect-ratio slot flow. All components (i.e., gas, liquid, and vapor) are represented by a homogeneous mixture approach. The cavitating fluid model is based on a thermodynamic-equilibrium assumption. Compressibility of all phases enables full resolution of collapse-induced pressure wave dynamics. The thermodynamic model is embedded into an implicit large-eddy simulation (LES) environment. The considered configuration follows the general setup of a reference experiment and is a generic reproduction of a scaled-up fuel injector or control valve as found in an automotive engine. Due to the experimental conditions, it operates, however, at significantly lower pressures. LES results are compared to the experimental reference for validation. Three different operating points are studied, which differ in terms of the development of cavitation regions and the jet break-up characteristics. Observed differences between experimental and numerical data in some of the investigated cases can be caused by uncertainties in meeting nominal parameters by the experiment. The investigation reveals that three main mechanisms promote primary jet break-up: collapse-induced turbulent fluctuations near the outlet, entrainment of free gas into the nozzle, and collapse events inside the jet near the liquid-gas interface.
Large-eddy simulation of flow through a plane, asymmetric diffuser
NASA Technical Reports Server (NTRS)
Kaltenbach, Hans-Jakob
1994-01-01
A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.
Chernyshov, Alexander A.; Karelsky, Kirill V.; Petrosyan, Arakel S.
2010-06-16
Using advantages of large eddy simulation method, nontrivial regime of compressible magnetohydrodynamic turbulence of space plasma when initially supersonic fluctuations become weakly compressible is studied. Establishment of weakly compressible limit with Kolmogorov-like density fluctuations spectrum is shown in present work. We use our computations results to study dynamics of the turbulent plasma beta and anisotropic properties of the magnetoplasma fluctuations in the local interstellar medium. An outstanding, as yet unexplained, observation is that density fluctuations in the local interstellar medium exhibit a Kolmogorov-like spectrum over an extraordinary range of scales with a spectral index close to -5/3. In spite of the compressibility and the presence of magnetic field in the local interstellar medium, density fluctuations nevertheless admit a Kolmogorov-like power law. Supersonic flows with high value of large-scale Mach numbers are characterized in interstellar medium, nevertheless, there are subsonic fluctuations of weakly compressible components of interstellar medium. These weakly compressible subsonic fluctuations are responsible for emergence of a Kolmogorov-type spectrum in interstellar turbulence which is observed from experimental data. It is shown that density fluctuations are a passive scalar in a velocity field in weakly compressible magnetohydrodynamic turbulence and demonstrate Kolmogorov-like spectrum.
Recent advances in large-eddy simulation of spray and coal combustion
NASA Astrophysics Data System (ADS)
Zhou, L. X.
2013-07-01
Large-eddy simulation (LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering. Spray and coal combustion is widely used in power, transportation, chemical and metallurgical, iron and steel making, aeronautical and astronautical engineering, hence LES of spray and coal two-phase combustion is particularly important for engineering application. LES of two-phase combustion attracts more and more attention; since it can give the detailed instantaneous flow and flame structures and more exact statistical results than those given by the Reynolds averaged modeling (RANS modeling). One of the key problems in LES is to develop sub-grid scale (SGS) models, including SGS stress models and combustion models. Different investigators proposed or adopted various SGS models. In this paper the present author attempts to review the advances in studies on LES of spray and coal combustion, including the studies done by the present author and his colleagues. Different SGS models adopted by different investigators are described, some of their main results are summarized, and finally some research needs are discussed.
Large-Eddy Simulation of Shock-Wave Boundary Layer Interaction and its Control Using Sparkjet
NASA Astrophysics Data System (ADS)
Yang, Guang; Yao, Yufeng; Fang, Jian; Gan, Tian; Lu, Lipeng
2016-06-01
Large-eddy simulation (LES) of an oblique shock-wave generated by an 8° sharp wedge impinging onto a spatially-developing Mach 2.3 turbulent boundary layer and their interactions has been carried out in this study. The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is 20,000. The detailed numerical approaches are described and the inflow turbulence is generated using the digital filter method to avoid artificial temporal or streamwise periodicity. Numerical results are compared with the available wind tunnel PIV measurements of the same flow conditions. Further LES study on the control of flow separation due to the strong shock-viscous interaction is also conducted by using an active control actuator “SparkJet” concept. The single-pulsed characteristics of the control device are obtained and compared with the experiments. Instantaneous flowfield shows that the “SparkJet” promotes the flow mixing in the boundary layer and enhances its ability to resist the flow separation. The time and spanwise averaged skin friction coefficient distribution demonstrates that the separation bubble length is reduced by maximum 35% with the control exerted.
NASA Astrophysics Data System (ADS)
Nouri Gheimassi, Arash; Givi, Peyman; Nik, Mehdi B.; Pope, Stephen B.
2015-11-01
A new model is developed which accounts for the effects of subgrid scale pressure in the context of the filtered density function (FDF) formulation. This results in a pressure-velocity-scalar filtered mass density function (PVS-FMDF), which is suitable for large eddy simulation of compressible turbulence. Following its mathematical definition, an exact transport equation is derived for the PVS-FMDF. This equation is modeled in a probabilistic manner by a system of stochastic differential equations (SDEs). The consistency and the predictive capability of the model are established by conducting LES of a three-dimensional compressible mixing layer, and comparison with direct numerical simulation (DNS) data.
Dynamic optimization methodology based on subgrid-scale dissipation for large eddy simulation
NASA Astrophysics Data System (ADS)
Yu, Changping; Xiao, Zuoli; Li, Xinliang
2016-01-01
A dynamic procedure based on subgrid-scale dissipation is proposed for large eddy simulation of turbulent flows. In the new method, the model coefficients are determined by minimizing the square error of the resolved dissipation rate based on the Germano identity. A dynamic two-term mixed model is tested and evaluated both a priori and a posteriori in simulations of homogeneous and isotropic turbulence. The new dynamic procedure proves to be more effective to optimize the model coefficients as compared with traditional method. The corresponding dynamic mixed model can predict the physical quantities more accurately than traditional dynamic mixed model.