Euler flow in a supersonic mixed‐compression inlet

Numerical simulation of a two‐dimensional mixed compression supersonic inlet is carried out by solving unsteady compressible Euler equations via a stabilized finite element method. The geometry of the inlet is similar to the one used by Anderson and Wong for experimental investigation for Mach 3 flow. The computations are capable of simulating the start‐up problems associated with the inlet. The critical back pressure for the successful operation of the inlet is computed. The effect of inlet back pressure on the total pressure recovery and the flow distortion level is analysed. Contrary to the popular belief, it is found that in addition to the throat to inlet capture area ratio, the ramp geometry close to the throat plays an important role in the start‐up dynamics. It is demonstrated via simulations that, everything else being same, the geometries of ramp upstream of the throat that are associated with a curvature higher than a certain threshold, result in unstarting the intake. Copyright © 2005 John Wiley & Sons, Ltd.     

Progress of the unified wave-particle methods for non-equilibrium flows from continuum to rarefied regimes

We utilize the nonlinear acoustic solver (NLAS) and Ffowcs-Williams/Hawkings (FW-H) equation to investigate the noise generation and radiation due to shock (wave) and boundary layer interaction (SBLI) in the inlet duct. A classical benchmark for SBLI is chosen to validate the flow features and numerical results show good agreement with experimental results. In the simulation of the noise generated by SBLI, the inlet buzz phenomenon is successfully observed. The oscillation of the normal shock is a kind of little buzz and the oscillation of inner shocks is a kind of big buzz with a frequency around 100 Hz. In the far-field, frequency spectrums show a dominant frequency close to the frequency of inner shocks oscillation. This indicates that the oscillation of inner shocks determines the magnitude of the overall sound pressure level (OASPL) of the far-field noise.

     

A simplified marker and cell method for unsteady flows on non-staggered grids

The SMAC (simplified marker and cell) time-advancing method for solving the unsteady incompressible Navier-Stokes equations on non-staggered grids is developed in generalized co-ordinate systems. The primitive variable formulation uses Cartesian velocities and pressure, all defined at the centre of the control volume, as the dependent variables. A special elliptic flux correction at the faces of the finite volume is utilized in discretizing the continuity equation to suppress pressure oscillations. The test flows considered are a polar cavity flow starting from rest and the flow around a circular cylinder. The numerical results are compared with experimental results and results obtained by the well-known SIMPLEC and PISO methods. The comparisons show that the elliptic flux correction technique works well in suppressing pressure oscillations and that the SMAC method is more efficient than the SIMPLEC and PISO methods for both steady and unsteady flows.     

Hall effects on unsteady flow of a viscous fluid due to non-coaxial rotation of a porous disk and a fluid at infinity

An exact solution of the unsteady hydromagnetic flow due to non-coaxial rotations of a porous disk and a fluid at infinity is obtained on taking Hall currents into account. An analytical solution of the problem is obtained for small and large times after the start by the Laplace transform method. It is found that for small values of time there is no inertial oscillations while for large time the steady state is reached through inertial oscillations. The frequency of these oscillations first increases, reaches a maximum and then decreases with increase in Hall parameter.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

Surface pressure and heat transfer measurements in a turbine cascade with unsteady oncoming wakes

Unsteady surface pressure and heat transfer have been measured on a blade of a linear turbine cascade exposed to unsteady oncoming wakes generated by moving cylinders on a squirrel cage device. The Reynolds number and the Strouhal number corresponded to the values in a real turbomachine. The periodic components of pressure and heat transfer showed clear response to the unsteady wakes. However, the distribution along the blade surface of both pressure and heat transfer coefficient changed very little from phase to phase. The heat transfer results have shown that the boundary layer on the pressure side remained laminar for all cases, but that the boundary layers on the suction side became transitional under the wake disturbance. With increasing wake-passing frequency, the start of the transition moved forward. Increasing the wake-passing frequency resulted in a significant increase in heat transfer along the whole blade surface including the portions where the boundary layers were nominally laminar.     

Acoustic resonance phenomena in air bleed channels in aviation engines

