Study of incident shock interaction with a two-dimensional boundary layer

This paper presents the results of an experimental study of a boundary layer disturbed by an incident shock for parameters which are characteristic of problems for flow about blade profiles in the final stages of high-power steam turbines.     

Optical vortices in the Ginzburg–Landau equation with cubic–quintic nonlinearity

In a dissipative system with cubic–quintic nonlinearity, the curious evolution of optical vortex beams characterized by different topological charges (TCs) is simulated numerically and presented their evolution profiles. We find that new vortices will be induced during propagation, and the behavior of vortices, as affected by the TC and the number of beads of the incident beam, as well as its size, is also discussed. Common rules associated with the initial conditions coming from various incident beams are developed to determine the number of induced vortices and the corresponding rotation direction. Attributed to the nonlinearity, during propagation we see the beams slowly expand to induce new vortices, which commonly appear in oppositely charged pairs, while the net topological charge of the vortex is conserved. Our results not only deepen the understanding of optical vortices, but also widen their potential applications.     

Three dimensional shock wave/boundary layer interactions

An investigation into a three-dimensional, curved shock wave interacting with a three-dimensional, curved boundary layer on a slender body is presented. Three different nose profiles mounted on a cylindrical body were tested in a supersonic wind tunnel and numerically simulated by solving the Navier–Stokes equations. The conical and hemispherical nose profiles tested were found to generate shock waves of sufficient strength to separate the boundary layer on the cylinder, while the shock wave generated by the ogival profile did not separate the boundary layer. For the separated flow, separation was found to occur predominantly on the windward side of the cylinder with the lee-side remaining shielded from the direct impact of the incident shock wave. A thickening of the boundary layer on the lee-side of all the profiles was observed, and in the conical and hemispherical cases this leads to the re-formation of the incident shock wave some distance away from the surface of the cylinder. A complex reflection pattern off the shock wave/boundary layer interaction (SWBLI) was also identified for the separated flow cases. For comparative purposes, an inviscid simulation was performed using the hemispherical profile. Significant differences between the viscous and inviscid results were noted including the absence of a boundary layer leading to a simplified shock wave reflection pattern forming. The behaviour of the incident shock wave on the lee-side of the cylinder was also affected with the shock wave amalgamating on the surface of the cylinder instead of away from the surface as per the viscous case. Test data from the wind tunnel identified two separation lines present on the cylindrical surface of the hemispherical SWBLI generator. The pair of lines were not explicitly evident in the original CFD simulations run, but were later identified in a high-resolution simulation.     

Short crested wave–current forces around a large vertical circular cylinder

An analytical solution for the diffraction of short crested incident wave along positive x-axis direction on a large circular cylinder with uniform current is derived. The important influences of currents on wave frequency, water run-up, wave force, inertia and drag coefficients on the cylinder profiles are investigated for short-crested incident wave. Based on the numerical results, we find wave frequency of short crested wave system is affected by incident angle and the strength of the currents. The wave frequency increases or decreases with increasing current speed following or opposing wave propagating direction. It shows that the effects of current speeds, current directions on water run-up on the circular cylinder with different radius for different wave numbers are very conspicuous when the incident wave changes from long crested plane waves to short-crested waves. With the increase of current speed, the water run-up on the cylinder becomes more and more high, and will exceed that of long crested plane wave and short crested wave case without currents even though the current speed is small. The total wave loads, inertia coefficient and drag coefficient exerted on a cylinder with currents would be larger compared to the wave loads exerted pure short-crested waves. Therefore, ocean engineers should consider the short crested wave–current load on marine constructs carefully.     

应用单色光的三维激光测速系统

本文发展了一种采用单色激光器的三维激光测速系统。四束入射光布置成等腰直角三角形,用一个前透镜聚焦。这种光学布置光能利用率高,并能适用于近壁处流速测量。当实验模型具有多层不同性质介面时,入射光仍能保持相交。为了分离三维速度分量,采用了频率分离和偏振分离技术。特别考虑了光轴分量的测量精度,提出了提高测量精度的一些方法。应用此系统测量了斜切尾钝体后部的三维平均速度和湍流度分布。     

The use of driver inserts to reduce non-ideal pressure variations behind reflected shock waves

