Investigation on the nonlinear optical properties of the DCNP/PEK-c guest-host polymer films

The polyetherketone (PEK-c) guest-host polymer thin films doped with 3-(1,1-dicyanothenyl)-1-phenyl-4,5-dihydro-1H-pryazole (DCNP) were prepared. The polymer films were investigated with in situ second-harmonic generation (SHG) measurement. The corona poling temperature was optimized by the temperature dependence of the in situ SHG signal intensity under the poling electric field applying. The temporal and temperature stability of the second-order properties of the poled polymer film were measured by the in situ SHG signal intensity probing. The second-order NLO coefficient chi ((2))(33) = 32.65 pm/V at lambda = 1064 nm was determined by using the Makel fringe method after poling under the optimal poling condition. The dispersion of the NLO coefficient of the guest-host polymer system was determined by the measured value of chi ((2))(33) at 1064 nm and the two-level model.

Electro-optical properties and temporal stability of the guest-host DCNP/PEK-c polymer thin film

The polyetherketone (PEK-c) guest-host system thin films in which the range of the weight percent of 3-(1,1-dicyanothenyl)-1-phenyl-4, 5- dihydro-1H-pryazole (DCNP) is from 20% to 50% were prepared. The predicted high value of electro-optical (EO) coefficient gamma(33) = 48.8 pm/V by using two-level model was obtained when the weight percent of DCNP in the polymer system is 40%, whereas EO coefficients are attenuated at higher chromophore loading then 40%. The temporal stability of the EO activity of the guest-host polymer was evaluated by probing the decay of the orientational order of the chromophores in the polymer system.

Photoluminescence studies of GaAs/AlGaAs single quantum well intermixed by Ga+ ion implantation

Ga(+)ion implantation followed by rapid thermal annealing (RTA) was used to enhance the interdiffusion in GaAs/AlGaAs single Quantum Wells(SQWs). The extent of intermixing was found to be dependent on the well depth, number of implanted ions and annealing time. A very fast interdiffusion process occurs at the initial annealing stage. After that, the enhanced diffusion coefficient goes back to the umimplanted value. We propose a two-step model to explain the diffusion process as a function of the annealing time : a fast diffusion process and a saturated diffusion process. The interdiffusion coefficient of the fast diffusion was found to be of well depth dependence and estimated to be in the range of 5.4x10(-16) similar to 1.5x10(-15)cm(2)s(-1). Copyright (C) 1996 Published by Elsevier Science Ltd

Temperature Distribution in Ridge Structure InGaN Laser Diodes and Its Influence on Device Characteristics

Time-dependent thermal simulation of ridge-geometry InGaN laser diodes is carried out with a two-dimensional model. A high temperature in the waveguide layer and a large temperature step between the regions under and outside the ridge are generated due to the poor thermal conductivity of the sapphire substrate and the large threshold current and voltage. The temperature step is thought to have a strong influence on the characteristics of the laser diodes. Time-resolved measurements of light-current curves,spectra, and the far-field pattern of the InGaN laser diodes under pulsed operation are performed. The results show that the thermal lensing effect improves the confinement of the higher order modes and leads to a lower threshold current and a higher slope efficiency of the device while the high temperature in the active layer results in a drastic decrease in the slope efficiency.

Study of turbulent heat transfer of aviation kerosene flows in a curved pipe at supercritical pressure

The now and heat transfer characteristics of China No. 3 aviation kerosene in a heated curved tube under supercritical pressure are numerically investigated by a finite volume method. A two-layer turbulence model, consisting of the RNG k-epsilon two-equation model and the Wolfstein one-equation model, is used for the simulation of turbulence. A 10-species kerosene surrogate model and the NIST Supertrapp software are applied to obtain the thermophysical and transport properties of the kerosene at various temperature under a supercritical pressure of 4 MPa. The large variation of thermophysical properties of the kerosene at the supercritical pressure make the flow and heat transfer more complicated, especially under the effects of buoyancy and centrifugal force. The centrifugal force enhances the heat transfer, but also increases the friction factors. The rise of the velocity caused by the variation of the density does not enhance the effects of the centrifugal force when the curvature ratios are less than 0.05. On the contrary, the variation of the density increases the effects of the buoyancy. (C) 2010 Elsevier Ltd. All rights reserved.

