Harmonic Millimeter Wave Generation and Frequency Up-Conversion Using Optical Injection Locking and Brillouin Selective Sideband Amplification

Harmonic millimeter wave (mm-wave) generation and frequency up-conversion are experimentally demonstrated using optical injection locking and Brillouin selective sideband amplification (BSSA) induced by stimulated Brillouin scattering in a 10-km single-mode fiber. By using this method, we successfully generate third-harmonic mm-wave at 27 GHz (f(LO) - 9 GHz) with single sideband (SSB) modulation and up-convert the 2GHz intermediate frequency signal into the mm-wave band with single mode modulation of the SSB modes. In addition, the mm-wave carrier obtains more than 23 dB power gain due to the BSSA. The transmission experiments show that the generated mm-wave and up-converted signals indicate strong immunity against the chromatic dispersion of the fibers.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

Simulating a typhoon storm surge in the East Sea of China using a coupled model

A coupled numerical model with a 2' x 2' resolution grid has been developed and used to simulate five typical typhoon storm surges (5612, 7413, 7910, 8114, and 9711) in the East Sea of China. Three main driving forces have been considered in this coupled model: wave radiation stress, combined wave-current bottom shear stress and wave-state-dependent surface wind stress. This model has then been compared with in situ measurements of the storm set-up. The effect of different driving force components on the total storm surge has also been investigated. This study has found that the coupled model with high resolution is capable of simulating the five typical typhoons better than the uncoupled models, and that the wave-dependent surface wind stress plays an important role in typhoon storm surge-wave coupling in this area and can increase the storm set-up by 1 m. The study of the five typhoon cases has shown that the general coupling effects could increase storm set-up by 20-32%. Thus, it is suggested that to predict typhoon storm surges in the East Sea of China, a storm surge-wave coupled model be adopted. (C) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

锦屏电站地下厂房卸荷岩体参数研究

It is a basic work to ascertain the parameters of rock mass for evaluation about stability of the engineering. Anisotropism、inhomogeneity and discontinuity characters of the rock mass arise from the existing of the structural plane. Subjected to water、weathering effect、off-loading, mechanical characters of the rock mass are greatly different from rock itself, Determining mechanical parameters of the rock mass becomes so difficult because of structure effect、dimension effect、rheological character, ‘Can’t give a proper parameter’ becomes one of big problems for theoretic analysis and numerical simulation. With the increment of project scale, appraising the project rock mass and ascertaining the parameters of rock mass becomes more and more important and strict. Consequently, researching the parameters of rock mass has important theoretical significance and actual meaning. The Jin-ping hydroelectric station is the first highest hyperbolic arch dam in the world under construction, the height of the dam is about 305m, it is the biggest hydroelectric station at lower reaches of Yalong river. The length of underground factory building is 204.52m, the total height of it is 68.83m, the maximum of span clearance is 28.90m. Large-scale excavation in the underground factory of Jin-ping hydroelectric station has brought many kinds of destructive phenomenon, such as relaxation、spilling, providing a precious chance for study of unloading parameter about rock mass. As we all know, Southwest is the most important hydroelectric power base in China, the construction of the hydroelectric station mostly concentrate at high mountain and gorge area, basically and importantly, we must be familiar with the physical and mechanical character of the rock mass to guarantee to exploit safely、efficiently、quickly, in other words, we must understand the strength and deformation character of the rock mass. Based on enough fieldwork of geological investigation, we study the parameter of unloading rock mass on condition that we obtain abundant information, which is not only important for the construction of Jin-ping hydroelectric station, but also for the construction of other big hydroelectric station similar with Jin-ping. This paper adopt geological analysis、test data analysis、experience analysis、theory research and Artificial Neural Networks （ANN） brainpower analysis to evaluate the mechanical parameter, the major production is as follows: (1)Through the excavation of upper 5-layer of the underground powerhouse and the statistical classification of the main joints fractures exposed, We believe that there are three sets of joints, the first group is lay fracture, the second group and the fourth group are steep fracture. These provide a strong foundation for the following calculation of and analysis; (2)According to the in-situ measurement about sound wave velocity、displacement and anchor stress, we analyses the effects of rock unloading effect，the results show a obvious time-related character and localization features of rock deformation. We determine the depth of excavation unloading of underground factory wall based on this. Determining the rock mass parameters according to the measurement about sound wave velocity with characters of low- disturbing、dynamic on the spot, the result can really reflect the original state, this chapter approximately the mechanical parameters about rock mass at each unloading area; (3)Based on Hoek-Brown experienced formula with geological strength index GSI and RMR method to evaluate the mechanical parameters of different degree weathering and unloading rock mass about underground factory, Both of evaluation result are more satisfied; (4)From the perspective of far-field stress, based on the stress field distribution ideas of two-crack at any load conditions proposed by Fazil Erdogan (1962)，using the strain energy density factor criterion (S criterion) proposed by Xue changming(1972)，we establish the corresponding relationship between far-field stress and crack tip stress field, derive the integrated intensity criterion formula under the conditions of pure tensile stress among two line coplanar intermittent jointed rock，and establish the corresponding intensity criterion for the exploratory attempt; (5)With artificial neural network, the paper focuses on the mechanical parameters of rock mass that we concerned about and the whole process of prediction of deformation parameters, discusses the prospect of applying in assessment about the parameters of rock mass，and rely on the catalog information of underground powerhouse of Jinping I Hydropower Station, identifying the rock mechanics parameters intellectually，discusses the sample selection, network design, values of basic parameters and error analysis comprehensively. There is a certain significance for us to set up a set of parameters evaluation system，which is in construction of large-scale hydropower among a group of marble mass.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

