Numerical Simulation of Two-Dimensional Shock/Boundary-Layer Interaction between a Rocket and Booster

A two-dimensional Reynolds-averaged Navier-Stokes solver is applied to analyze the aerodynamic behavior of the Shock/Boundary-Layer interaction of rocket with a boosted The K-ε turbulence model and a finite volume method in a unstructured body-fitted curvilinear coordinates have been used. The results indicate that the separation and the reattachment occur in the Boundary-Layer of the main rocket because of the shock interaction. The shape of the booster nose effects the flow field obviously. In the case of the hemisphere booster nose the pressure has complicate distributions and the separation is very clear. The distance between the booster and main rocket has the evident effect on the flow field. If the distance is smaller the pressure coefficient is bigger the separation zone even the separation bubble occurs.

Boundary-layer receptivity under interaction of free-stream turbulence and localized wall roughness

The boundary-layer receptivity under the interaction of free-stream turbu- lence （FST） and localized wall roughness is studied by the direct numerical simulation （DNS） and the fast Fourier transform. The results show that the Tollmien-Schlichting （T-S） wave packets superposed by a group of stability, neutral, and instability T-S waves are generated in the boundary layer. The propagation speeds of the T-S wave packets are calculated. The relation among the boundary-layer receptivity response, the amplitude of the FST, the roughness height, and the roughness width is determined. The results agree well with Dietz＇s experiments. The effect of the roughness geometries on the receptivity is also studied.

Boundary-Layer Wind Structure in a Landfalling Tropical Cyclone

In this study, a slab boundary layer model with a constant depth is used to analyze the boundary-layer wind structure in a landfalling tropical cyclone. Asymmetry is found in both the tangential and radial components of horizontal wind in the tropical cyclone boundary layer at landfall. For a steady tropical cyclone on a straight coastline at landfall, the magnitude of the radial component is greater in the offshoreflow side and the tangential component is greater over the sea, slightly offshore, therefore the greater total wind speed occurs in the offshore-flow side over the sea. The budget analysis suggests that： （1） a greater surface friction over land produces a greater inflow and the nonlinear effect advects the maximum inflow downstream, and （2） a smaller surface friction over the sea makes the decrease of the tangential wind component less than that over land. Moreover, the boundary layer wind structures in a tropical cyclone are related to the locations of the tropical cyclone relative to the coastline due to the different surface frictions. During tropical cyclone landfall, the impact of rough terrain on the cyclone increases, so the magnitude of the radial component of wind speed increases in the offshore-flow side and the tangential component outside the radius of maximum wind speed decreases gradually.

Series Solution of Non-Similarity Boundary-Layer Flow in Porous Medium

This paper aims to present complete series solution of non-similarity boundary-layer flow of an incompressible viscous fluid over a porous wedge. The corresponding nonlinear partial differential equations are solved analytically by means of the homotopy analysis method (HAM). An auxiliary parameter is introduced to ensure the convergence of solution series. As a result, series solutions valid for all physical parameters in the whole domain are given. Then, the effects of physical parameters γ and Prandtl number Pr on the local Nusselt number and momentum thickness are investigated. To the best of our knowledge, it is the first time that the series solutions of this kind of non-similarity boundary-layer flows are reported.

Controlling secondary flow in high-lift low-pressure turbine using boundary-layer slot suction

The design of high-lift Low-Pressure Turbines(LPTs)causes the separation of the boundary layer on the suction side of the blade and leads to a strong secondary flow.This present study aims to minimize secondary losses through endwall slot suction and incoming wakes in a front-loaded high-lift LPT cascade with Zweifel of 1.58 under low Reynolds number of 25000.Two slotted schemes for the boundary layer of the endwall were designed(Plan A and Plan B),and the effects of suction mass flow on secondary flow were studied.The underlying physics of the endwall boundary layer of the suction and secondary flow under unsteady wakes was discussed.The results show that slot suction at the endwall boundary layer can significantly suppress the secondary flow by removing low-momentum fluids.Plans A and B significantly reduced the secondary kinetic energy by 44.2%and 36.9%,respectively,compared with the baseline cascade at the suction mass flow ratios of 1%.With an increase in the mass flow ratio of suction,the secondary flow was gradually reduced in both Plans A and B.It is more beneficial to control the secondary flow to destroy the intersection of the pressure side and suction side of the horseshoe vortex before it develops into a passage vortex.Under unsteady wakes,the combined effects of incoming wakes and endwall boundary layer suction can further suppress the secondary flow at the suction mass flow ratios of 2%for Plan A,because the positive and negative vorticity inside upstream wakes accelerated the mixing of the main flow and secondary flow and thus increased the energy of secondary vortices.

Experimental study of bluntness effects on hypersonic boundary-layer transition over a slender cone using surface mounted pressure sensors