The existence of axial–radial acoustic resonance oscillations of the basic air flow in bleed channels of aviation engines is demonstrated theoretically and experimentally. Numerical and analytical methods are used to determine the frequency of acoustic resonance oscillations for the lowest modes of open and closed bleed channels of the PS-90A engine. Experimental investigations reveal new acoustic resonance phenomena arising in the air flow in bleed channel cavities in the core duct of this engine owing to instability of the basic air flow. The results of numerical, analytical, and experimental studies of the resonance frequencies reached in the flow in bleed channel cavities in the core duct of the PS-90A engine are found to be in reasonable agreement. As a result, various types of resonance oscillations in bleed channels can be accurately described.     

Unsteady RANS simulation of oscillating mould flows

Mould flow oscillations are of major importance for the performance of the continuous casting process. They are suspected to promote entrainment of slag and other unwanted secondary phases into the melt pool. These oscillating turbulent flows are investigated by means of numerical simulations. The numerical model is based on the equation of continuity and the unsteady Reynolds averaged Navier–Stokes equations. The system of flow equations is closed by a Reynolds stress turbulence model in combination with non‐equilibrium wall functions. The unsteady simulation resolves low‐frequency oscillations of the flow field. These frequencies and numerically resolved mean values are in agreement with results of corresponding model experiments. The proposed model should be advantageous in order to investigate the mechanisms of the oscillations and the process of slag entrainment in more detail. Copyright © 2006 John Wiley & Sons, Ltd.     

Limit cycle oscillation behavior of transonic control surface buzz considering free-play nonlinearity

The limit cycle oscillation (LCO) behaviors of control surface buzz in transonic flow are studied. Euler equations are employed to obtain the unsteady aerodynamic forces for Type B and Type C buzz analyses, and an all-movable control surface model, a wing/control surface model and a three-dimensional wing with a full-span control surface are adopted in the study. Aerodynamic and structural describing functions are used to deal with aerodynamic and structural nonlinearities, respectively. Then the buzz speed and buzz frequency are obtained by V-g method. The LCO behavior of the transonic control surface buzz system with linear structure exhibits subcritical or supercritical bifurcation at different Mach numbers. For nonlinear structural model with a free-play nonlinearity in the control surface deflection stiffness, the double LCO phenomenon is observed in certain range of flutter speed. The free-play nonlinearity changes the stability of LCOs at small amplitudes and turns the unstable LCO into a stable one. The LCO behavior is dominated by the aerodynamic nonlinearity for the case with large control surface oscillation amplitude but by the structural nonlinearity for the case with small amplitude. Good agreements between LCO behaviors obtained by the present method and available experimental data show that our study may help to explain the experimental observation in wind tunnel tests and to understand the physical mechanism of transonic control surface buzz.     

Oscillations of a viscous fluid in a thin layer between two coaxial disks rotating together

A study is made of the problem of the motion of an incompressible viscous fluid in the space between two coaxial disks rotating together with constant angular velocity under the assumption that the pressure changes in time in accordance with a harmonic law. The problem is solved using the equations of unsteady motion of an incompressible viscous fluid in a thin layer. It is shown that the velocity field in this case is a superposition on a steady field of damped oscillations with cyclic frequency equal to twice the angular velocity of the disks and forced oscillations with cyclic frequency equal to the cyclic frequency of the oscillations of the pressure field. It is shown that the amplitude of the forced oscillations of the velocity field depends strongly on the ratio of the cyclic frequency of the oscillations of the pressure field to the angular velocity of the disks. It is shown that there is a certain value of the ratio at which the amplitude of the forced oscillations has a maximal value (resonance). It is shown that even for very small amplitudes of the pressure oscillations the amplitude of the oscillations of the relative velocity at resonance may reach values comparable with the mean velocity of the main flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 166–169, January–February, 1984.     

On the stabilization of unphysical pressure oscillations in MPS method simulations

