Shock wave structures ahead of nonuniform fan cascades

We study shock wave structures (SWS), consisting of shock waves and expansion waves between them, that occur in supersonic flow past nonuniform fan cascades when the velocity component normal to their front (“axial” component) is subsonic. The cascade nonuniformity is due to the scatter in the setting angles of identical blades, either sharp or blunt. A result of the uniformity is the generation of combined noise, whose frequencies are much smaller than the fundamental frequency of the uniform cascade, and slower nonlinear SWS attenuation. The accurate and fast “simple wave method” and “nonlinear acoustics approximation”, together with numerical algorithms for integrating Euler equations on overlapping grids (in calculating flow past blunt edges) and on SWS-adapted grids, are applied to determine the “guiding” action of nonuniform cascades and to describe the SWS evolution. The application of the Fourier analysis gives the sound field spectrum. The use of blades with rectilinear initial regions of the “backs” for reducing supersonic fan blade noise is efficient only at small (less than 0.25°) scatter in the setting angles. The shock wave structures attenuate more rapidly ahead of nonuniform cascades composed of blunt blades than ahead of those with sharp blades. For uniform cascades the blade bluntness effect is not large.     

Level set method for numerical simulation of a cavitation bubble,its growth,collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier–Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving interface is captured by the level set function, and the interface velocity is resolved by “one-side” velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the “shock wave”-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating. The project supported by the National Natural Science Foundation of China (10272032 and 10672043). The English text was polished by Keren Wang.     

Turbulent boundary layer velocity profiles over a smooth surface

Mean velocity measurements made for two-dimensional incompressible turbulent boundary layers, for 3 pressure gradients on a smooth wall have been used to evaluate the pressure gradient parameter, the skin-friction and the roughness function by the method of “curve-fit” to the entire profile. Local skin-friction coefficients so obtained do not agree with the experimental results but the integral parameters are well represented.     

Vortex structure simulation for supersonic mixing layers using nonlinear PSE method

The method of nonlinear parabolized stability equations (PSE) is applied in the simulation of vortex structures in compressible mixing layer. The spatially-evolving unstable waves, which dominate the vortex structure, are investigated through spatial marching method. The instantaneous flow field is obtained by adding the harmonic waves to basic flow. The results show that T-S waves do not keep growing exponentially as the linear evolution, the energy transfer to high order harmonic modes, and that finally all harmonic modes get saturated due to nonlinear interaction. The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones, makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up, vortex pairing, and Λ vortex.     

Laminar–turbulent transition in oscillating boundary layer: experimental and numerical analysis using continuous wavelet transform

The laminar–turbulent transition of a forced oscillating boundary layer with a varying pressure gradient is experimentally and numerically investigated for two Strouhal numbers. Time-dependent characterization of the natural instability modes is carried out using continuous wavelet analysis of velocity signals. The periodic evolution of the total growth rates of the most unstable disturbances are measured and compared to the results of the linear stability theory. The “Tollmien–Schlichting” (TS) and “convective” transition modes are identified. It is shown that they correspond to the extrema of opposite signs of the skewness factor of unstable wavelet transform of the most unstable frequencies.     

Artificial Boundary Conditions of Pressure Type for Viscous Flows in a System of Pipes

In a three-dimensional domain Ω with J cylindrical outlets to infinity the problem is treated how solutions to the stationary Stokes and Navier–Stokes system with pressure conditions at infinity can be approximated by solutions on bounded subdomains. The optimal artificial boundary conditions turn out to have singular coefficients. Existence, uniqueness and asymptotically precise estimates for the truncation error are proved for the linear problem and for the nonlinear problem with small data. The results include also estimates for the so called “do-nothing” condition.     

Spatial direct numerical simulation of high-speed boundary-layer flows Part II: Transition on a cone in Mach 8 flow

