scholarly journals Analysis of Aspect Ratio Effects on Transonic Rotor Performance Using a 3D Viscous Solver

1998 ◽  
Author(s):  
Cecil R. Buchanan ◽  
Paul R. Emmerson ◽  
Michael Spruce
Keyword(s):  
Aspect Ratio ◽  
Numerical Study ◽  
Navier Stokes ◽  
Test Cases ◽  
Compressor Rotor ◽  
Aspect Ratios ◽  
Radial Profiles ◽  
Ratio Change ◽  
The One ◽  
Good Agreement

This paper presents the results of a numerical study into the effects of aspect ratio on compressor rotor performance. The test cases studied are NASA rotors 37 and 38, which have aspect ratios of 1.19 and 1.63 respectively. A 3D, single-passage steady flow Navier-Stokes solver was used to predict complete performance characteristics, including the numerical instability point, for both rotors. The predictions are generally in good agreement with the test data (characteristics, radial profiles and rotor over tip measurements) at all conditions modelled for rotor 37. The performance for rotor 38 is overpredicted, with slightly less than half of the measured performance difference between the two rotors being captured. The effect of a pure aspect ratio change (divorced from the rotor inlet to exit area changes present in the rotor 37/38 comparisons) was also investigated, and a case with an aspect ratio double that of rotor 37 was also modelled. The results indicated that the code predicted little effect on rotor performance due to an aspect ratio change alone (from 1.19 to 2.38). This is surprising and it raises doubts about the ability of current codes (or at least the one used in the study) to predict this important aspect of a compressor design adequately.

Numerical Prediction of Turbulent Wakes Behind a Square Cylinder

1998 ◽  
Vol 14 (1) ◽  
pp. 23-29
Author(s):  
Robert R. Hwang ◽  
Sheng-Yuh Jaw
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Navier Stokes ◽  
Wall Region ◽  
Square Cylinder ◽  
Mean Velocity ◽  
Navier Stokes Equations ◽  
Turbulent Vortex Shedding ◽  
Model Equations ◽  
Good Agreement

ABSTRACTThis paper presents a numerical study on turbulent vortex shedding flows past a square cylinder. The 2D unsteady periodic shedding motion was resolved in the calculation and the superimposed turbulent fluctuations were simulated with a second-order Reynolds-stress closure model. The calculations were carried out by solving numerically the fully elliptic ensemble-averaged Navier-Stokes equations coupled with the turbulence model equations together with the two-layer approach in the treatment of the near-wall region. The performance of the computations was evaluated by comparing the numerical results with data from available experiments. Results indicate that the present study gives good agreement in the shedding frequency and mean drag as well as in some phase profiles of the mean velocity.

Modular Turbulence Modeling Applied to an Engine Intake

2013 ◽  
Vol 136 (5) ◽  
Author(s):  
Ugochukwu R. Oriji ◽  
Paul G. Tucker
Keyword(s):  
Surface Roughness ◽  
Turbulence Modeling ◽  
Fluid Dynamic ◽  
Navier Stokes ◽  
Test Cases ◽  
Laminar Separation ◽  
Streamline Curvature ◽  
Favorable Pressure Gradient ◽  
Flow Features ◽  
The One

The one equation Spalart–Allmaras (SA) turbulence model in an extended modular form is presented. It is employed for the prediction of crosswind flow around the lip of a 90 deg sector of an intake with and without surface roughness. The flow features around the lip are complex. There exists a region of high streamline curvature. For this, the Richardson number would suggest complete degeneration to laminar flow. Also, there are regions of high favorable pressure gradient (FPG) sufficient to laminarize a turbulent boundary layer (BL). This is all terminated by a shock and followed by a laminar separation. Under these severe conditions, the SA model is insensitive to capturing the effects of laminarization and the reenergization of eddy viscosity. The latter promotes the momentum transfer and correct reattachment prior to the fan face. Through distinct modules, the SA model has been modified to account for the effect of laminarization and separation induced transition. The modules have been implemented in the Rolls-Royce HYDRA computational fluid dynamic (CFD) solver. They have been validated over a number of experimental test cases involving laminarization and also surface roughness. The validated modules are finally applied in unsteady Reynolds-averaged Navier–Stokes (URANS) mode to flow around an engine intake and comparisons made with measurements. Encouraging agreement is found and hence advances made towards a more reliable intake design framework.

