scholarly journals Finite Element Analysis of Turbulent Flows Using LES and Dynamic Subgrid-Scale Models in Complex Geometries

2011 ◽  
Vol 2011 ◽  
pp. 1-20 ◽  
Author(s):  
Wang Wenquan ◽  
Zhang Lixiang ◽  
Yan Yan ◽  
Guo Yakun

An innovative computational model is presented for the large eddy simulation (LES) of multidimensional unsteady turbulent flow problems in complex geometries. The main objectives of this research are to know more about the structure of turbulent flows, to identify their three-dimensional characteristic, and to study physical effects due to complex fluid flow. The filtered Navier-Stokes equations are used to simulate large scales; however, they are supplemented by dynamic subgrid-scale (DSGS) models to simulate the energy transfer from large scales toward subgrid-scales, where this energy will be dissipated by molecular viscosity. Based on the Taylor-Galerkin schemes for the convection-diffusion problems, this model is implemented in a three-dimensional finite element code using a three-step finite element method (FEM). Turbulent channel flow and flow over a backward-facing step are considered as a benchmark for validating the methodology by comparing with the direct numerical simulation (DNS) results or experimental data. Also, qualitative and quantitative aspects of three-dimensional complex turbulent flow in a strong 3D blade passage of a Francis turbine are analyzed.

1987 ◽  
Vol 109 (1) ◽  
pp. 71-76 ◽  
Author(s):  
J. O. Medwell ◽  
D. T. Gethin ◽  
C. Taylor

The performance of a cylindrical bore bearing fed by two axial grooves orthogonal to the load line is analyzed by solving the Navier-Stokes equations using the finite element method. This produces detailed information about the three-dimensional velocity and pressure field within the hydrodynamic film. It is also shown that the method may be applied to long bearing geometries where recirculatory flows occur and in which the governing equations are elliptic. As expected the analysis confirms that lubricant inertia does not affect bearing performance significantly.


Author(s):  
Jose´ Manuel Franco-Nava ◽  
Oscar Dorantes-Go´mez ◽  
Erik Rosado-Tamariz ◽  
Jose´ Manuel Ferna´ndez-Da´vila ◽  
Reynaldo Rangel-Espinosa

The stress analysis of the runner due to different loading is one of the most important tools that contribute its structural integrity evaluation. Finite element method has shown to be a strong numerical technique to provide good engineering accuracy. In this paper, the flow induced stresses in a Francis turbine runner is presented by using the finite element analysis. The runner geometry considered within the computational domain was modelled by using a three-dimensional laser triangulation scanner coupled with a portable coordinate measurement system. The runner geometry was generated by a number of 3D sub models, one for each of the main components of the runner, crown, band and a blade. In order to obtain a blade geometry a portable coordinate measurement system based on optical digitalization technology (scanner technology) was used. Because of symmetry, only a section of the runner domain was used for the finite element analysis. The runner was modeled with twenty-node solid elements. Loads due to pressure on the blade were derived from CFD computations for the runner at different power conditions (100%, 85% and 75%) for a medium head hydro power plant. CFD computations were carried out using the Finite Volume Method implement within FINE™/Turbo by NUMECA. The turbulence mathematical model used for the CFD computation was the Sparlart-Allmaras. The mesh of the turbine runner included different computational domains. For the runner blades the computational domain (mesh block) was defined in order to capture the complete blade row. All mesh blocks were structured hexahedral. Centrifugal force based on the rotational speed was considered. Also, a combined type loading analysis was computed including both pressure and centrifugal force. Appropriate boundary conditions were set in order to obtain the results due to the different type of analysis. The number of finite elements included in the FEM model was able to capture the pressure gradients on the blade surfaces obtained from the CFD results, which were investigated by application of a three dimensional Navier-Stoke commercial turbomachinery oriented CFD code. Analysis of the flow through the spiral case and stay vanes was carried out so as to include appropriate flow effects induced by these components and boundary conditions at the inlet of the wicket. A CFD analysis for the wicket and runner was carried out to generate the so called CFD reference solution. The analysis presented in this paper represents an initial characterization in order to increase understanding about combined loads acting on blades and to establish a reference state of stresses further comparison after refurbishments or optimization of the runner blades for a medium head hydroelectric power station.


Author(s):  
Leo A. Carrilho ◽  
Jamil Khan ◽  
Michael E. Conner ◽  
Abdel Mandour ◽  
Milorad B. Dzodzo

The effects of artificial roughness for the purpose of thermal performance improvement in pressurized water nuclear reactors are investigated. The artificial roughness consists of two-dimensional ribs parallel to the turbulent flow. The fuel rod bundle subchannel is preliminarily modeled as an annulus using the finite element method in ANSYS/FLOTRAN. The Navier-Stokes equations are solved from the SST (Shear Stress Transport) turbulence model for the simulated annulus thermal-flow. The analyses are performed for ribs dimensions and pitch provided by published previous work. It is found that, heat transfer and differential pressure have similar behavior with highest heat transfer occurring at the reattachment point. The finite element model describes well the characteristics of turbulent flow in smooth and rough rod when compared to previous semi-empirical models. Next paper extends the analysis by comparing numerical results with experimental test data and sensitivity analyses for different roughness configurations.