Non-ideal shock tube facility effects, such as incident shock wave attenuation, can cause variations in the pressure histories seen in reflected shock wave experiments. These variations can be reduced, and in some cases eliminated, by the use of driver inserts. Driver inserts, when designed properly, act as sources of expansion waves which can counteract or compensate for gradual increases in reflected shock pressure profiles. An algorithm for the design of these inserts is provided, and example pressure measurements are presented that demonstrate the success of this approach. When these driver inserts are employed, near- ideal, constant-volume performance in reflected shock wave experiments can be achieved, even at long test times. This near-ideal behavior simplifies the interpretation of shock tube chemical kinetics experiments, particularly in experiments which are highly sensitive to temperature and pressure changes, such as measurements of ignition delay time of exothermic reactions.     

On flows behind shock waves reflected from a solid wall

The flow fields associated with the interaction of a normal shock wave with a plane wall kept at a constant temperature were studied based on kinetic theory which can describe appropriately the shock structure and its reflection process. With the use of a difference scheme, the time developments of the distributions of the fluid dynamic quantities (velocity, temperature, pressure and number density of the gas) were obtained numerically from the BGK model of the Boltzmann equation subject to the condition of diffusive-reflection at the wall for several cases of incident Mach number:M ₁=1.2, 1.5, 2.0, 3.0, 4.0, 5.0 and 6.0. The reflection process of the shocks is shown explicitly together with the resulting formation of the flow fields as time goes on. The nonzero uniform velocity toward the wall occurring between the viscous boundary layer and the reflected shock wave is found to be fairly large, the magnitude of which is of the order of several percent of the velocity induced behind the incident shock, decreasing as the incident Mach number increases. It is also seen that a region of positive velocity (away from the wall) within the viscous boundary layer manifests itself in the immediate vicinity of the wall, which is distinct for larger incident Mach numbers. Some of the calculated density profiles are compared with available experimental data and also with numerical results based on the Navier-Stokes equations. The agreement between the three results is fairly good except in the region close to the wall, where the difference in the conditions of these studies and the inappropriateness of the Navier-Stokes equations manifest themselves greatly in the gas behavior.This article was processed using Springer-Verlag TEX Shock Waves macro package 1990.     

Free surface Stokes flows obstructed by multiple obstacles

Gravity‐driven Stokes flow down an inclined plane over and around multiple obstacles is considered. The flow problem is formulated in terms of a boundary integral equation and solved using the boundary element method. A Hermitian radial basis function (RBF) is used for the interpolation of the free surface, generation of the unit normal and curvature, and to prescribe the far‐field conditions. For flow over an obstacle, hemispheres are taken. For flow around an obstacle, circular cylinders are modelled and the contact angle condition on the obstacle/free surface intersection specified using the RBF formulation. Explicit profiles are produced for flow over and around two obstacles placed in various locations relative to one another. Interaction due to two obstacles is given by comparisons made with the profiles for flow over and around individual obstacles. In general, when the obstacles are separated by a sufficiently large distance the flow profiles are identical to a single obstacle analysis. For flow over and around two obstacles in‐line with the incident flow, effects of the governing parameters are examined, with variations in plane inclination angle, Bond number, obstacle size, and in the case of obstacles intersecting the free surface, static contact angle is considered. Finally flows over and around three obstacles are modelled. Copyright © 2009 John Wiley & Sons, Ltd.     

The interactions between wave-currents and offshore structures with consideration of fluid viscosity

Study of the flow field around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the flow field around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averagedk-ɛ turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds numberRe=10⁵, a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies. The project supported by the National Natural Science Foundation of China and Foundation of State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University.     

Shock-induced mixing zone characterization: an attempt by hot-wire diagnostic

An attempt to extract velocity and molar fraction from a single hot-wire trace within a turbulent mixing zone induced by a shock accelerated gaseous interface has been proposed. Experiments have been conducted for negative and positive density jumps across the interface. The hot-wire signals clearly show interfaces between mixed and unmixed regions and the locations of incident and reflected shocks. With some hypotheses on the temperature, velocity and molar fraction profiles within the turbulent mixing zone have been obtained solving an inverse problem. Results show that if the molar fraction profiles follow physically coherent evolutions, those of the local velocity are strongly correlated with the choice of its variation range. So, we reasonably think that the results obtained from single wire have to remain limited to interface and shock locations. And it is only by coupling the present technique with the laser Doppler velocimetry, which we will be able to possibly obtain reliable estimates of the variations of quantities in the turbulent mixing zone.      

