Mechanisms of the sidewall facet evolution in lateral epitaxial overgrowth of GaN by MOCVD

The lateral epitaxial overgrowth of GaN was carried out by low-pressure metalorganic chemical vapor deposition, and the cross section shape of the stripes was characterized by scanning electron microscopy. Inclined {11-2n} facets (n approximate to 1-2.5) were observed in the initial growth, and they changed gradually into the vertical {11-20} sidewalls in accordance with the process of the lateral overgrowth. A model was proposed utilizing diffusion equations and boundary conditions to simulate the concentration of the Ga species constituent throughout the concentration boundary layer. Solutions to these equations are found using the two-dimensional, finite element method. We suggest that the observed evolution of sidewall facets results from the variation of the local V/III ratio during the process of lateral overgrowth induced by the lateral supply of the Ga species from the SiNx mask regions to the growing GaN regions.

Annealing behavior of hexagonal phase content in cubic GaN thin films grown on GaAs (001) by MOCVD

The annealing behavior of the hexagonal phase content in cubic GaN (c-GaN) thin films grown on GaAs (001) by MOCVD is reported. C-GaN thin films are grown on GaAs (001) substrates by metalorganic chemical vapor deposition (MOCVD). High temperature annealing is employed to treat the as-grown c-GaN thin films. The characterization of the c-GaN films is investigated by photoluminescence (PL) and Raman scattering spectroscopy. The change conditions of the hexagonal phase content in the metastable c-GaN are reported. There is a boundary layer existing in the c-GaN/GaAs film. When being annealed at high temperature, the intensity of the TOB and LOB phonon modes from the boundary layer weakens while that of the E-2 phonon mode from the hexagonal phase increases. The content change of hexagonal phase has closer relationship with annealing temperature than with annealing time period.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

Performance of supersonic model combustors with staged injections of supercritical aviation kerosene

Supersonic model combustors using two-stage injections of supercritical kerosene were experimentally investigated in both Mach 2.5 and 3.0 model combustors with stagnation temperatures of approximately 1,750 K. Supercritical kerosene of approximately 760 K was prepared and injected in the overall equivalence ratio range of 0.5-1.46. Two pairs of integrated injector/flameholder cavity modules in tandem were used to facilitate fuel-air mixing and stable combustion. For single-stage fuel injection at an upstream location, it was found that the boundary layer separation could propagate into the isolator with increasing fuel equivalence ratio due to excessive local heat release, which in turns changed the entry airflow conditions. Moving the fuel injection to a further downstream location could alleviate the problem, while it would result in a decrease in combustion efficiency due to shorter fuel residence time. With two-stage fuel injections the overall combustor performance was shown to be improved and kerosene injections at fuel rich conditions could be reached without the upstream propagation of the boundary layer separation into the isolator. Furthermore, effects of the entry Mach number and pilot hydrogen on combustion performance were also studied.

MHD控制超声速边界层的理论研究和数值分析

对MHD(mechanisms of magnetohy drodynamics)控制超声速平板湍流边界层的机理进行了理论研究和数值模拟.理论上,采用等离子体低频近似碰撞频率模型,建立等离子体中电子和离子的力平衡方程,得到等离子体速度、极化电场以及边界层速度.数值上,通过空间HLLE格式、LU-SGS时间推进求解时均磁流体动力学湍流方程,其中湍流模型采用sst-kω双方程模型.研究结果表明:(1)边界层速度的理论结果和数值结果误差在7%范围内;(2)只有磁场而电场为零时,洛仑兹力起到减小摩阻的作用.施加电场后,洛仑兹力能够加速边界层低速区流体;(3)在边界层外层,越靠近壁面,作用参数越小;而在边界层近壁区黏性底层,虽然惯性力减小,但黏性力却迅速增加,因此越靠近壁面,作用参数反而越大,加速低速流的代价增加.

MHD控制激波诱导湍流边界层分离的机理分析

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程.其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解.计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除.

Diagnostic experiments for transport mechanisms of suspended sediment discharged from the Yellow River in the Bohai Sea

