Application of a New Hybrid Explicit-Implicit Flow Solver to 1D Unsteady Flows with Shock Waves

Author(s):  
E. Timofeev ◽  
F. Norouzi
Author(s):  
C B Allen

A grid adaptation procedure suitable for use during unsteady flow computations is described. Transfinite interpolation is used to generate structured grids for the computation of steady and unsteady Euler flows past aerofoils. This technique is well suited to unsteady flows, since instantaneous grid positions and speeds required by the flow solver are available directly from the algebraic mapping. A different approach to grid adaptation is described, wherein adaptation is performed by redistributing the interpolation parameters, instead of the physical grid positions. This results in the adapted grid positions, and hence speeds, still being available algebraically. Grid adaptation during an unsteady computation is performed continuously by imposing an ‘adaptation velocity’ on grid points, thereby applying the adaptation over several time steps and avoiding the interpolation of the solution from one grid to another, which is associated with instantaneous adaptation. For both steady and unsteady flows the adapted grid technique is shown to produce sharper shock resolution for a very small increase in CPU (central processing unit) requirements.


2003 ◽  
Vol 125 (1) ◽  
pp. 25-32 ◽  
Author(s):  
W. Ning ◽  
Y. S. Li ◽  
R. G. Wells

A multistage frequency domain (time-linearized/nonlinear harmonic) Navier-Stokes unsteady flow solver has been developed for predicting unsteady flows induced by bladerow interactions. In this paper, the time-linearized option of the solver has been used to analyze unsteady flows in a subsonic turbine test stage and the DLR transonic counter-rotating shrouded propfan. The numerical accuracy and computational efficiency of the time-linearized viscous methods have been demonstrated by comparing predictions with test data and nonlinear time-marching solutions for these two test cases. It is concluded that the development of efficient frequency domain approaches enables unsteady flow predictions to be used in the design cycles to tackle aeromechanics problems.


2000 ◽  
Vol 123 (3) ◽  
pp. 680-685 ◽  
Author(s):  
L. He ◽  
K. Sato

A three-dimensional incompressible viscous flow solver of the thin-layer Navier-Stokes equations was developed for the unsteady turbomachinery flow computations. The solution algorithm for the unsteady flows combines the dual time stepping technique with the artificial compressibility approach for solving the incompressible unsteady flow governing equations. For time accurate calculations, subiterations are introduced by marching the equations in the pseudo-time to fully recover the incompressible continuity equation at each real time step, accelerated with a multi-grid technique. Computations of test cases show satisfactory agreements with corresponding theoretical and experimental results, demonstrating the validity and applicability of the present method to unsteady incompressible turbomachinery flows.


Author(s):  
Pengcheng Du ◽  
Fangfei Ning

Time periodic unsteady flows are often encountered in turbomachinery. Simulating such flows using conventional time marching approach is very time-consuming and hence expensive. To handle this problem, several Fourier-based reduced order models have been developed recently. Among these, the time-domain harmonic balance method solves the governing equations purely in the time domain and there is also no need for the turbulence model to be linearized, making it easy to be implemented in an existing RANS code. Thus, the time-domain harmonic balance method was chosen and incorporated into an in-house Navier-Stokes flow solver. Several test cases were performed for the validations of the developed code. They cover standard unsteady test cases such as the low speed vortex shedding cylinder flow and the Sajben transonic diffuser under periodically oscillating back pressure. Further, two different practical turbomachinery unsteady flows were considered. One is a transonic fan under circumferential inlet distortion and the other is the rotor-stator interactions in a single stage compressor. The results illustrate the capability of the harmonic balance method in capturing the dominant nonlinear effects. The number of harmonics should be retained in the harmonic balance method is depend on the strength of the nonlinear unsteady effects and differs from case to case. With appropriate number of harmonics retained, it can resolve the unsteady flow field satisfactory, meanwhile, reducing the computational time significantly. In a word, the harmonic balance method promise to be an effective way to simulate time periodic unsteady flows.


