Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

Reflection and transmission properties of double submerged rectangular blocks

An analytic method is used to study the reflection and transmission coefficients of the double submerged rectangular blocks (DSRBs) in oblique waves.. The scattering potentials are obtained by means of the eigenfunction expansion method, and expressions for the reflection and transmission coefficients are determined. The boundary element method is employed to verify the correctness of the present analytical method. The DSRBs have better performance than the single submerged rectangular block (SSRB) in certain cases. The reflection and transmission properties of the DSRBs are investigated for some specific cases, and the influences of the geometric parameters are also presented.

Theoretical design or single-polarization single-mode microstructured polymer optical fibres

Using the full-vector plane-wave expansion method, a kind of PMMA-based polarization-maintaining microstructured optical fibre (PM-mPOF) is theoretically studied. Dependence of the cutoff wavelengths of the two orthogonal polarization states (polarized along the two principal axes of PM-mPOF) on the structure parameters of the fibre is investigated in detail. A single-polarization single-mode (SPSM) PM-mPOF working in the visible region is designed and optimized with the result of the maximum SPSM bandwidth of 140 nm.

Joint statistical distribution of two-point sea surface elevations in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, a joint statistical distribution of two-point sea surface elevations is derived by using the characteristic function expansion method. It is found that the joint distribution depends on five parameters. These five parameters can all be determined by the water depth, the relative position of two points and the wave-number spectrum of ocean waves. As an illustrative example, for fully developed wind-generated sea, the parameters that appeared in the joint distribution are calculated for various wind speeds, water depths and relative positions of two points by using the Donelan and Pierson spectrum and the nonlinear effects of sea waves on the joint distribution are studied. (C) 2003 Elsevier B.V. All rights reserved.

Statistical distribution of depth-integrated local horizontal momentum for second-order random ocean waves in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth- integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all determined by the water depth and the wave-number spectrum of ocean waves. As an illustrative example, a fully developed wind-generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.

Statistical distribution of wave-surface elevation for second-order random directional ocean waves in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, a statistical distribution of the wave-surface elevation is derived by using the characteristic function expansion method. It is found that the distribution, after normalization of the wave-surface elevation, depends only on two parameters. One parameter describes the small mean bias of the surface produced by the second-order wave-wave interactions. Another one is approximately proportional to the skewness of the distribution. Both of these two parameters can be determined by the water depth and the wave-number spectrum of ocean waves. As an illustrative example, we consider a fully developed wind-generated sea and the parameters are calculated for various wind speeds and water depths by using Donelan and Pierson spectrum. It is also found that, for deep water, the dimensionless distribution reduces to the third-order Gram-Charlier series obtained by Longuet-Higgins [J. Fluid Mech. 17 (1963) 459]. The newly proposed distribution is compared with the data of Bitner [Appl. Ocean Res. 2 (1980) 63], Gaussian distribution and the fourth-order Gram-Charlier series, and found our distribution gives a more reasonable fit to the data. (C) 2002 Elsevier Science B.V. All rights reserved.

Probability distribution of random wave-current forces

Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave-wave and wave-current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan-Pierson-Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. (c) 2006 Elsevier Ltd. All rights reserved.

Effective thermal conductivity of graded nonlinear composites with heat contact resistance

The perturbation expansion method is used to find the effective thermal conductivity of graded nonlinear composites having thermal contact resistance on the inclusion surface. As an example, we have studied the graded composites with cylindrical inclusions immersed in a homogeneous matrix. The thermal conductivity of the cylindrical inclusion is assumed to have a power-law profile of the radial distance r measured from its origin. For weakly nonlinear constitutive relations between the heat flow density q and the temperature field T, namely, q = -mu del T - chi vertical bar del T vertical bar(2) del T, in both the inclusion and the matrix regions, we have derived the temperature distributions using the perturbation expansion method. A nonlinear effective medium approximation of graded composites is proposed to estimate the effective linear and nonlinear thermal conductivities. by considering the temperature singularity on the inclusion surface due to the heat contact resistance. (c) 2006 Elsevier B.V. All rights reserved.

基于波扩散方法的工业机器人路径规划

为避免繁琐的机器人示教过程,提出一种离线的基于波扩散方法的工业机器人路径规划算法．首先对机器人的工作空间离散化,针对工作空间中的障碍点和自由点进行二值标记;然后用波扩散方法对自由点进一步标记,并进行了路径搜索;最后,对波扩散法与深度优先算法路径搜索进行了比较．将该算法用于6-自由度工业机器人的仿真实验,得到了满意的效果．

基于MCX314as的运动控制器设计

文章介绍了用C8051F120单片机和MCX314as组成运动控制器的方案。该方案应用存储器扩展的方法为运动控制芯片MCX314as和液晶显示模块MT22G06B的内部寄存器分配地址空间,从而实现C8051F120单片机对运动控制芯片和液晶模块的控制作用。文章还给出了MCX314和MT22G06B的驱动程序架构。

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

An eigenfunction expansion-variational method in prediction of the transverse thermal conductivity of fiber reinforced composites considering interfacial characteristics

An eigenfunction expansion-variational method based on a unit cell is developed to deal with the steady-state heat conduction problem of doubly-periodic fiber reinforced composites with interfacial thermal contact resistance or coating. The numerical results show a rapid convergence of the present method. The present solution provides a unified first-order approximation formula of the effective thermal conductivity for different interfacial characteristics and fiber distributions. A comparison with the present high-order results, available experimental data and micromechanical estimations demonstrates that the first-order approximation formula is a good engineering closed-form formula. An engineering equivalent parameter reflecting the overall influence of the thermal conductivities of the matrix and fibers and the interfacial characteristic on the effective thermal conductivity, is found. The equivalent parameter can greatly simplify the complicated relation of the effective thermal conductivity to the internal structure of a composite. (c) 2010 Elsevier Ltd. All rights reserved.

Variational method for adiabatic compression of plasma with poloidal and toroidal rotation

A variational method is developed for adiabatic compression of plasma with both poloidal and toroidal rotation.

An application of adjoint variational method to the capillary instability in liquid

In this paper, applying the direct variational approach of first-order approximation to the capillary instability problem for the eases of rotating liquid column, toroid and films on both sides of cylinder, we have obtained the necessary and sufficient conditions for motion stability of the "cylindrical coreliquid-liquid-cylindrical shell" systems. The results obtained before are found to be special cases of the present investigation. At the same time, we have explained physical essence of rotating instability and settled a few disputes in previous investigations.

Research on the influence of fabrication parameters on the performance of buried ion-exchanged glass waveguides using the variational method

The variational method is proposed to analyze the influence of the fabrication parameters on the performance of buried K+-Na+ ion-exchanged Er3+-Yb3+ ions co-doped glass waveguide. The unknown parameters of the Hermite-Gaussian functions as the trial field distribution are determined based on the scalar variational principle. It is demonstrated that the results calculated in this paper agree with those measured in the experiment. The mode dimensions, the effective refractive index, and the overlap factor as the functions of the fabrication parameters are investigated. These results of the variational analysis are useful for the design and optimization of Er3+-Yb3+ ions co-doped waveguides.