In this work,we studied the bluntness effect on the hypersonic boundary-layer transition over a slender cone at Mach 6 with interchangeable tips in a noisy Ludwieg tube tunnel before the so-called“transition reversal”phenomenon occurs.The evolution of instability waves is characterized using surface flush-mounted pressure sensors deployed along the streamwise direction within unit Reynolds number from 4E+6/m≤Reunit≤10E+6/m,and the bluntness of the cone nose ranges from 0.1 mm to 5 mm.Power spectral density(PSD)of pressure fluctuation indicates that small nose bluntness(ReR≤2000)has little influence on the evolution of instability waves along the hypersonic boundary-layer,whereas with a moderate nose size(2000≤ReR≤5000),the hypersonic boundary layer transition is delayed monotonically as the nose radius increases before the boundary-layer turns into fully laminar without instability waves.The delaying effect can be attributed to the increased entropy-layer swallowing distance with a large tip radius.Instability wave characterization reveals that the second mode instability wave plays a dominant role before the transition reversal happens.The quadratic phase locking of second mode instabilities can be identified by bispectral analysis,and it attenuates as the nose tip radius increases.

Boundary-layer transition prediction using a simplified correlation-based model

This paper describes a simplified transition model based on the recently developed correlation-based γ - Reot transition model. The transport equation of transition momentum thick- ness Reynolds number is eliminated for simplicity, and new transition length function and critical Reynolds number correlation are proposed. The new model is implemented into an in-house com- putational fluid dynamics （CFD） code and validated for low and high-speed flow cases, including the zero pressure flat plate, airfoils, hypersonic flat plate and double wedge. Comparisons between the simulation results and experimental data show that the boundary-layer transition phenomena can be reasonably illustrated by the new model, which gives rise to significant improvements over the fully laminar and fully turbulent results. Moreover, the new model has comparable features of accuracy and applicability when compared with the original 3＇ - Reot model. In the meantime, the newly proposed model takes only one transport equation of intermittency factor and requires fewer correlations, which simplifies the original model greatly. Further studies, especially on separation- induced transition flows, are required for the improvement of the new model.

Application of focused laser differential interferometer to hypersonic boundary-layer instability study

Accurate prediction of hypersonic boundary-layer transition plays an important role in thermal protection system design of hypersonic vehicles.Restricted by the capability of spatial diagnostics for hypersonic boundary-layer study,quite a lot of problems of hypersonic boundary-layer transition,such as nonlinearity and receptivity,remain outstanding.This work reports the application of focused laser differential interferometer to instability wave development across hypersonic boundary-layer on a flared cone model.To begin with,the focused laser differential interferometer is designed and set up in a Mach number 6 hypersonic quiet wind tunnel with the focal point in the laminar boundary-layer of a 5 degree half-angle flared cone model.Afterwards,instability experiments are carried out by traversing the focal point throughout the hypersonic boundary-layer and the density fluctuation along the boundary-layer profile is measured and analyzed.The results show that three types of instability waves ranging from 10 k Hz to over 1 MHz are co-existing in the hypersonic boundary-layer,indicating the powerful capability of focused laser differential interferometer in dynamic response resolution for instability wave study in hypersonic flow regime;furthermore,quantitative analyses including spectra and bicoherence analysis of instability waves throughout the hypersonic boundary-layer for both cold and heated cone models are performed.

Measuring boundary-layer height under clear and cloudy conditions using three instruments

Boundary-layer height （BLH） under clear, altostratus and low stratus cloud conditions were measured by GPS sounding, wind profiler radar, and micro-pulse lidar during the atmospheric radiation measurement experiment from Sep. to Dec. 2008 in Shouxian, Anhui, China. Results showed that during daytime or nighttime, regardless of cloud conditions, the GPS sounding was the most accurate method for measuring BLH. Unfortunately, because of the long time gap between launchings, sounding data did not capture the diurnal evolution of the BLH. Thus, wind profile radar emerged as a promising instrument for direct and continuous measurement of the mixing height during the daytime, accurately determining BLH using the structure parameter of the electromagnetic refractive index. However, during nighttime, radar was limited by weak signal extraction and did not work well for determining the BLH of the stable boundary layer, often recording the BLH of the residual layer. While micro-pulse lidar recorded the evolution of BLH, it overestimated the BLH of the stable boundary layer. This method also failed to work under cloudy conditions because of the influence of water vapor. Future work needs to develop a method to determine BLH that combines the complimentary features of all three algorithms.

Mach reflection in steady supersonic flow considering wedge boundary-layer correction