The blast resistance of sandwich composites with stepwise graded cores

Shock tube experiments were performed to study the dynamic response of sandwich panels with E-Glass Vinyl Ester (EVE) composite face sheets and stepwise graded styrene foam cores. Two types of core configurations, with identical areal density, were subjected to the shock wave loading. The core layers were arranged according to the density of the respective foam; configuration 1 consisted of low/middle/high density foams and configuration 2 consisted of middle/low/high density foams. The method to calculate the incident and reflected energies of the shock wave, as well as the deformation energy of the specimen, were proposed based on the shock wave pressure profiles and the high speed deflection images that were obtained. The experimental results showed that configuration 1 outperformed configuration 2 in regards to their blast resistance. Significant core material compression was observed in configuration 1, while in configuration 2 the core layers disintegrated and the front skin (blast side) fractured into two pieces along the midsection. The estimated energies were then calculated for both configurations. The total energy difference between the incident and reflected energies was almost identical, even though the deformation energy for configuration 2 was larger.     

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

Free-field wave solutions in a half-plane exhibiting a special-type of continuous inhomogeneity

This work presents closed-form solutions for free-field motions in a continuously inhomogeneous half-plane that include contributions of incident waves as well as of waves reflected from the traction-free horizontal surface. Both pressure and vertically polarized shear waves are considered. Furthermore, two special types of material inhomogeneity are studied, namely (a) a shear modulus that varies quadratically with respect to the depth coordinate and (b) one that varies exponentially with the said coordinate. In all cases, Poisson’s ratio is fixed at one-quarter, while both shear modulus and material density profiles vary proportionally. Next, a series of numerical results serve to validate the aforementioned models, and to show the differences in the wave motion patterns developing in media that are inhomogeneous as compared to a reference homogeneous background. These results clearly show the influence of inhomogeneity, as summarized by a single material parameter, on the free-field motions that develop in the half-plane. It is believed that this type of information is useful within the context of wave propagation studies in non-homogeneous continua, which in turn find applications in fields as diverse as laminated composites, geophysical prospecting, oil exploration and earthquake engineering.     

Loading and unloading split hopkinson pressure bar pulse-shaping techniques for dynamic hysteretic loops

Pulse-shaping techniques are developed for both the loading and unloading paths of a split Hopkinson pressure bar (SHPB) experiment to obtain valid dynamic stress-strain loops for engineering materials. Front and rear pulse-shapers, in association with a momentum trap, are used to precisely control the profiles of the loading and unloading portions of the incident pulse. The modifications, ensure that the specimen deforms at the same constant strain rate under dynamic stress equilibrium during both loading and unloading stages of an experiment so that dynamic stress-strain loops can be accurately determined. Dynamic stress-strain loops with a constant strain rate for a nickel-titanium shape memory alloy and polymethyl methacrylate are determined using the modified SHPB. The modified momentum trap prevents repeated loading on a specimen without affecting the amplitude of the desired loading pulse and without damaging the bar at high stress levels.     

连续分层流体中垂直薄板的水动力特性

研究了在线性连续分层流体中水波与半潜式刚性垂直薄板相互作用的问题. 在 Boussinesq近似下，基于分离变量法，导出了具有自由面的平面前进波的色散关系，建立 了半潜式刚性垂直薄板的反射与透射能量、水平波浪力的计算方法. 对给定的频率，当它大 于浮力频率时，流场中只有一种模态的平面前进波，当它小于浮力频率时，流场 中有无数多个模态的平面前进波，并证明了对每一种模态的入射波，其它每个模态水波的反射与透射能量是 相等的. 对水面漂浮和座底半潜式薄板的反射与透射能量，以及作用在薄板上的水平波浪力 进行了数值计算分析，表明了流体的线性连续分层效应对这些水动力的影响是不可忽视 的. 特别地，在入射波频率小于浮力频率时，与第1模态入射波的能量转化量及其对薄板产 生的水平波浪力相比，其它模态入射波的能量转化量及其对薄板产生的水平波浪力都要大得 多.     

Effect of surroundings on the transient energy transfer in a layer of radiating gas