Five diagnostic experiments with a 3D baroclinic hydrodynamic and sediment transport model ECOMSED in couple with the third generation wave model SWAN and the Grant-Madsen bottom boundary layer model driven by the monthly sediment load of the Yellow River, were conducted to separately diagnose effects of different hydrodynamic factors on transport of suspended sediment discharged from the Yellow River in the Bohai Sea. Both transport and spatio-temporal distribution of suspended sediment concentration in the Bohai Sea were numerially simulated. It could be concluded that suspended sediment discharged from the Yellow River cannot be delivered in long distance under the condition of tidal current. Almost all of sediments from the Yellow River are deposited outside the delta under the condition of wind-driven current, and only very small of them are transported faraway. On the basis of wind forcing, sediments from the Yellow River are mainly transported north-northwestward, and others which are first delivered to the Laizhou Bay are continuously moved northward. An obvious 3D structure characteristic of sediment transport is produced in the wind-driven and tide-induced residual circulation condition. Transport patterns at all layers are generally consistent with circulation structure, but there is apparent deviation between the depth-averaged sediment flux and the circulation structure. The phase of temporal variation of sediment concentration is consistent with that of the bottom shear stress, both of which are proved to have a ten-day cycle in wave and current condition.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.

Characteristics of low-frequency components of the near-bottom current in the Chinese Pioneer Area

Two deep-sea moorings were deployed respectively in the east area and the west area of Chinese Pioneer Area (CPA) in the tropic east Pacific to monitor the regional deep-sea dynamics below 600 meters above bottom (mab) from July 1997 to Oct. 1999. Results of statistics, spectral estimate and correlation analysis of the low-passed velocity data show that time scales of low-frequency components of the near-bottom currents are 25similar to120 days, in which 51-day period dominates the lower band of the frequency domain. Topographic features have obvious effect on low-frequency currents below 50 mab; modulations of the bottom-intensified sheared mean flow to the low-frequency currents are the dynamic mechanism of the frequency shift that occurs in both the east-area and the west-area.

Ship-board measurements and estimations of air-sea fluxes in the western tropical Pacific during TOGA COARE

Direct air-sea flux measurements were made on RN Kexue #1 at 40 degrees S, 156 degrees E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-alpha hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme, There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results, Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heal flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.

STUDY ON THE FORMATIONAL MECHANISM OF THE NORTHERN YELLOW (HUANGHAI) SEA COLD WATER MASS(I)

A hydrodynamic-thermodynamic equation set was set up to reflect the formational mechanism and evolution of the Northern Yellow (Huanghai) Sea cold water mass (NYSCWM) and its density circulation. Appropriate mathematical physical models were established by using some physical postulations. An approximate analytic solution to expound the distributions of temperature and three-dimensional current velocity, which can be used to expound the formational mechanism of the NYSCWM and its density circulation is obtained by using the theory of boundary layer and perturbational analyses.

Fine sediment resuspension dynamics in Moreton Bay

A comprehensive field study has been undertaken to investigate sediment resuspension dynamics in the Moreton Bay, a large semi-enclosed bay situated in South East Queensland, Australia. An instrumented tripod, which housed three current meters, three OBS sensors and one underwater video camera, was used to collect the field data on tides, currents, waves and suspended sediment concentrations at four sites (Sites 1, 2, 4, and 5) in the bay. Site I was located at the main entrance, Site 2 at the central bay in deep water, and Sites 4 and 5 at two small bays in shallow water. The bed sediment was fine sand (d(50) = 0. 2 mm) at Site 1, and cohesive sediment at the other three sites. Based on the collected field data, it is found that the dominant driving forces for sediment resuspension are a combination of ocean swell and tidal currents at Site 1, tidal currents at Site 2, and wind-waves at Sites 4 and 5. The critical bed shear stress for cohesive sediment resuspension is determined as 0. 079 Pa in unidirectional flow at Site 2, and 0. 076 Pa in wave-induced oscillatory flow at Site 5.

Wind-vector estimation for RADARSAT-1 SAR images: Validation of wind-direction estimates based upon geometry diversity

