Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Analysis of critical excavation depth for a jointed rock slope using a face-to-face discrete element method

The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D-s and the natural slope angle alpha, the joint inclination angle theta as well as the strength parameters of the joints c(r) ,phi(r) is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Distribution, stability, bifurcations and catastrophe of steady rotation of a symmetrical heavy gyroscope with viscous-liquid-filled cavity

The number, the angles of orientation and the stability in Rumyantsev Movchan's sense of oblique steady rotations of a symmetric heavy gyroscope with a cavity completely filled with a uniform viscous liquid, possessing a fixed point 0 on its symmetric axis. are given for various values of the parameters. By taking the square of the upright component of the angular momentum M2 as a control parameter, three types of bifurcation diagrams of the steady rotations, two types of jumps and two kinds of local catastrophes, one being the symmetric reduced cusp type and the other being of the symmetric reduced butterfly type, are obtained. By taking account of the M2-damping owing to the moment of unavoidable faint friction, two different modes for the gyroscope, initially in a stable quasi-steady upright rotation with a nutation angle theta(s) equal to zero, to topple over are found.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Experimental Study and Numerical Simulation of Cellular Structures and Mach Reflection of Gaseous Detonation Wave

In this paper the Deflagration to Detonation Transition (DDT) process of gaseous H-2-O-2 mixture and Mach reflection of gaseous detonation wave on a wedge have been conducted experimentally. The cellular pattern of DDT process and Mach reflection were obtained from experiments with wedge angle theta = 10(0) similar to 40(0) and initial pressure of gaseous mixture 16kPa similar to 26.7kPa. The 2-D numerical simulations of DDT process and Mach reflection of detonation wave were performed by using the simplified ZND model and improved space-time conservation element and solution element (CE/SE) method. The numerical cellular structures were compared with the cellular patterns of soot track. Compared results were shown that it is satisfactory. The characteristic comparisons on Mach reflection of air shock wave and detonation wave were carried also out and their differences were given.

On the contact width in generalized plain strain JKR model for isotropic solids

Adhesive contact model between an elastic cylinder and an elastic half space is studied in the present paper, in which an external pulling force is acted on the above cylinder with an arbitrary direction and the contact width is assumed to be asymmetric with respect to the structure. Solutions to the asymmetric model are obtained and the effect of the asymmetric contact width on the whole pulling process is mainly discussed. It is found that the smaller the absolute value of Dundurs' parameter beta or the larger the pulling angle theta, the more reasonable the symmetric model would be to approximate the asymmetric one.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ &theta; ≤ 30° (&theta; = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ &theta; ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ &theta; ≤ 30° (&theta; = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ &theta; ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Angular distribution and conversion of multi-keV L-shell X-ray sources produced from nanosecond laser irradiated thick-foil targets

An experimental study on the angular distribution and conversion of multi-keV X-ray sources produced from 2 ns-duration 527nm laser irradiated thick-foil targets on Shenguang II laser facility (SG-II) is reported. The angular distributions measured in front of the targets can be fitted with the function of f(theta) = alpha+ (1- alpha)cos(beta) theta (theta is the viewing angle relative to the target normal), where alpha = 0.41 +/- 0.014, beta = 0.77 +/- 0.04 for Ti K-shell X-ray Sources (similar to 4.75 keV for Ti K-shell), and alpha = 0.085 +/- 0.06, beta = 0.59 +/- 0.07 for Ag/Pd/Mo L-shell X-ray Sources (2-2.8 keV for Mo L-shell, 2.8-3.5 keV for Pd L-shell, and 3-3.8 keV for Ag L-shell). The isotropy of the angular-distribution of L-shell emission is worse than that of the K-shell emission at larger viewing angle (>70 degrees), due to its larger optical depth (stronger self-absorption) in the cold plasma side lobe Surrounding the central emission region, and in the central hot plasma region (emission region). There is no observable difference in the angular distributions of the L-shell X-ray emission among Ag, Pd, and Mo. The conversion efficiency of Ag/Pd/Mo L-shell X-ray sources is higher than that of the Ti K-shell X-ray sources, but the gain relative to the K-shell emission is not as high as that by using short pulse lasers. The conversion efficiency of the L-shell X-ray sources decrease, with increasing atomic numbers (or X-ray photon energy), similar to the behavior of the K-shell X-ray Source.

