基于菲索干涉仪的边界层测风激光雷达研究

分析了多光束菲索(Fizeau)干涉仪的光学特性，以及影响干涉光谱性能的参数，并研究了激光雷达系统接收信号能力，提出了一套基于Fizeau干涉仪和CCD探测器的边界层测风激光雷达系统。利用标准大气参数和切合实际的激光雷达系统参数，对干涉仪进行优化设计模拟计算，得到了合理的参数，满足边界层1m／s风速测量精度要求。

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.

太阳耀斑和日食期间电离层行为的统计和模拟研究

The ionosphere is the ionized component of the Earth's upper atmosphere. Solar EUV radiation is the source of ionospheric ionization. Thus the ionosphere is affected strongly by the variations in solar radiation. Solar flares and solar eclipses can induce remarkable short time changes in solar radiation: the solar radiation would increase suddenly during solar flares and decrease significantly during solar eclipses. Solar flare and eclipse events not only affect directly the photochemical processes, but also affect the dynamic processes, and even affect the neutral atmosphere, which is strongly coupled with the ionosphere. The study on the ionospheric response to solar flares and eclipses can advance our knowledge on the ionosphere and its photochemical and dynamic processes and help us to evaluate the ionospheric parameters (such as ion loss coefficients). In addition, the study on the ionospheric responses to solar flares and eclipses is an important part of the ionospheric space weather, which can provide guides for space weather monitoring. This thesis devotes to the study on the ionospheric responses to solar flares and solar eclipses. I have developed two models to simulate the variations of solar EUV radiation during solar flares and solar eclipses, and involved in developing a 2D mid- and low-latitude ionospheric model. On the basis of some observed data and the ionospheric model, I study the temporal and spatial variations of the ionosphere during solar flares and eclipses, and investigate the influences of solar activity, solar zenith angle, neutral gas density, and magnetic dip angle on the ionospheric responses to solar flares and solar eclipses. The main points of my works and results are summarized as follows. 1. The ionospheric response to the X17.2 solar flare on October 28, 2003 was modeled via using a one-dimension theoretical ionospheric model. The simulated variation of TEC is in accordance with the observations, though there are some differences in the amplitude of the variation. Then I carried out a series of simulations to explore the local time and seasonal dependences of the ionospheric responses to solar flares. These calculations show that the ionospheric responses are largely related with the solar zenith angle (SZA). During the daytime (small SZA), most of the increases in electron density occur at altitudes below 300 km with a peak at around 115 km; whereas around sunrise and sunset (SZA>90°), the strongest ionospheric responses occur at much higher altitudes. The TEC increases slower at sunrise than at sunset, which is caused by the difference in the evolution of SZA at sunrise and sunset: SZA decreases with time at sunrise and increase with time at sunset. The ionospheric response is largest in summer and smallest in winter, which is also related to the seasonal difference of SZA. 2. Based on the observations from the ionosondes in Europe and the ionospheric model, I investigated the differences of the ionosphere responses to solar eclipses between the E-layer and F1-layer. Both the observation and simulation show that the decrease in foF1 due to the solar eclipses is larger than that in foE. This effect is due to that the F1 region locates at the transition height between the atomic ion layer and the molecular ion layer. With the revised model of solar radiation during solar flares, our model calculates the radiations from both the inside and outside of photosphere. Large discrepancy can be found between the observations and the calculations with an unrevised model, while the calculations with the revised model consist with the observations. 3. I also explore the effects of the F2-layer height, local time, solar cycle, and magnetic dip angle on the ionospheric responses to solar eclipses via using an ionospheric model and study on the solar zenith angle and the dip dependences by analyzing the data derived from 23 ionosonde stations during seven eclipse events. Both the measured and simulated results show that these factors have significant effect on the ionospheric response. The larger F2-layer height causes the smaller decrease in foF2, which is because that the electron density response decreases with height. The larger dip results in the smaller eclipse effect on the F2 layer, because the larger dip would cause the more diffusion from the top ionosphere which can make up for the plasma loss. The foF2 response is largest at midday and decreases with the increasing SZA. The foF2 response is larger at high solar activity than at low solar activity. The simulated results show that the local time and solar activity discrepancy of the eclipse effect mainly attribute to the difference of the background neutral gas density. 4. I carried out a statistical study on the latitudinal dependence of the ionospheric response to solar eclipses and modeled this latitudinal dependence by the ionospheric model. Both the observations and simulations show that the foF2 and TEC responses have the same latitudinal dependence: the eclipse effects on foF2 and TEC are smaller at low latitudes than at middle latitudes; at the middle latitudes (>40°), the eclipse effect decreases with increasing latitude. In addition, the simulated results show the change in electron temperature at the heights of above 300 km of low latitudes is much smaller than that at the same heights of middle latitudes. This is due to the smaller decrease in photoelectron production rate at its conjugate low heights. 5. By analyzing the observed data during the October 3, 2005 solar eclipse, I find some significant disturbances in the conjugate region of the eclipse region, including a decrease in Te, an increase in foF2 and TEC, and an uprising in hmF2. I also simulated the ionosphere behavior during this eclipse using a mid-low latitude ionospheric model. The simulations reproduce the measured ionospheric disturbances mentioned above in the conjugated hemisphere. The simulations show that the great loss of arriving photoelectron heat from the eclipse region is the principal driving source for the disturbances in the conjugate hemisphere.