In this paper, pressure stability through the suppression of high‐frequency pressure oscillations in the moving particle semi‐implicit (MPS) method is presented. To obtain a stable pressure field, we improve the free‐surface particle search algorithm. Pressure stability follows from the suppression of high‐frequency pressure oscillations due to a correction in the Laplacian operator of the Poisson pressure equation and from the correction of the pressure gradient operator. The three proposed modifications are applied gradually and compared with the MPS method to show the improvements in the hydrostatic pressure and dam‐breaking problems. To validate the suppression of the high‐frequency numerical pressure oscillations, modified MPS methods with and without a removable wall are compared with published dam‐breaking experiment pressure measurements. Copyright © 2016 John Wiley & Sons, Ltd.     

Computations of two passing‐by high‐speed trains by a relaxation overset‐grid algorithm

This paper presents a relaxation algorithm, which is based on the overset grid technology, an unsteady three‐dimensional Navier–Stokes flow solver, and an inner‐ and outer‐relaxation method, for simulation of the unsteady flows of moving high‐speed trains. The flow solutions on the overlapped grids can be accurately updated by introducing a grid tracking technique and the inner‐ and outer‐relaxation method. To evaluate the capability and solution accuracy of the present algorithm, the computational static pressure distribution of a single stationary TGV high‐speed train inside a long tunnel is investigated numerically, and is compared with the experimental data from low‐speed wind tunnel test. Further, the unsteady flows of two TGV high‐speed trains passing by each other inside a long tunnel and at the tunnel entrance are simulated. A series of time histories of pressure distributions and aerodynamic loads acting on the train and tunnel surfaces are depicted for detailed discussions. Copyright © 2004 John Wiley & Sons, Ltd.     

The influence of boundary layers on supersonic inlet flow unstart induced by mass injection

A transverse jet is injected into a supersonic model inlet flow to induce unstart. Planar laser Rayleigh scattering from condensed CO₂ particles is used to visualize flow dynamics during the unstart process, while in some cases, wall pressure traces are simultaneously recorded. Studies conducted over a range of inlet configurations reveal that the presence of turbulent wall boundary layers strongly affect the unstart dynamics. It is found that relatively thick turbulent boundary layers in asymmetric wall boundary layer conditions prompt the formation of unstart shocks; in symmetric boundary conditions lead to the propagation of pseudo-shocks; and in both cases facilitate fast inlet unstart, when compared with thin, laminar boundary layers. Incident shockwaves and associated reflections are found to affect the speed of pressure disturbances. These disturbances, which induce boundary layer separation, are found to precede the formation of unstart shocks. The results confirm the importance of and need to better understand shock-boundary layer interactions in inlet unstart dynamics.     

Turbulence model and numerical scheme assessment for buffet computations

The prediction of shock‐induced oscillations over transonic rigid airfoils is important for a better understanding of the buffeting phenomenon. The unsteady resolution of the Navier–Stokes equations is performed with various transport‐equation turbulence models in which corrections are added for non‐equilibrium flows. The lack of numerical efficiency due to the CFL stability condition is circumvented by the use of a wall law approach and a dual time stepping method. Moreover, various numerical schemes are used to try and be independent of the numerical discretization. Comparisons are made with the experimental results obtained for the supercritical RA16SC1 airfoil. They show the interest in using the SST correction or realizability conditions to get correct predictions of the frequency, amplitude and pressure fluctuations over the airfoil. Copyright © 2004 John Wiley & Sons, Ltd.     

Starting flow in a rotating parallel plate channel

The starting flow due to a suddenly applied pressure gradient in a parallel plate channel which is rotating as a system is studied. Exact analytic series solutions to the unsteady Navier-Stokes equations are found by both the Laplace transform method and the separation of parameters method, the latter is shown to be superior. Rotation not only induces a secondary transverse flow but also alters the character of the transient flow rate and velocity profiles. Back flow and inertial oscillations occur, especially at higher rotation rates.     

Adaptive strategy of transonic flows over vibrating blades with interblade phase angles