In this letter, a new wind-vector algorithm is presented that uses radar backscatter sigma(0) measurements at two adjacent subscenes of RADARSAT-1 synthetic aperture radar (SAR) images, with each subscene having slightly different geometry. Resultant wind vectors are validated using in situ buoy measurements and compared with wind vectors determined from a hybrid wind-retrieval model using wind directions determined by spectral analysis of wind-induced image streaks and observed by colocated QuikSCAT measurements. The hybrid wind-retrieval model consists of CMOD-IFR2 [applicable to C-band vertical-vertical (W) polarization] and a C-band copolarization ratio according to Kirchhoff scattering. The new algorithm displays improved skill in wind-vector estimation for RADARSAT-1 SAR data when compared to conventional wind-retrieval methodology. In addition, unlike conventional methods, the present method is applicable to RADARSAT-1 images both with and without visible streaks. However, this method requires ancillary data such as buoy measurements to resolve the ambiguity in retrieved wind direction.

Tropospheric aerosol detection over Southern Ireland using a backscatter Raman lidar

Lidar is an optical remote sensing instrument that can measure atmospheric parameters. A Raman lidar instrument (UCLID) was established at University College Cork to contribute to the European lidar network, EARLINET. System performance tests were carried out to ensure strict data quality assurance for submission to the EARLINET database. Procedures include: overlap correction, telecover test, Rayleigh test and zero bin test. Raman backscatter coefficients, extinction coefficients and lidar ratio were measured from April 2010 to May 2011 and February 2012 to June 2012. Statistical analysis of the profiles over these periods provided new information about the typical atmospheric scenarios over Southern Ireland in terms of aerosol load in the lower troposphere, the planetary boundary layer (PBL) height, aerosol optical density (AOD) at 532 nm and lidar ratio values. The arithmetic average of the PBL height was found to be 608 ± 138 m with a median of 615 m, while average AOD at 532 nm for clean marine air masses was 0.119 ± 0.023 and for polluted air masses was 0.170 ± 0.036. The lidar ratio showed a seasonal dependence with lower values found in winter and autumn (20 ± 5 sr) and higher during spring and winter (30 ± 12 sr). Detection of volcanic particles from the eruption of the volcano Eyjafjallajökull in Iceland was measured between 21 April and 7 May 2010. The backscatter coefficient of the ash layer varied between 2.5 Mm-1sr-1 and 3.5 Mm-1sr-1, and estimation of the AOD at 532 nm was found to be between 0.090 and 0.215. Several aerosol loads due to Saharan dust particles were detected in Spring 2011 and 2012. Lidar ratio of the dust layers were determine to be between 45 and 77 sr and AOD at 532 nm during the dust events range between 0.84 to 0.494.

Uniformly convergent finite element methods for singularly perturbed parabolic partial differential equations

This thesis is concerned with uniformly convergent finite element methods for numerically solving singularly perturbed parabolic partial differential equations in one space variable. First, we use Petrov-Galerkin finite element methods to generate three schemes for such problems, each of these schemes uses exponentially fitted elements in space. Two of them are lumped and the other is non-lumped. On meshes which are either arbitrary or slightly restricted, we derive global energy norm and L2 norm error bounds, uniformly in the diffusion parameter. Under some reasonable global assumptions together with realistic local assumptions on the solution and its derivatives, we prove that these exponentially fitted schemes are locally uniformly convergent, with order one, in a discrete L∞norm both outside and inside the boundary layer. We next analyse a streamline diffusion scheme on a Shishkin mesh for a model singularly perturbed parabolic partial differential equation. The method with piecewise linear space-time elements is shown, under reasonable assumptions on the solution, to be convergent, independently of the diffusion parameter, with a pointwise accuracy of almost order 5/4 outside layers and almost order 3/4 inside the boundary layer. Numerical results for the above schemes are presented. Finally, we examine a cell vertex finite volume method which is applied to a model time-dependent convection-diffusion problem. Local errors away from all layers are obtained in the l2 seminorm by using techniques from finite element analysis.