Short and long waves over a muddy seabed

The available experimental results have shown that in time-periodic motion the rheology of fluid mud displays complex viscoelastic behaviour. Based on the measured rheology of fluid mud from two field sites, we study the interaction of water waves and fluid mud by a two-layered model in which the water above is assumed to be inviscid and the mud below is viscoelastic. As the fluid-mud layer in shallow seas is usually much thinner than the water layer above, the sharp contrast of scales enables an approximate analytical theory for the interaction between fluid mud and small-amplitude waves with a narrow frequency band. It is shown that at the leading order and within a short distance of a few wavelengths, wave pressure from above forces mud motion below. Over a Much longer distance, waves are modified by the accumulative dissipation in mud. At the next order, infragravity waves owing to convective inertia (or radiation stresses) are affected indirectly by mud motion through the slow modulation of the short waves. Quantitative predictions are made for mud samples of several concentrations and from two different field sites.

Physical mechanism and numerical simulations of surface layer temperature inversion in tropical ocean

The one-dimensional Kraus-Turner mixed layer model improved by Liu is developed to consider the effect of salinity and the equations of temperature and salinity under the mixed layer. On this basis, the processes of growth and death of surface layer temperature inversion is numerically simulated under different environmental parameters. At the same time, the physical mechanism is preliminarily discussed combining the observations at the station of TOGA-COARE 0 degrees N, 156 degrees E. The results indicate that temperature inversion sensitively depends on the mixed layer depth, sea surface wind speed and solar shortwave radiation, etc., and appropriately meteorological and hydrological conditions often lead to the similarly periodical occurrence of this inversion phenomenon.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

太平洋西部边界及边界流对赤道暖水库形成和演化的影响

A model of equatorial ocean is used to study the roles of the Pacific western boundary and the Mindanao Current (MC) in the evolution of the equatorial warm pool. The model consists of the single baroclinic mode of a two-layer ocean, with the parameterization of the anomalous increment of the interface representing the SST difference from its long-term-space-mean. The ocean is driven by a wind path in the middle ocean with a real or an artificial geometry assigned at the western and eastern boundaries. In order to test the role of the MC, the western boundary current is introduced into the model by a boundary condition at a position, real and unreal, respectively. The model experiments show that the warm pool, which is insensitive to the longitudinal width of the wind band in middle ocean, results mainly from the accumulation o the eastly-drifted warm water in the equatorial western Pacific. It is the dominant factor for the formation of the warm pool that, at a very low latitude, the Papua New Guinea coast intersects the longitudinally lined Philippine Islands at an obtuse angle. In contrast, the western Atlantic boundary, which inclines poleward from the equator at some 135 degrees, could guide the warm water there moving to a higher latitude. On the other hand, the equatorial warm pool in the western equatorial Pacific is very sensitive to the assignment of th Mindanao Current at 7.5°N and displaces southward, with a stronger southern branch than the northern one. We attribute this asymmetry to the combined effect of the western boundary and the MC upon the equatorial warm away from the equator. A by-product of our solutions is the possible mechanism of the "secondary warm pool" in the eastern Pacific north of the equator. It is suggested that, mainly or partly, the "secondary warm pool" results from the cooperation of the southeast monsoon in eastern Pacific and the eastern boundary hindering the propagation of the Kelvin wave poleward alongshore.

中太平洋海区富钴结壳地球化学特征及成因机制