Author(s):  
W. Ning ◽  
Y. S. Li ◽  
R. G. Wells

A multistage frequency domain (time-linearized/nonlinear harmonic) Navier-Stokes unsteady flow solver has been developed for predicting unsteady flows induced by bladerow interactions. In this paper, the time-linearized option of the solver has been used to analyze unsteady flows in a subsonic turbine test stage and the DLR transonic counter-rotating shrouded propfan. The numerical accuracy and computational efficiency of the time-linearized viscous methods have been demonstrated by comparing predictions with test data and nonlinear time-marching solutions for these two test cases. It is concluded that the development of efficient frequency domain approaches enables unsteady flow predictions to be used in the design cycles to tackle aeromechanics problems.


Author(s):  
Mohagna Pandya ◽  
Neal Frink ◽  
Khaled Abdol-Hamid ◽  
James Chung
Keyword(s):  

2008 ◽  
Vol 9 (2) ◽  
Author(s):  
R.A.B. Queiroz ◽  
V.G. Ferreira ◽  
R.G. Cuenca

Author(s):  
В.И. Пинчуков

Работа посвящена поиску новых автоколебательных течений и их численному изучению. Эти поиски проводились путем расчетных исследований течений с максимальным количеством контактных разрывов и точек пересечения разрывов ударных волн с ударными волнами или ударных волн с контактными разрывами. Рассмотрены два класса течений: натекание звуковой недорасширенной струи на цилиндрическое тело, расположенное в открытой трубе, и взаимодействие сверхзвукового однородного потока с системой открытый канал вращения (с переменным сечением) цилиндрическое тело на оси. В обоих случаях найдены автоколебательные режимы. Двумерные осесимметричные уравнения сжимаемого газа решены с помощью неявной схемы РунгеКутты третьего порядка. Применяется алгебраическая турбулентная вязкость, основанная на использовании обобщенной формулы Кармана. This paper addresses a search for new selfoscillatory compressible flows and numerical studies of these flows. These searches are carried out by computational modelling of currents with the maximum number of contact discontinuities and points of intersection of discontinuities shock waves with shock waves or shock waves with contact discontinuities. Two families of unsteady flows are considered. The first one contains flows near underextended sonic jets, impinging on cylindrical bodies placed in open tubes. The second family corresponds to interactions of uniform supersonic streams with pairs containing the open channel of rotation (with transient crossection) and a cylindrical body on the axis. Selfoscillatory regimes are found in both cases. Twodimensional axysimmetrical compressble flow equations are solved by an implicit RungeKutta scheme of the third order. Algebraic turbulent viscosity is assumed which is based on the implementation of the generalized Karman formulae. Numerical results allow concluding that unsteady flows, which take place when sonic jets impinge on a pair containing of cylinders and open tubes are typical for jets impinging on obstacles. Flows, which take place when uniform streams interact with these pairs comprise a new original class of selfoscillatory flows.


Author(s):  
M.A. Mogilevsky ◽  
L.S. Bushnev

Single crystals of Al were loaded by 15 to 40 GPa shock waves at 77 K with a pulse duration of 1.0 to 0.5 μs and a residual deformation of ∼1%. The analysis of deformation structure peculiarities allows the deformation history to be re-established.After a 20 to 40 GPa loading the dislocation density in the recovered samples was about 1010 cm-2. By measuring the thickness of the 40 GPa shock front in Al, a plastic deformation velocity of 1.07 x 108 s-1 is obtained, from where the moving dislocation density at the front is 7 x 1010 cm-2. A very small part of dislocations moves during the whole time of compression, i.e. a total dislocation density at the front must be in excess of this value by one or two orders. Consequently, due to extremely high stresses, at the front there exists a very unstable structure which is rearranged later with a noticeable decrease in dislocation density.


Sign in / Sign up

Export Citation Format

Share Document