An error indicator and a locally implicit scheme with anisotropic dissipation model on dynamic quadri‐ lateral–triangular mesh are developed to study transonic flows over vibrating blades with interblade phase angles. In the Cartesian co‐ordinate system, the unsteady Euler equations with moving domain effects are solved. The error indicator, in which unified magnitudes of dynamic grid speed, substantial derivative of pressure, and substantial derivative of vorticity magnitude are incorporated to capture the unsteady wave behaviours and vortex‐shedding phenomena due to unsteadiness. To assess the accuracy of the locally implicit scheme with anisotropic dissipation model on quadrilateral–triangular mesh, two flow calculations are performed. Based on the comparison with the related numerical and experimental data, the accuracy of the present approach is confirmed. According to the high‐resolutional result on the adaptive mesh, the unsteady pressure wave, shock and vortex‐shedding behaviours are clearly demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Buffeting in transonic flow prediction using time‐dependent turbulence model

In transonic flow conditions, the shock wave/turbulent boundary layer interaction and flow separations on wing upper surface induce flow instabilities, ‘buffet’, and then the buffeting (structure vibrations). This phenomenon can greatly influence the aerodynamic performance. These flow excitations are self‐sustained and lead to a surface effort due to pressure fluctuations. They can produce enough energy to excite the structure. The objective of the present work is to predict this unsteady phenomenon correctly by using unsteady Navier–Stokes‐averaged equations with a time‐dependent turbulence model based on the suitable (k–ε) turbulent eddy viscosity model. The model used is based on the turbulent viscosity concept where the turbulent viscosity coefficient (Cμ) is related to local deformation and rotation rates. To validate this model, flow over a flat plate at Mach number of 0.6 is first computed, then the flow around a NACA0012 airfoil. The comparison with the analytical and experimental results shows a good agreement. The ONERA OAT15A transonic airfoil was chosen to describe buffeting phenomena. Numerical simulations are done by using a Navier–Stokes SUPG (streamline upwind Petrov–Galerkin) finite‐element solver. Computational results show the ability of the present model to predict physical phenomena of the flow oscillations. The unsteady shock wave/boundary layer interaction is described. Copyright © 2005 John Wiley & Sons, Ltd.     

Oscillations of a wedge in a free stream with particular reference to parachutes

The phenomenon of coupled breathing and longitudinal oscillations of a wedge-mass system in a free stream is examined. As a first step the unsteady pressure distribution on the surface of the oscillating wedge is calculated. For dynamic equilibrium of the wedge-mass system, the moment about the apex of the wedge must be zero. This condition establishes the amplitude and phase relation between breathing and longitudinal oscillations. As a final step the equation of motion of the store is used to calculate the frequency of the breathing oscillations. This frequency is shown to be dependent on four parameters. These parameters include the Froude number, the rigging line length to wedge breadth ratio and the rigging line stiffness and damping. Current results are compared with Hume and Stevens [1] experimental results.     

Measurements of pulsating and oscillating flows of non-Newtonian fluids through concentric and eccentric cylinders

The unsteady flow of non-Newtonian fluids through concentric and eccentric cylinders was investigated experimentally. Two experiments were carried out; one was pulsating flow and the other was flow under a constant pressure gradient with the inner cylinder oscillating longitudinally. The flow enhancement was examined and its dependence on the frequency of the oscillations and the eccentricity of the apparatus was determined.     

带喷流激波针流动特性实验研究

采用动态测力、动态测压和纹影等风洞实验技术,对加装了带喷流激波针的钝头体的绕流特性、稳定和非稳模态的形成条件和机理进行了研究.结果表明:带喷流激波针流场存在稳态和非稳态两种模态,超声速喷流的压比大于临界压比时流动处于稳定模态,反之则为非稳模态;增大激波针长度可减小钝头体阻力,但达到一定长度后,进一步减阻的效果不再显著;增大喷流压比能够有效减弱再附激波强度,有利于缓解单独激波针的肩部热斑问题;非稳模态下波系自激振荡对再附激波在钝头体表面所围的区域影响剧烈,振荡是周期性的,且存在确定的主导频率,主导频率随喷流压力比增大而减小;自激振荡的产生是由于喷流出口周围的反压在喷流压比小于临界压比时无法获得持续的平衡而导致.