井间地震技术研究和应用

Cross well seismic technique is a new type of geophysical method, which observes the seismic wave of the geologic body by placing both the source and receiver in the wells. By applying this method, it averted the absorption to high-frequency component of seismic signal caused by low weathering layers, thus, an extremely high-resolution seismic signal can be acquired. And extremely fine image of cross well formations, structure, and reservoir can be achieved as well. An integrated research is conducted to the high-frequency S-wave and P-wave data and some other data to determine the small faults, small structure and resolving the issues concerning the thin bed and reservoir's connectivity, fluid distribution, steam injection and fracture. This method connects the high-resolution surface seismic, logging and reservoir engineering. In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a research is conducted on cross well seismic technology in general and its important issues in cross well seismic technology in particular. A technological series of integrated field acquisition, data processing and interpretation and its integrated application research were developed and this new method can be applied to oilfield development and optimizing oilfield development scheme. The contents and results in this paper are as listed follows: An overview was given on the status quo and development of the cross well seismic method and problems concerning the cross well seismic technology and the difference in cross well seismic technology between China and international levels; And an analysis and comparison are given on foreign-made field data acquisition systems for cross-well seismic and pointed out the pros and cons of the field systems manufactured by these two foreign companies and this is highly valuable to import foreign-made cross well seismic field acquisition system for China. After analyses were conducted to the geometry design and field data for the cross well seismic method, a common wave field time-depth curve equation was derived and three types of pipe waves were discovered for the first time. Then, a research was conducted on the mechanism for its generation. Based on the wave field separation theory for cross well seismic method, we believe that different type of wave fields in different gather domain has different attributes characteristics, multiple methods (for instance, F-K filtering and median filtering) were applied in eliminating and suppressing the cross well disturbances and successfully separated the upgoing and downgoing waves and a satisfactory result has been achieved. In the area of wave field numerical simulation for cross well seismic method, a analysis was conducted on conventional ray tracing method and its shortcomings and proposed a minimum travel time ray tracing method based on Feraiat theory in this paper. This method is not only has high-speed calculation, but also with no rays enter into "dead end" or "blinded spot" after numerous iterations and it is become more adequate for complex velocity model. This is first time that the travel time interpolation has been brought into consideration, a dynamic ray tracing method with shortest possible path has been developed for the first arrivals of any complex mediums, such as transmission, diffraction and refraction, etc and eliminated the limitation for only traveling from one node to another node and increases the calculation accuracy for minimum travel time and ray tracing path and derives solution and corresponding edge conditions to the fourth-order differential sonic wave equation. The final step is to calculate cross well seismic synthetics for given source and receivers from multiple geological bodies. Thus, real cross-well seismic wave field can be recognized through scientific means and provides important foundation to guide the cross well seismic field geometry designing. A velocity tomographic inversion of the least square conjugated gradient method was developed for cross well seismic velocity tomopgraphic inversion and a modification has been made to object function of the old high frequency ray tracing method and put forward a thin bed oriented model for finite frequency velocity tomographic inversion method. As the theory model and results demonstrates that the method is simple and effective and is very important in seismic ray tomographic imaging for the complex geological body. Based on the characteristics of the cross well seismic algorithm, a processing flow for cross well seismic data processing has been built and optimized and applied to the production, a good section of velocity tomopgrphic inversion and cross well reflection imaging has been acquired. The cross well seismic data is acquired from the depth domain and how to interprets the depth domain data and retrieve the attributes is a brand new subject. After research was conducted on synthetics and trace integration from depth domain for the cross well seismic data interpretation, first of all, a research was conducted on logging constraint wave impedance of cross well seismic data and initially set up cross well seismic data interpretation flows. After it applied and interpreted to the cross well seismic data and a good geological results has been achieved in velocity tomographic inversion and reflection depth imaging and a lot of difficult problems for oilfield development has been resolved. This powerful, new method is good for oilfield development scheme optimization and increasing EOR. Based on conventional reservoir geological model building from logging data, a new method is also discussed on constraining the accuracy of reservoir geological model by applying the high resolution cross well seismic data and it has applied to Fan 124 project and a good results has been achieved which it presents a bight future for the cross well seismic technology.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Experiments on interaction between current-induced vibration and scour of submarine pipelines on sandy bottom