The laminar breakdown of the boundary-layer flow of an axisymmetric sharp cone in a Mach 8 flow is simulated by a synergistic approach that combines the parabolized stability equation (PSE) method and spatial direct numerical simulation (DNS). The transitional state is triggered by a symmetric pair of oblique second-mode disturbances whose nonlinear interactions generate strong streamwise vorticity, which leads in turn to severe spanwise variations in the flow and eventual laminar breakdown. The PSE method is used to compute the weakly and moderately nonlinear initial stages of the transition process and, thereby, to derive a harmonically rich inflow condition for the DNS. The strongly nonlinear and laminar-breakdown stages of transition are then computed by well-resolved DNS, with a highly accurate algorithm that exploits spectral collocation and high-order compact-difference methods. Evolution of the flow is presented in terms of modal energies, mean quantities (e.g., skin friction), Reynolds stresses, turbulent kinetic energy, and flow visualization. The numerical test case is an approximate computational analog of one of the few stability experiments performed for hypersonic boundary-layer flows. Comparisons and contrasts are drawn between the experimental and the computational results. Rope-like waves similar to those observed in schlieren images of high-speed transitional flows are also observed in the numerical experiment and are shown to be visual manifestations of second-mode instability waves.This research was supported under NASA Contracts NAS1-19831 and NAS1-20059 for the first and second authors, respectively.     

Studies of quality of geometrically nonlinear elasticity theory for small strains and arbitrary displacements