Mach reflection in steady supersonic flow is an important phenomenon having received extensive studies,among which simplified theoretical models to predict the size of Mach stem and other flow structure are of particular interest.Past efforts for such models were based on inviscid assumption while in real cases the flow is viscous.Here in this paper we consider the influence of wedge boundary layer on the Mach stem height.This is done by including a simplified boundary layer model into a recently published inviscid model.In this viscous model,the wedge angle and the trailing edge height,which control the Mach stem height,are replaced by their equivalent ones accounting for the displacement effect of the wedge boundary layer,with the boundary layer assumed to be laminar or fully turbulent.This viscous model is shown to compare well with numerical results by computational fluid dynamics and gives a Mach stem height as function of the Reynolds number and Mach number.It is shown that due to the viscous effect,the Mach stem height is increased,through increasing the effective wedge angle.

Series solution of non-similarity natural convection boundary-layer flows over permeable vertical surface

The non-similarity solution for natural convection from a permeable isothermal vertical wall is considered. The governing boundary-layer equations for non-similarity flow and temperature fields are solved using the homotopy analysis method. The homotopy-Pade’ technique is applied to accelerate the convergence of the homotopy-series solution. The influence of physical parameters on the non-similarity flows is investigated in detail. Different from the previous analytic results,the homotopy-series solutions are convergent and valid for all physical parameters in the whole domain 0 x 【 ∞ and 0 y 【 ∞.

Wing aeroelasticity analysis based on an integral boundary-layer method coupled with Euler solver

An interactive boundary-layer method, which solves the unsteady flow, is developed for aeroelastic computation in the time domain. The coupled method combines the Euler solver with the integral boundary-layer solver （Euler/BL） in a ＂semi-inverse＂ manner to compute flows with the inviscid and viscous interaction. Unsteady boundary conditions on moving surfaces are taken into account by utilizing the approximate small-perturbation method without moving the computational grids. The steady and unsteady flow calculations for the LANN wing are presented. The wing tip displacement of high Reynolds number aero-structural dynamics （HIRENASD） Project is simulated under different angles of attack. The flutter-boundary predictions for the AGARD 445.6 wing are provided. The results of the interactive boundary-layer method are compared with those of the Euler method and experimental data. The study shows that viscous effects are significant for these cases and the further data analysis confirms the validity and practicability of the coupled method.

THE ROLES OF BOUNDARY-LAYER DYNAMICS ON THE DEVELOPMENT OF MADDEN-JULIAN OSCILLATION

In this paper,a tropical atmospheric model of relevance to shorts-term climate variations(Wang and Li 1993)is util- ized for study of the development of Madden-Julian oscillation.The model contains an interactive process of boundary-layer Ekman convergence and precipitation heating.The model is solved by expanding dependent variables in terms of parabolic cylindrical functions in the meridional direction and truncating three meridional modes n=0,2,4 for equatorial symmetric solutions.The free wave solutions obtained under long-wave approximation are induced as a Kelvin wave and two Rossby waves.After considering the effect of boundary-layer dynamic process,the modified Kelvin wave becomes unstable in long-wave bands with a typical growth rate on an order of 10^(-6) s^(-1)and an eastward phase speed of 10 m s^(-1);the most unstable mode is wavenumber one.These theoretical results are consistent with the ob- served Madden-Julian oscillation in equatorial area.For the two modified Rossby waves,one with a smaller meridional scale(n=4)decays except for extra long-waves;the other with a larger meridional scale(n=2)grows in short-wave bands.This may be relevant to explaining the westward propagation of super cloud clusters in the Madden-Julian oscillation.The theory suggests that the boundary-layer dynamic process is an important mechanism in the develop- ment of the Madden-Julian oscillation.

MHD BOUNDARY-LAYER FLOW AND HEAT TRANSFER OVER PERMEABLE PLATE WITH CONVECTIVE SURFACE BOUNDARY CONDITION

A numerical analysis has been carried out to investigate the problem of magnetohydrodynamic(MHD)boundary-layer flow and heat transfer of a viscous incompressible fluid over a fixed plate.Convective surface boundary condition is taken into account for thermal boundary condition.A problem formulation is developed in the presence of thermal radiation,magnetic field and heat source/sink parameters.A similarity transformation is used to reduce the governing boundary-layer equations to couple higher-order nonlinear ordinary differential equations.These equations are numerically solved using Keller–Box method.The effect of the governing parameters such as radiation,Prandtl number,Hartman number,heat source/sink parameter on velocity and temperature profile is discussed and shown by plotting graphs.It is found that the temperature is an increasing function of convective parameter A,radiation and heat source parameters.Besides,the numerical results for the local skin friction coefficient and local Nusselt number are computed and presented in tabular form.Finally a comparison with a previously published results on a special case of the problem has done and shows excellent agreement.

Nonlinear Violence in Nonlinear Oscillations at High Energy