scholarly journals Modeling of Tip Clearance Flows Through an Improved 3-D Pressure Correction Navier-Stokes Solver

1993 ◽  
Author(s):  
N. Lymberopoulos ◽  
K. Giannakoglou ◽  
I. Nikolaou ◽  
K. D. Papailiou ◽  
A. Tourlidakis ◽  
...  
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Three Dimensional ◽  
Tip Clearance ◽  
Navier Stokes ◽  
Special Treatment ◽  
Pressure Correction ◽  
Compressor Rotor ◽  
Numerical Predictions ◽  
Correction Technique

Mechanical constraints dictate the existence of tip clearances in rotating cascades, resulting to a flow leakage through this clearance which considerably influences the efficiency and range of operation of the machine. Three-dimensional Navier-Stokes solvers are often used for the numerical study of compressor and turbine stages with tip-clearance. The quality of numerical predictions depends strongly on how accurately the blade tip region is modelled; in this respect the accurate modelling of tip region was one of the main goals of this work. In the present paper, a 3-D Navier-Stokes solver is suitably adapted so that the flat tip surface of a blade and its sharp edges could be accurately modelled, in order to improve the precision of the calculation in the tip region. The adapted code solves the fully elliptic, steady, Navier-Stokes equations through a space-marching algorithm and a pressure correction technique; the H-type topology is retained, even in cases with thick leading edges where a special treatment is introduced herein. The analysis is applied to two different cases, a linear cascade and a compressor rotor, and comparisons with experimental data are provided.

scholarly journals On the Importance of the Influence of both Velocity and Aspect Ratios on the Occurrence of Bifurcation Phenomena within a Two-Sided Lid-Driven Enclosure

2015 ◽  
Vol 55 ◽  
pp. 160-172
Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Numerical Results ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Left Wall ◽  
Driven Cavity ◽  
Aspect Ratios ◽  
The Stability ◽  
The Right

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.

Numerical Study on Supersonic Combustor with an Allotype Cavity

2014 ◽  
Vol 694 ◽  
pp. 187-192
Author(s):  
Jin Xiang Wu ◽  
Jian Sun ◽  
Xiang Gou ◽  
Lian Sheng Liu
Keyword(s):  
Alternative Model ◽  
Numerical Study ◽  
Reverse Flow ◽  
Three Dimensional ◽  
Flow Characteristics ◽  
Experimental Result ◽  
Navier Stokes ◽  
Cold Flow ◽  
Hypersonic Inlet ◽  
Good Agreement

The three-dimensional coupled explicit Reynolds Averaged Navier–Stokes (RANS) equations and the two equation shear-stress transport k-w (SST k-w) model has been employed to numerically simulate the cold flow field in a special-shaped cavity-based supersonic combustor. In a cross-section shaped rectangular, hypersonic inlet with airflow at Mach 2.0 chamber, shock structures and flow characteristics of a herringbone-shaped boss and a herringbone-shaped cavity models were discussed, respectively. The results indicate: Firstly, according to the similarities of bevel-cutting shock characteristics between the boss case and the cavity case, the boss structure can serve as an ideal alternative model for shear-layer. Secondly, the eddies within cavity are composed of herringbone-spanwise vortexes, columnar vortices in the front and main-spanwise vortexes in the rear, featuring tilting, twisting and stretching. Thirdly, the simulated bottom-flow of cavity is in good agreement with experimental result, while the reverse flow-entrainment resulting from herringbone geometry and pressure gradient. However, the herringbone-shaped cavity has a better performance in fuel-mixing.

Numerical Study of Deterministic Fluxes in Compressor Passages

2018 ◽  
Vol 140 (10) ◽  
Author(s):  
Feng Wang ◽  
Mauro Carnevale ◽  
Luca di Mare
Keyword(s):  
High Speed ◽  
Numerical Study ◽  
Navier Stokes ◽  
Design Stage ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Radial Profiles ◽  
And Performance ◽  
Stage Matching ◽  
Unsteady Effects

