Computation of complex turbulent flow using matrix-free implicit dual time-stepping scheme and LRN turbulence model on unstructured grids

2004 ◽  
Vol 33 (1) ◽  
pp. 119-136 ◽  
Author(s):  
Yong Zhao
2009 ◽  
Vol 23 (03) ◽  
pp. 325-328 ◽  
Author(s):  
SHULING TIAN ◽  
YIZHAO WU ◽  
JIAN XIA

A parallel Navier-Stokes solver based on dynamic overset unstructured grids method is presented to simulate the unsteady turbulent flow field around helicopter in forward flight. The grid method has the advantages of unstructured grid and Chimera grid and is suitable to deal with multiple bodies in relatively moving. Unsteady Navier-Stokes equations are solved on overset unstructured grids by an explicit dual time-stepping, finite volume method. Preconditioning method applied to inner iteration of the dual-time stepping is used to speed up the convergence of numerical simulation. The Spalart-Allmaras one-equation turbulence model is used to evaluate the turbulent viscosity. Parallel computation is based on the dynamic domain decomposition method in overset unstructured grids system at each physical time step. A generic helicopter Robin with a four-blade rotor in forward flight is considered to validate the method presented in this paper. Numerical simulation results show that the parallel dynamic overset unstructured grids method is very efficient for the simulation of helicopter flow field and the results are reliable.


AIAA Journal ◽  
1997 ◽  
Vol 35 ◽  
pp. 1548-1550 ◽  
Author(s):  
S. DeRango ◽  
D. W. Zingg

Volume 4 ◽  
2004 ◽  
Author(s):  
Branislav Basara ◽  
Ales Alajbegovic ◽  
Decan Beader

The paper presents calculations of flow in a mixing vessel stirred by a six-blade Rushton impeller. Mathematical model used in computations is based on the ensemble averaged conservation equations. An efficient finite-volume method based on unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements is used together with various turbulence models. Besides the standard k-ε model which served as a reference, k-ε-v2 model (Durbin, 1995) and the recently proposed hybrid EVM/RSM turbulence model (Basara & Jakirlic, 2003) were used in the calculations. The main aim of the paper is to investigate if more advanced turbulence models are needed for this type of CFD applications. The results are compared with the available experimental data.


Author(s):  
Loi¨c Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-Harmonic Balance Method (mHBM). The Dynamic Lagrangian Frequency-Time method is used to calculate contact forces in the frequency domain. A new strategy for solving non-linear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.


2007 ◽  
Vol 11 (4) ◽  
pp. 171-178
Author(s):  
Khalid Alammar

Using the standard k-e turbulence model, an incompressible, axisymmetric turbulent flow with a sudden expansion was simulated. Effect of Prandtl number on heat transfer characteristics downstream of the expansion was investigated. The simulation revealed circulation downstream of the expansion. A secondary circulation (corner eddy) was also predicted. Reattachment was predicted at approximately 10 step heights. Corresponding to Prandtl number of 7.0, a peak Nusselt number 13 times the fully-developed value was predicted. The ratio of peak to fully-developed Nusselt number was shown to decrease with decreasing Prandtl number. Location of maximum Nusselt number was insensitive to Prandtl number.


Sign in / Sign up

Export Citation Format

Share Document