This paper focuses on the characteristics of solutions of nonlinear oscillatory systems in the limit of very high oscillation energy, E;specifically, systems, in which the nonlinear driving force grows with energy much faster for x(t) close to the turning point, a(E), than at any position, x(t), that is not too close to a(E). This behavior dominates important aspects of the solutions. It will be called “nonlinear violence”. In the vicinity of a turning point, the solution of a nonlinear oscillatory systems that is affected by nonlinear violence exhibits the characteristics of boundary-layer behavior (independently of whether the equation of motion of the system can or cannot be cast in the traditional form of a boundary-layer problem.): close to a(E), x(t) varies very rapidly over a short time interval (which vanishes for E → ∞). In traditional boundary layer systems this would be called the “inner” solution. Outside this interval, x(t) soon evolves into a moderate profile (e.g. linear in time, or constant)—the “outer” solution. In (1 + 1)-dimensional nonlinear energy-conserving oscillators, if the solution is reflection-invariant, nonlinear violence determines the characteristics of the whole solution. For large families of nonlinear oscillatory systems, as E → ∞, the solutions for x(t) tend to common, indistinguishable profiles, such as periodic saw-tooth profiles or step-functions. If such profiles are observed experimentally in high-energy oscillations, it may be difficult to decipher the dynamical equations that govern the motion. The solution of motion in a central field with a non-zero angular momentum exhibits extremely fast rotation around a turning point that is affected by nonlinear violence. This provides an example for the possibility of interesting phenomena in (1 + 2)-dimensional oscillatory systems.

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.

A Numerical Tackling on Sakiadis Flow with Thermal Radiation

Momentum and energy laminar boundary layers of an incompressible fluid with thermal radiation about a moving plate in a quiescent ambient fluid are investigated numerically. Also, it has been underlined that the analysis of the roles of both velocity and temperature gradient at infinity is of key relevance for our results.

Large Eddy Simulation of Turbulent Channel Flow with 3D Roughness Using a Roughness Element Model

Large eddy simulation of turbulent channel flow with dense and small 3D roughness elements is carried out using a roughness element model Profiles of mean Reynolds stress, mean velocity and rms velocity as well as turbulent structures near the wall are obtained. The shear stress in the rough wall is larger than that in the smooth wall side and the rough wall has a larger influence on the channel flow. Profiles of mean streamwise velocity near the wall have logarithmic velocity distributions for both smooth and roughness walls, while there is a velocity decrease for the rough wall due to larger fractional drag. All the three components of rms velocities in the rough wall region are larger than that in the smooth wall region, and the roughness elements on the wall increase turbulent intensity in all directions. The s^reak spacing and average diameter of near wall quasi-s^reamwise vortices increase with the presence of roughness elements on the wall and it is shown that the rough wall induces complex and strong streamwise vortices. Results of dense and small 3D roughness elements in both turbulent statistics and structure, obtained with a relatively simple method, are found to be comparable to related experiments.

Analytical solution to stagnation-point flow and heat transfer over a stretching sheet based on homotopy analysis

This paper is concerned with two-dimensional stagnation-point steady flow of an incompressible viscous fluid towards a stretching sheet whose velocity is proportional to the distance from the slit. The governing system of partial differential equations is first transformed into a system of dimensionless ordinary differential equations. Analytical solutions of the velocity distribution and dimensionless temperature profiles are obtained for different ratios of free stream velocity and stretching velocity, Prandtl number, Eckert number and dimensionality index in series forms using homotopy analysis method（HAM）. It is shown that a boundary layer is formed when the free stream velocity exceeds the stretching velocity, and an inverted boundary layer is formed when the free stream velocity is less than the stretching velocity. Graphs are presented to show the effects of different parameters.

Studies on the Evaporation Regulation Mechanisms of Crude Oil and Petroleum Products

Various concepts for oil evaporation prediction are summarized. Models can be divided into those models that use the basis of air-boundary-regulation or those that do not. Experiments were conducted to determine if oil and petroleum evaporation is regulated by the saturation of the air boundary layer. Experiments included the examination of the evaporation rate with and without wind, in which case it was found that evaporation rates were similar for all wind conditions and no-wind conditions. Experiments where the area and mass varied showed that boundary-layer regulation was not governing for petroleum products. Under all experimental and environmental conditions, oils or petroleum products were not found to be boundary-layer regulated. Experiments on the rate of evaporation of pure compounds showed that compounds larger than Decane were not boundary-layer regulated. Many oils and petroleum products contain few compounds smaller than decane, and this explains why their evaporation is not air boundary-layer limited. Comparison of the air saturation levels of various oils and petroleum products shows that the saturation concentration of water, which is strongly air boundary-regulated, is significantly less than that of several petroleum hydrocarbons. Lack of air boundary-layer regulation for oils is shown to be a result of both this higher saturation concentration as well as a low (below boundary-layer value) evaporation rate.