Computational fluid dynamics (CFD) has been widely adopted in the compressor design process, but it remains a challenge to predict the flow details, performance, and stage matching for multistage, high-speed machines accurately. The Reynolds Averaged Navier-Stokes (RANS) simulation with mixing plane for bladerow coupling is still the workhorse in the industry and the unsteady bladerow interaction is discarded. This paper examines these discarded unsteady effects via deterministic fluxes using semi-analytical and unsteady RANS (URANS) calculations. The study starts from a planar duct under periodic perturbations. The study shows that under large perturbations, the mixing plane produces dubious values of flow quantities (e.g., whirl angle). The performance of the mixing plane can be considerably improved by including deterministic fluxes into the mixing plane formulation. This demonstrates the effect of deterministic fluxes at the bladerow interface. Furthermore, the front stages of a 19-blade row compressor are investigated and URANS solutions are compared with RANS mixing plane solutions. The magnitudes of divergence of Reynolds stresses (RS) and deterministic stresses (DS) are compared. The effect of deterministic fluxes is demonstrated on whirl angle and radial profiles of total pressure and so on. The enhanced spanwise mixing due to deterministic fluxes is also observed. The effect of deterministic fluxes is confirmed via the nonlinear harmonic (NLH) method which includes the deterministic fluxes in the mean flow, and the study of multistage compressor shows that unsteady effects, which are quantified by deterministic fluxes, are indispensable to have credible predictions of the flow details and performance of compressor even at its design stage.

Three-Dimensional Thermocapillary Convection in Open Cylindrical Annuli

2000 ◽  
Author(s):  
Bok-Cheol Sim ◽  
Abdelfattah Zebib
Keyword(s):  
Reynolds Number ◽  
Aspect Ratio ◽  
Critical Frequency ◽  
Thermocapillary Convection ◽  
Three Dimensional ◽  
Time Dependent ◽  
Infinite Layer ◽  
Aspect Ratios ◽  
Temperature Oscillations ◽  
Good Agreement

Abstract Three-dimensional, time-dependent thermocapillary convection in open cylindrical containers is investigated numerically. Results for aspect ratios (Ar) of 1, 2.5, 8, and 16 and a Prandtl number of 6.84 are obtained to compare the results of numerical simulations with ongoing experiments. Convection is steady and axisymmetric at sufficiently low values of the Reynolds number (Re). Transition to oscillatory states occurs at critical values of Re which depend on Ar. With Ar = 1.0 and 2.5, we observe, respectively, 5 and 9 azimuthal wavetrains travelling clockwise at the free surface near the critical Re. With Ar = 8.0 and 16.0, there are substantially more, but pulsating waves near the critical Re. In the case of Ar = 16.0, which approaches the conditions in an infinite layer, our results are in good agreement with linear theory. While the critical Reynolds number decreases with increasing aspect ratio in the case of azimuthal rotating waves, it increases with increasing aspect ratio in the case of azimuthal pulsating waves. The critical frequency of temperature oscillations is found to decrease linearly with increasing Ar. We have also computed supercritical time-dependent states and find that while the frequency increases with increasing Re near the critical region, the frequency of supercritical convection decreases with Re.

Numerical Performance of Transition Models in Different Turbomachinery Flow Conditions: A Comparative Study

2000 ◽  
Author(s):  
Jiasen Hu ◽  
Torsten H. Fransson
Keyword(s):  
Reynolds Number ◽  
Pressure Gradient ◽  
Numerical Study ◽  
Adverse Pressure Gradient ◽  
Navier Stokes ◽  
Test Cases ◽  
Pressure Gradients ◽  
Flow Conditions ◽  
Transition Models ◽  
High Reynolds

A numerical study has been performed to compare the overall performance of three transition models when used with an industrial Navier-Stokes solver. The three models investigated include two experimental correlations and an integrated eN method. Twelve test cases in realistic turbomachinery flow conditions have been calculated. The study reveals that all the three models can work numerically well with an industrial Navier-Stokes code, but the prediction accuracy of the models depends on flow conditions. In general, all the three models perform comparably well to predict the transition in weak or moderate adverse pressure-gradient regions. The two correlations have the merit if the transition starts in strong favorable pressure-gradient region under high Reynolds number condition. But only the eN method works well to predict the transition controlled by strong adverse pressure gradients. The three models also demonstrate different capabilities to model the effects of turbulence intensity and Reynolds number.