本文利用ICP-AES、EPMA、X-ray衍射等测试技术以及聚类分析和因子分析等多种数理统计方法，系统地对中太平洋海区富钴结壳的元素地球化学特征、矿物组成和微观组构进行了研究，并探讨了其成因机制，获得以下主要认识： 1中太平洋海区富钴结壳类型多样，均为水成成因，其矿物相主要由锰矿物相、铁矿物相和非金属矿物相组成；富钴结壳壳层发育多种原生和次生构造类型。 2磷酸盐化作用不仅强烈改变富钴结壳元素初始含量，而且造成富钴结壳某些元素间的相关性发生改变，这些敏感型元素对可用于指示富钴结壳是否发生磷酸盐化。在不同水深段内，未磷酸盐化型富钴结壳的主要元素随经向、纬向的变化趋势相似，表明其受水体化学障、表面生产力和物质来源等环境参数控制；而随水深的变化则具有区域一致性，表明水体化学具有区域成层性。 3未磷酸盐化富钴结壳稀土元素含量和轻重稀土分馏程度随水深发生规律变化，这种变化不仅与它们在海洋中的含量和行为有关，也与海洋背景颗粒的吸附有关；铈(Ce)在富钴结壳中基本上呈4价，且动力学因素控制了其富集过程，因此Ce异常不能用于指示富钴结壳形成环境的氧化程度。 4基于富钴结壳微层呈锯齿状且同一微层生长速率不同，提出了富钴结壳在各种基底表面生长以及后继发育过程受固液界面双电层控制的发育模式。在富钴结壳整个发育过程中，经历了从贫氧环境向富氧环境的转变，但微环境则呈富氧-低氧过程的交替。 关键词： 富钴结壳；中太平洋；元素地球化学；界面双电层

双模式相关技术研究与实现

文中研究的相关跟踪技术主要应用于飞航导弹的末制导。现阶段，激光制导技术、GPS制导技术、合成孔径雷达（SAR）制导技术在我国的实际应用还不成熟，传统的电视、红外制导技术仍然具有很强的生命力，因而，基于可见光的视频相关跟踪技术具有重要的研究价值。本文主要论述的两种相关算法是：多灰度点相关（MPC）、区域模板相关（RTC）。其中多点相关（MPC）算法的跟踪灵敏度高，定位精度好，硬件实现比较方便，实时性能好；区域模板相关算法（RTC），在图像的匹配过程中，不仅考虑了目标区域的灰度特征，而且兼顾了区域里多灰度层次的位置特征、面积特征，算法具有很好的鲁棒性。文中深入研究了两种相关跟踪算法，并针对它们在实际应用中的不足，提出了有效的改善措施。最后，本文对两种相关跟踪算法进行了初步融合，一是：通过粗匹配、精匹配过程来选取目标跟踪点；二是：提出了一种度量模板更新的能量准则函数。大量的仿真实验结果表明：改进后的两种相关跟踪技术可以较好地完成一些复杂背景下的目标跟踪任务，两种算法的有效结合又进一步提高了目标跟踪的稳定性能和可靠性能。本文研究的一些相关跟踪技术已经运用到实际工程项目中。

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

电离层垂直TEC映射函数的实验观测与统计特性