In order to understand the dynamic behavior of submarine pipelines exposed to current and the mechanism of the interaction between current-induced vibration and scour of pipelines on a sandy bottom, an experimental investigation is conducted with a small scale model A test model which can be tested in the flume is set up by taking into account the typical working conditions of the pipelines and by applying the similarity theory. The interactions between the shape of the scour hole and the behavior of the pipeline as well as the flow patterns of the current are detailed, and the interaction mechanism outlined. The effect of vibration of the pipeline on the development of dynamic scour at different stages is found out. The proposed experimental method and test results provide an effective means for design of marine pipelines against scouring.

重力对异向聚集过程影响的探讨

自90年代初,分散体系聚集过程的微重力研究开始受到重视,并有一系列论文发表[1－4].Folkersma[3,4]等人在最近发表的文章中,公布了他们的最新研究结果.在探空火箭实现的微重力条件下,他们发现聚集速率要比地面实验快11倍多.为了解释这一现象,他们把重力的影响分为两部分,即沉降和自然对流.用密度匹配法,检验了在无沉降时的聚集过程.得到的结果是,聚集速率在无重力沉降影响时会比有沉降时快2.7倍.即便如此,微重力实验的聚集速率仍高出4.2倍.他们推断,这一巨大差异是由在地面上无法避免的自然对流造成的.

A study of the influence of free convection on coagulation process

It was assumed [1, 2] that gravity affects the coagulation process in two ways: free convection, which is hard to be avoided on the ground and sedimentation, which can be greatly reduced by the density-matching method. We present a ground-based experiment set-up to study the influence of convection on the perikinetic coagulation for aqueous polystyrene (PS) dispersions. The turbidity measurement was used to evaluate the relative coagulation rate and convection-driven flows in the solution were checked with a visual-magnification system. The pattern of flow field temperature profile in the sample cell is given. Our experiments show that there was no noticeable difference of coagulation rate observed no matter whether convection flows exist (with the flow speed up to 180 mu m/s) or not.

Characterising Of Internal Stresses In Duplex Coating By Fem

The internal stresses in a duplex coating involving a prequenched layer are believed to change if it is exposed to thermal loading. To characterise the internal stresses in such a duplex coating, a gradient model of finite element method is set up. The initial stress within the substrate developed in as quenching and the internal stresses due to the tempering of the prequenched layer ( TPQL) in such a duplex coating are calculated. The synthetical internal stresses in coating can be estimated by superposing uniform initial stresses developed during plating. The results indicate that the residual tensile stresses due to fabrication in coating will be decreased greatly, or even synthetical compressive internal stresses may arise in the coating.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

Distribution of Temperature in GlidArc Discharge Plasma

The GlidArc discharge is one of the main generation methods of non-equilibrium plasma near atmospheric pressures. In general, Gliding Arc discharge is driven by gas flow [1] in axial direction or by magnetic field in circumferential direction. [2] In this paper, a GlidArc discharge driven by rotating-gas-flow in circumferential direction is presented. The principle of the plasma generator is analyzed. The distribution of the temperature in axial direction is measured by a digital thermometer for three different gases. The experimental set-up of the GlidArc plasma is shown in Fig.1. It consists of a center electrode, an outside electrode, a power supply and a gas supply. The shortest distance between the electrodes is 2-3 mm. When a power supply with 10000 volts is attached to the electrodes, the arc will be ignited at the shortest distance. The small plasma column is rotated by the rotating gas flow in circumferential direction and then the rotating arc is driven towards the exit of the setup by the gas flow.