2005 ◽  
Vol 127 (3) ◽  
pp. 467-483 ◽  
Author(s):  
Ramin K. Rahmani ◽  
Theo G. Keith ◽  
Anahita Ayasoufi

In many branches of processing industries, viscous liquids need to be homogenized in continuous operations. Consequently, fluid mixing plays a critical role in the success or failure of these processes. Static mixers have been utilized over a wide range of applications such as continuous mixing, blending, heat and mass transfer processes, chemical reactions, etc. This paper describes how static mixing processes of single-phase viscous liquids can be simulated numerically, presents the flow pattern through a helical static mixer, and provides useful information that can be extracted from the simulation results. The three-dimensional finite volume computational fluid dynamics code used here solves the Navier-Stokes equations for both laminar and turbulent flow cases. The turbulent flow cases were solved using k-ω model and Reynolds stress model (RSM). The flow properties are calculated and the static mixer performance for different Reynolds numbers (from creeping flows to turbulent flows) is studied. A new parameter is introduced to measure the degree of mixing quantitatively. Furthermore, the results obtained by k-ω and RSM turbulence models and various numerical details of each model are compared. The calculated pressure drop is in good agreement with existing experimental data.


1984 ◽  
Vol 148 ◽  
pp. 413-442 ◽  
Author(s):  
S. Biringen

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 1500. Three-dimensional incompressible Navier–Stokes equations are numerically integrated to obtain the time evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 × 51 × 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and characteristic of the various stages of transition that lead to final breakdown. In particular, evidence points to the formation of a A-shaped vortex and the subsequent system of horsehoe vortices inclined to the main flow direction as the primary elements of transition. Details of the resulting flow field after breakdown indicate the evolution of streaklike formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), the introduction of subgrid-scale terms seems necessary to obtain fully developed turbulence statistics.


Author(s):  
Hiroshi Dan ◽  
◽  
Hiroshi Kanayama

In this paper, viscous shallow-water equations are derived from three-dimensional Navier-Stokes equations under the hydrostatic assumption. The viscous shallow-water equations are approximated by the finite element method based on our numerical scheme developed in 1978. This approach is used to simulate a tsunami in Hakata Bay. Results show a reasonable estimate of the tsunami arrival time.


Author(s):  
Raza A. Saeed

This paper presents the results of modelling of the complete three-dimensional fluid flow through the spiral casing, stay vanes, guide vanes, and then through the Francis turbine runner to the draft tube of the Derbendikhan power station. To investigate the flow in the Francis turbine and also to compute stress distribution in the runner blades, a three-dimensional model was prepared according to specifications provided. The two topics discussed in this study are: (i) the simulation of the 3D fluid flow through the inter blade channels for the Francis turbine runner by using Computational Fluid Dynamics (CFD) and, (ii) the simulation of the stress analysis of the turbine runner by using Finite Element Analysis (FEA). In this study, the water pressure obtained from the CFD analysis for different boundary conditions are incorporated into a Finite Element model to calculate stress distributions in the runner.


Author(s):  
Nurullah Türker ◽  
Hümeyra Tercanlı Alkış ◽  
Steven J Sadowsky ◽  
Ulviye Şebnem Büyükkaplan

An ideal occlusal scheme plays an important role in a good prognosis of All-on-Four applications, as it does for other implant therapies, due to the potential impact of occlusal loads on implant prosthetic components. The aim of the present three-dimensional (3D) finite element analysis (FEA) study was to investigate the stresses on abutments, screws and prostheses that are generated by occlusal loads via different occlusal schemes in the All-on-Four concept. Three-dimensional models of the maxilla, mandible, implants, implant substructures and prostheses were designed according to the All-on-Four concept. Forces were applied from the occlusal contact points formed in maximum intercuspation and eccentric movements in canine guidance occlusion (CGO), group function occlusion (GFO) and lingualized occlusion (LO). The von Mises stress values for abutment and screws and deformation values for prostheses were obtained and results were evaluated comparatively. It was observed that the stresses on screws and abutments were more evenly distributed in GFO. Maximum deformation values for prosthesis were observed in the CFO model for lateral movement both in the maxilla and mandible. Within the limits of the present study, GFO may be suggested to reduce stresses on screws, abutments and prostheses in the All-on-Four concept.


1990 ◽  
Vol 18 (4) ◽  
pp. 216-235 ◽  
Author(s):  
J. De Eskinazi ◽  
K. Ishihara ◽  
H. Volk ◽  
T. C. Warholic

Abstract The paper describes the intention of the authors to determine whether it is possible to predict relative belt edge endurance for radial passenger car tires using the finite element method. Three groups of tires with different belt edge configurations were tested on a fleet test in an attempt to validate predictions from the finite element results. A two-dimensional, axisymmetric finite element analysis was first used to determine if the results from such an analysis, with emphasis on the shear deformations between the belts, could be used to predict a relative ranking for belt edge endurance. It is shown that such an analysis can lead to erroneous conclusions. A three-dimensional analysis in which tires are modeled under free rotation and static vertical loading was performed next. This approach resulted in an improvement in the quality of the correlations. The differences in the predicted values of various stress analysis parameters for the three belt edge configurations are studied and their implication on predicting belt edge endurance is discussed.


Sign in / Sign up

Export Citation Format

Share Document