A Numerical Study of Labyrinth Seal Flutter

2010 ◽  
Vol 132 (2) ◽  
Author(s):  
L. di Mare ◽  
M. Imregun ◽  
J. S. Green ◽  
A. I. Sayma
Keyword(s):  
Numerical Study ◽  
Large Diameter ◽  
Labyrinth Seal ◽  
Navier Stokes ◽  
Time Step ◽  
Modal Model ◽  
Flutter Analysis ◽  
Time Histories ◽  
Structural Motion ◽  
Good Agreement

A numerical study of a labyrinth-type turbine seal flutter in a large turbofan engine is described. The flutter analysis was conducted using a coupled fluid-structure interaction code, which was originally developed for turbomachinery blade applications. The flow model is based on an unstructured, implicit Reynolds-averaged Navier–Stokes solver. The solver is coupled to a modal model for the structure obtained from a standard structural finite element code. During the aeroelasticity computations, the aerodynamic grid is moved at each time step to follow the structural motion, which is due to unsteady aerodynamic forces applied onto the structure by the fluid. Such an integrated time-domain approach allows the direct computation of aeroelastic time histories from which the aerodynamic damping, and hence, the flutter stability, can be determined. Two different configurations of a large-diameter aeroengine labyrinth seal were studied. The first configuration is the initial design with four fins, which exhibited flutter instability during testing. The second configuration is a modified design with three fins and a stiffened ring. The steady-state flow was computed for both configurations, and good agreement was reached with available reference data. An aeroelasticity analysis was conducted next for both configurations, and the model was able to predict the observed flutter behavior in both cases. A flutter mechanism is proposed, based on the matching of the structural frequencies to the frequencies of waves traveling in the fluid, in the interfin cavities and in the high- and low-pressure cavities.

scholarly journals Fracture Shape and Orientation Contributions to P-Wave Velocity and Anisotropy of Alpine Fault Mylonites

2021 ◽  
Vol 9 ◽  
Author(s):  
Jirapat Charoensawan ◽  
Ludmila Adam ◽  
Michael Ofman ◽  
Virginia Toy ◽  
Jonathan Simpson ◽  
...  
Keyword(s):  
Aspect Ratio ◽  
Shear Zone ◽  
Numerical Study ◽  
Numerical Models ◽  
Propagation Model ◽  
P Wave ◽  
Micro Ct ◽  
Alpine Fault ◽  
Aspect Ratios ◽  
Ct Data

P-wave anisotropy is significant in the mylonitic Alpine Fault shear zone. Mineral- and texture-induced anisotropy are dominant in these rocks but further complicated by the presence of fractures. Electron back-scattered diffraction and synchrotron X-ray microtomography (micro-CT) data are acquired on exhumed schist, protomylonite, mylonite, and ultramylonite samples to quantify mineral phases, crystal preferred orientations, microfractures, and porosity. The samples are composed of quartz, plagioclase, mica and accessory garnet, and contain 3–5% porosity. Based on the micro-CT data, the representative pore shape has an aspect ratio of 5:2:1. Two numerical models are compared to calculate the velocity of fractured rocks: a 2D wave propagation model, and a differential effective medium model (3D). The results from both models have comparable pore-free fast and slow velocities of 6.5 and 5.5 km/s, respectively. Introducing 5% fractures with 5:2:1 aspect ratio, oriented with the longest axes parallel to foliation decreases these velocities to 6.3 and 5.0 km/s, respectively. Adding both randomly oriented and foliation-parallel fractures hinders the anisotropy increase with fracture volume. The anisotropy becomes independent of porosity when 80% of fractures are randomly oriented. Modeled anisotropy in 2D and 3D are different for similar fracture aspect ratios, being 30 and 15%, respectively. This discrepancy is the result of the underlying assumptions and limitations. Our numerical results explain the effects that fracture orientations and shapes have on previously published field- and laboratory-based studies. Through this numerical study, we show how mica-dominated, pore-free P-wave anisotropy compares to that of fracture volume, shape and orientation for protolith and shear zone rocks of the Alpine Fault.

Sign in / Sign up

Export Citation Format

Share Document