Earlier it was shown in [1, 2] that the equations of classical nonlinear elasticity constructed for the case of small strains and arbitrary displacements are ill posed, because their use in specific problems may result in the appearance of “spurious” bifurcation points. A detailed analysis of these equations and the construction, in their stead, of consistent equations of geometrically nonlinear theory of elasticity can be found in [3]. Certain steps in this direction were also made in [4, 5]. In [3], it was also stated that the methods and applied program packages (APPs) based on the use of the classical relations of nonlinear elasticity require some revision and correction. In the present paper, this conclusion is justified and confirmed by numerical finite-element solutions of several three-dimensional geometrically nonlinear deformation problems and linearized problems on the stability of equilibrium of a rectilinear beam. These solutions were obtained by using two APPs developed by the authors and the well-known APP “ANSYS.” It is shown that the classical equations of the geometrically nonlinear theory of elasticity, which underly the first of the developed APP and the well-known APP “ANSYS,” often lead to overestimated buckling loads for structural members as compared with the consistent equations proposed in [1–3].     

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations ( PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub- characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic., respectively . The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories , the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time , the methods of removing the remained ellipticity are further obtained from the nonlinear PSE .     

Nonlinear evolution analysis of T-S disturbance wave at finite amplitude in nonparallel boundary layers

ＩｎｔｒｏｄｕｃｔｉｏｎＮｏｎｌｉｎｅａｒｓｔａｂｉｌｉｔｙｐｒｏｂｌｅｍｉｎｎｏｎｐａｒａｌｌｅｌｂｏｕｎｄａｒｙｌａｙｅｒｓｎｅｅｄｓｍａｉｎｌｙｃｏｎｓｉｄｅｒｔｈｅｅｖｏｌｕｔｉｏｎｏｆｕｎｓｔａｂｉｌｉｔｙＴ＿Ｓｗａｖｅａｎｄｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｗａｖｅｓｏｆｄｉｆｆｅｒｅｎｔｆｒｅｑｕｅｎｃｉｅｓ.Ｉｔｉｓｖｅｒｙｉｍｐｏｒｔａｎｔｔｏｓｉｍｕｌａｔｅａｃｃｕｒａｔｅｌｙｓｐａｔｉａｌｅｖｏｌｕｔｉｏｎｒｏｌｅｏｆｄｉｓｔｕｒｂａｎｃｅｉｎｂｏｕｎｄａｒｙｌａｙｅｒｓ,…     

Applications of parabolized stability equation for predicting transition position in boundary layers

The phenomenon of laminar-turbulent transition exists universally in nature and various engineering practice.The prediction of transition position is one of crucial theories and practical problems in fluid mechanics due to the different characteristics of laminar flow and turbulent flow.Two types of disturbances are imposed at the entrance,i.e.,identical amplitude and wavepacket disturbances,along the spanwise direction in the incompressible boundary layers.The disturbances of identical amplitude are consisted of one two-dimensional(2D) wave and two three-dimensional(3D) waves.The parabolized stability equation(PSE) is used to research the evolution of disturbances and to predict the transition position.The results are compared with those obtained by the numerical simulation.The results show that the PSE method can investigate the evolution of disturbances and predict the transition position.At the same time,the calculation speed is much faster than that of the numerical simulation.     

Mathematical programming analysis of bilinear rheological lubrication mechanics

Bilinear rheological lubrication mechanics provides an important basis for the designs of recently developed electrorheological (ER) “smart”journal bearings and those lubricated by mixed fluid-solid lubricants. But there is not yet a reliable and efficient numerical method for such a problem of non-Newtonian fluid mechanics. In the present paper, a finite element method (FEM) together with mat hematical programming solution is successfully used to solve such a problem. A reliable and generalized numerical method for the designs of electrorheological “smart” journal bearings and the bearings lubricated by mixed fluid-solid lubricant is presented.     

Numerical modeling of the dynamics of a viscous-fluid droplet

A numerical model developed on the basis of the level set method is proposed. This makes it possible to describe both the nonlinear oscillations of a single viscous-fluid droplet and the fragmentation and coalescence processes. The Navier-Stokes equations written in “velocity-pressure” variables on a rectangular uniform grid in cylindrical coordinates are solved using the method of splitting into physical processes. Non-oscillating solutions for two-phase media with a characteristic density ratio of less than 10⁻³ and Re > 1000 are obtained. The possibilities of the approach proposed are demonstrated with the reference to the problem of a droplet falling from a capillary (detachment from the capillary, formation of a “Plato ball”, droplet motion, collision with a plane wall, droplet oscillations on the wall, and droplet spreading). A comparison of the numerical results with the known calculation models and experimental data shows satisfactory agreement with respect to both the phases and the shape of the droplet.     

Lower bound limit analysis of three-dimensional elastoplastic structures by boundary element method

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

Imaging spectroscopy of the nonequilibrium shock front radiation in air

The nonequilibrium radiation of shock fronts in air is experimentally investigated by means of the imaging spectroscopy technique. Shock velocity ranges from 9.7 to 11.6 km/s and initial pressure from 13.3 to 41.6 Pa. The spectral diagnostic system consists of an imaging spectrograph, a streak camera, a gated image-intensified CCD camera and a personal computer for data acquisition/processing. This spectral diagnostic system is capable of simultaneous wavelength-, intensity- and time-resolved spectroscopic measurements in the nanosecond order. The image processing of the streak images includes a combined smoothing/deconvolution process in the time direction to diminish experimental noise effects and the temporal broadening due to the streak camera entrance slit. Wavelength range is chosen to investigate the first negative band of . “Large” and “slim” streak image types are observed. In the “large” streak images greater contribution from (1-)(1,0) behind the radiation peak is observed. Experimental data are compared with a streak image numerically simulated. The numerical simulation agrees better with the “slim” streak image. Received 7 July 1995 / Accepted 10 January 1996     

Verification of parabolized stability equations for its application to compressible boundary layers

Parabolized stability equations (PSE) were used to study the evolution of disturbances in compressible boundary layers.The results were compared with those ob- tained by direct numerical simulations (DNS),to check if the results from PSE method were reliable or not.The results of comparison showed that no matter for subsonic or supersonic boundary layers,results from both the PSE and DNS method agreed with each other reasonably well,and the agreement between temperatures was better than those between velocities.In addition,linear PSE was used to calculate the neutral curve for small amplitude disturbances in a supersonic boundary layer.Compared with those obtained by linear stability theory (LST),the situation was similar to those for incom- pressible boundary layer.     

The “rebel travelling” of general nonlinear evolutional equation and discussion on related problems

This paper is a part of series works for discussing the “auto-destruction effects” of general nonlinear evolutional equations. The blown-up of Navier-Stockes equation is discussed in references [1,2]. Some expansion is made in this paper, and the blown-up of order-1 or 2 models and the “rebel travelling” of complex model of poly-order are discussed. The results indicate that “semi-rupture” appears for some models on specific condition: the blown-up appears during the whole evolution. For fluid, however, the weakly-nonlinear model is of more artificiality and there is much room for arguing about the smoothing scheme of the numerical integral on the basis of continuous thinking and so on.     

Stability of multilayered shells of revolution with allowance for geometric nonlinearity of the subcritical state

The finite-difference method and the Trefftz-Reissner variational principle are used to obtain a system of equations in mixed from to describe the stability and geometric nonlinearity of composite shells of revolution. Methods are developed and an algorithm is proposed to calculate the components of the geometrically nonlinear subcritical stress-strain state and to use those components to determine the “upper” critical values for shells with zero Gaussian curvature loaded by uniform external pressure, an axisymmetric load, or a combination of these loads. The stability of cylindrical, conical, and compound shells under uniform pressure is examined for different support conditions. Linear and nonlinear methods of determining the subcritical stress-strain state are compared and their effect on the critical loads is estimated. Ukrainian Transportation Institute and the Ukrainian Academy of Water Management, Kiev, Ukraine. Translated from Prikladnaya Mekhanika, Vol. 35, No. 6, pp. 60–66, June, 1999.     

Simulations of Spatially Developing Two-Dimensional Shear Layers and Jets

A computational study of spatially evolving two-dimensional free shear flows has been performed using direct numerical simulation of the Navier–Stokes equations in order to investigate the ability of these two-dimensional simulations to predict the overall flow-field quantities of the corresponding three-dimensional “real” turbulent flows. The effects of inflow forcing on these two-dimensional flows has also been studied. Simulations were performed of shear layers, as well as weak (large co-flow and relatively weak shear) and strong (small co-flow and relatively strong shear) jets. Several combinations of discrete forcing with and without a broadband background spectrum were used. Although spatially evolving direct simulations of shear layers have been performed in the past, no such simulations of the plane jet have been performed to the best of our knowledge. It was found that, in the two-dimensional shear layers, external forcing led to a strong increase in the initial growth of the shear-layer thickness, followed by a region of decreased growth as in physical experiments. The final downstream growth rate was essentially unaffected by forcing. The mean velocity profile and the naturally evolving growth rate of the shear layer in the case of broadband forcing compare well with experimental data. However, the total and transverse fluctuation intensities are larger in the two-dimensional simulations with respect to experimental data. In the weak-jet simulations it was found that symmetric forcing completely overwhelms the natural tendency to transition to the asymmetric jet column mode downstream. It was observed that two-dimensional simulations of “strong” jets with a low speed co-flow led to a fundamentally different flow with large differences even in mean velocity profiles with respect to experimental data for planar jets. This was a result of the dominance of the two-dimensional mechanism of vortex dipole ejection in the flow due to the lack of spanwise instabilities. Experimental studies of planar jets do not show vortex dipole formation and ejection. A three-dimensional “strong”-jet simulation showed the rapid evolution of three-dimensionality effectively preventing this two-dimensional mechanism, as expected from experimental results. Received: 25 November 1996 and accepted 17 April 1997     

Nonparallel stability of the flow in an axially rotating pipe

The linear stability of the developing flow in an axially rotating pipe is analyzed using parabolized stability equations (PSE). The results are compared with those obtained from a near-parallel stability approximation that only takes into account the axial variation of the basic flow. Though the PSE results obviously coincide with the near-parallel ones far downstream, when the flow has reached a Hagen-Poiseuille axial velocity profile with superimposed solid-body rotation, they differ significantly in the developing region. Therefore, the onset of instability strongly depends on the axial evolution of the perturbations. The PSE results are also compared with experimental data from Imao et al. [Exp. Fluids 12 (1992) 277], showing a good agreement in the frequencies and wavelengths of the unstable disturbances, that take the form of spiral waves. Finally, a simple method for detecting one of the conditions to characterize the onset of absolute instability using PSE is given.