Summary An analysis is presented for the transient cooling or heating of a stagnant layer of hot radiating gas surrounded by a cold gas capable of absorbing and emitting radiation. Scattering of radiation is neglected, and energy transfer by conduction and convection is considered to be negligible compared with radiation. The gas is assumed to be perfect and in local thermodynamic equilibrium. The heating of a cold gas by a diffuse and a collimated radiation flux incident on the boundary of the gas from some external source is considered, and the dependence of physical and radiative properties on temperature is taken into account. The problem is formulated exactly using radiative transfer theory. A scheme is developed for the numerical solution of the nonlinear integrodifferential equations of energy conservation. Starting with arbitrary, but given, initial temperature profiles, temperature distributions and local radiative fluxes are predicted as a function of time for a wide range of physically interesting conditions.     

Flow regimes and frequency selection of a cylinder oscillating in an upstream cylinder wake

This paper describes flow around a pair of cylinders in tandem arrangement with a downstream cylinder being fixed or forced to oscillate transversely. A sinusoidal parietal velocity is applied to simulate cylinder oscillation. Time-dependent Navier-Stokes equations are solved using finite element method. It is shown that there exist two distinct flow regimes: ‘vortex suppression regime’ and ‘vortex formation regime’. Averaged vortex lengths between the two cylinders, pressure variations at back and front stagnant points as well as circumferential pressure profiles of the downstream cylinder are found completely different in the two regimes and, thus, can be used to identify the flow regimes. It is shown that frequency selection in the wake of the oscillating cylinder is a result of non-linear interaction among vortex wakes upstream and downstream of the second cylinder and its forced oscillation. Increasing cylinder spacing results in a stronger oscillatory incident flow upstream of the second cylinder and, thus, a smaller synchronization zone.     

The interaction of a shock wave with the boundary layer in a reflected shock tunnel

The influence of a nontotal reflection on the interaction of a reflected shock wave with the boundary layer in a reflected shock tunnel has been investigated. The calculating method of the velocity, the temperature and the Mach number profiles in the boundary layer in reflected shock fixed coordinates has been obtained. To account for equilibrium real gas effects of nitrogen, the numerical results show that the minimum Mach number in the boundary layer has been moved from the wall into the boundary layer with the increasing of the incident shock Mach number. The minimum Mach number, the shock angle in the bifurcated foot and the jet velocity along the wall to the end plate are reduced owing to the increasing of the area of nozzle throat. The numerical results are in good agreement with measurements.     

Reflection and transmission of regular water waves by a thin,floating plate

Measurements of the wave fields reflected and transmitted by a thin floating plastic plate are reported for regular incident waves over a range of incident periods (producing wavelengths comparable to the plate length) and steepnesses (ranging from mild to storm-like). Two different plastics are tested, with different densities and mechanical properties, and three different configurations are tested. The configurations include freely floating plates, loosely moored plates (to restrict drift), and plates with edge barriers (to restrict waves overwashing the plates). The wave fields reflected and transmitted by plates without barriers are shown to become irregular, as the incident waves become steeper, particularly for the denser plastic and the moored plate. Further, the proportion of energy transmitted by the plates without barriers is shown to decrease as the incident wave becomes steeper, and this is related to wave energy dissipation.     

Lamb wave scattering by a symmetric pair of surface-breaking cracks in a plate

The scattering problem of a Lamb wave incident on a symmetric pair of surface-breaking transverse cracks in a plate is considered. The Lamb wave is assumed to be obliquely incident on the crack plane. Since the cracks are part-through, the scattered field will contain reflected as well as transmitted waves. The energy of the incoming wave is partitioned into reflected and transmitted wave modes. Energy coefficients of the reflected and transmitted waves are calculated as a function of incident frequency and crack depth. The incidence angle of the incoming wave is also treated as a parameter. Both the reflected and transmitted wave fields are considered as linear superpositions of all real and complex wave modes in the plate. Decomposition of modes is achieved with the help of an orthogonality condition based on the principle of reciprocal work. Continuity of displacement and stress fields is imposed at the crack plane. Energy coefficients for reflection and transmission are obtained from the mode amplitudes. Energy coefficients are shown to be a strong function of incident frequency and crack depth. Experiments are conducted with a PZT transducer network interacting with a symmetric pair of machined cracks in an aluminum plate. Trends predicted by the analysis are reflected in the experimental results.