When used in the determining the total electron content (TEC), which may be the most important ionospheric parameter, the worldwide GPS observation brings a revolutionary change in the ionospheric science. There are three steps in the data processing to retrieve GPS TEC: (1) to estimate slant TEC from the measurements of GPS signals; (2) to map the slant TEC into vertical; and (3) to interpolate the vertical TEC into grid points. In this scientific dissertation we focus our attention on the second step, the mapping theory and method to convert slant TEC into vertical. This is conventionally done by multiplying on the slant TEC a mapping function which is usually determined by certain models of electron density profile. Study of the vertical TEC mapping function is of significance in GPS TEC measurement. This paper first reviews briefly the three steps in GPS TEC mapping process. Then we compare the vertical TEC mapping function which were respectively calculated from the electron density profiles of the ionospheric model and retrieved from the observation of worldwide GPS TEC. We also perform the statistical analysis on the observational mapping functions. The main works and results are as follows: 1. We calculated the vertical TEC mapping functions for both SLM and Chapman models, and discussed the modulation of the ionosphere height to the mapping functions. We use two simple models, single layer model (SLM) and Chapman models, of the ionospheric electron density profiles to calculate the vertical TEC mapping function. In the case of the SLM, we discuss the control of the ionospheric altitude, i.e., the layer height hipp, to the mapping function. We find that the mapping function decreases rapidly as hipp increases. For the Chapman model we study also the control mapping function by both ionospheric altitude indicated by the peak electron density height hmF2, and the scale height, H, which present the thickness of the ionosphere. It is also found that the mapping function decreases rapidly as hmF2 increases. and it also decreases as H increases. 2. Then we estimate the mapping functions from the GPS observations and compare them with those calculated from the electron density models. We first, proposed a new method to estimate the mapping functions from GPS TEC data. This method is then used to retrieve the observational mapping function from both the slant TEC (TECS) provided by International GPS Service (IGS）and vertical TEC provide by JPL Global Ionospheric Maps (GIMs). Then we compare the observational mapping function with those calculated from the electron density models, SLM and Chapman. We find that the values of the observational mapping functions are much smaller than that from the model mapping functions, when the zenith angle is large enough. We attribute this to the effect of the plasmasphere which is above about 1000 km. 3. We statistically analyze the observational mapping functions and reveal their climatological changes. Observational mapping functions during 1999-2007 are used in our statistics. The main results are as follows. (1) The observational mapping functions decrease obviously with the decrement of the solar activity which is represented by the F10.7 index; (2) In annual variations of the observational mapping functions, the semiannual component is found at low-latitudes, and the remarkable seasonal variations at mid- and high-latitudes. (3) The diurnal variation of the observational mapping functions is that they are large in daytime and small at night, they become extremely small in the early morning before sunrise. (4) The observational mapping functions change with latitudes that they are smaller at lower latitudes and larger at higher. All of the above variations of the observational mapping functions are explained by the existence of the plasmasphere, which changes more slowly with time and more rapidly with latitude than the ionosphere does . In summary, our study on the vertical TEC mapping function imply that the ionosphere height has a modulative effect on the mapping function. We first propose the concept of the 'observational mapping functions' , and provide a new method to calculate them. This is important in improving the TEC mapping. It may also possible to retrieving the plasmaspheric information from GPS observations.

青藏高原东部电性结构特征及地球动力学意义

As the largest and highest plateau on the Earth, the Tibetan Plateau has been a key location for understanding the processes of mountain building and plateau formation during India-Asia continent-continent collision. As the front-end of the collision, the geological structure of eastern Tibetan Plateau is very complex. It is ideal as a natural laboratory for investigating the formation and evolution of the Tibetan Plateau. Institute of Geophysics, Chinese Academy of Sciences (CAS) carried out MT survey from XiaZayii to Qingshuihe in the east part of the plateau in 1998. After error analysis and distortion analysis, the Non-linear Conjugate Gradient inversion(NLCG), Rapid Relaxation Inversin (RRI) and 2D OCCAM Inversion algorithms were used to invert the data. The three models obtained from 3 algorithms provided similar electrical structure and the NLCG model fit the observed data better than the other two models. According to the analysis of skin depth, the exploration depth of MT in Tibet is much more shallow than in stable continent. For example, the Schmucker depth at period 100s is less than 50km in Tibet, but more than 100km in Canadian Shield. There is a high conductivity layer at the depth of several kilometers beneath middle Qiangtang terrane, and almost 30 kilometers beneath northern Qiangtang terrane. The sensitivity analysis of the data predicates that the depth and resistivity of the crustal high conductivity layer are reliable. The MT results provide a high conductivity layer at 20~40km depth, where the seismic data show a low velocity zone. The experiments show that the rock will dehydrate and partially melt in the relative temperature and pressure. Fluids originated from dehydration and partial melting will seriously change rheological characteristics of rock. Therefore, This layer with low velocity and high conductivity layer in the crust is a weak layer. There is a low velocity path at the depth of 90-110 km beneath southeastern Tibetan Plateau and adjacent areas from seismology results. The analysis on the temperature and rheological property of the lithosphere show that the low velocity path is also weak. GPS measurements and the numerical simulation of the crust-mantle deformation show that the movement rate is different for different terranes. The regional strike derived from decomposition analysis for different frequency band and seismic anisotropy indicate that the crust and upper mantle move separately instead of as a whole. There are material flow in the eastern and southeastern Tibetan Plateau. Therefore, the faults, the crustal and upper mantle weak layers are three different boundaries for relatively movement. Those results support the "two layer wedge plates" geodynamic model on Tibetan formation and evolution.