大气湍流对斜程传输准单色高斯-谢尔光束空间相干性的影响

由湍流大气中斜程传输时准单色高斯-谢尔(GSM)光束互相干函数的解析式导出了该光束的复相干度.然后,利用表征光束空间相干性的横向相干长度,研究了斜程传输时大气湍流对准单色GSM光束空间相干性的影响.研究结果表明:1)当传输路径偏离水平方向较大(即θ≤88°)时,准单色GSM光束横向相干长度随传输距离均为先迅速增加,后缓慢增加,最后基本保持不变.2)当传输路径接近水平方向(即θ≥89°)时,准单色GSM光束横向相干长度随传输距离均为先增大,达到一个最大值后开始下降并持续减小.3)空间相干性不同的准单色GSM光束在湍流大气中传输一定距离后空间相干性趋于相同.最后,对以上结果在物理上给予了解释.

Super-resolution with annular binary phase filter in the 4Pi-confocal system

In order to realize super-resolution in the 4Pi-confocal systems, the annular binary pure phase filter is designed with the vector diffraction theory. The relations between the super-resolved parameters, such as S, G(T), G(A), and the radial position theta(i) of each zone, are obtained. For simple illumination of the design procedure, three-zone binary pure phase filters are studied, and several numerical simulation results show that in the 4Pi-confocal system with the properly designed binary pure phase filter the super-resolution can be realized with low sidelobes.

高功率钕玻璃激光三倍频脉冲时间波形的研究

基于目前国内规模最大的激光驱动器——“神光Ⅱ”八路基频光已经实现功率平衡运行，通过改变其中若干路三倍频系统各调谐量的偏离，对输出三倍频波形进行束与束之间的横向对比研究．研究发现，对于Ⅱ类-Ⅱ类偏振失配三倍频系统，在影响转换效率的三个调谐量中，偏振分配角失配△θp，对三倍频波形影响最大；在入射基频功率密度约为1．0GW／cm^2情况下，当三倍频系统三个调谐量都处在最佳匹配时，三倍频波形半峰全宽τ最小。研究工作为最终实现“神光Ⅱ”八路光束三倍频功率平衡输出提供了晶体调试的方法。

偏振分配角θp对高功率激光三倍频波形的影响

使用四阶龙格-库塔算法数值求解Ⅱ类-Ⅱ类偏振失配方案的三倍频稳态耦合波方程组。针对“神光Ⅱ”激光驱动器的三倍频系统，在二倍频匹配角吼和三倍频匹配角只均处于最佳匹配情况下，以1°为单位改变偏振分配角郇的偏离量，得到对应的三倍频波形特性的变化。研究表明，偏振分配角θp处于最佳理论匹配位时，对应的三倍频波形半峰全宽最窄，且波形上升沿抬起也最快；θp从最佳匹配位增加比从最佳匹配位减小对应的半峰全宽变宽更快，对应的上升沿抬起更慢。根据所得理论结果，在实验中调节偏振分配角θp，达到了改变三倍频输出波形特性，进而实现八

Strong visible up-conversion emissions and thermometric applications of Er3+-Yb3+ codoped Al2O3 prepared by the sol-gel method

The Er3+-Yb3+ codoped Al2O3 has been prepared by the sol-gel method using the aluminium isopropoxide [Al(OC3H7)(3)]-derived Al2O3 sols with addition of the erbium nitrate [Er(NO3)(3) center dot 5H(2)O] and ytterbium nitrate [Yb(NO3)(3) center dot 5H(2)O]. The phase structure, including only two crystalline types of doped Al2O3 phases, theta and gamma, was obtained for the 1 mol% Er3+ and 5 mol% Yb3+ codoped Al2O3 at the sintering temperature of 1,273 K. By a 978 nm semiconductor laser diodes excitation, the visible up-conversion emissions centered at about 523, 545, and 660 nm were obtained. The temperature dependence of the green up-conversion emissions was studied over a wide temperature range of 300-825 K, and the reasonable agreement between the calculated temperature by the fluorescence intensity ratio (FIR) theory and the measured temperature proved that Er3+-Yb3+ codoped Al2O3 plays an important role in the application of high temperature sensor.