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.

Dielectric response of spherically anisotropic graded piezoelectric composites

A graded piezoelectric composite consisting of a spherically anisotropic graded piezoelectric inclusion imbedded in an infinite nonpiezoelectric matrix, with the physical properties of the graded spherical inclusion having a power-law profile with respect to the radial variable r, is studied theoretically. Under an external uniform electric field, the electric displacement field and the elastic stress tensor field of this spherically anisotropic graded piezoelectric composite are derived exactly by means of displacement separation technique, based on the governing equations in the dilute limit. A piezoelectric response mechanism, in which the effective piezoelectric response vanishes along the z direction (or x,y directions), is revealed in this kind of graded piezoelectric composites. Furthermore, it is found that the effective dielectric constant decreases (or increases) with the volume fraction p of the inclusions if the exponent parameter k of the grading profile is larger (or smaller) than a critical value. (C) 2007 American Institute of Physics.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.

Recovery of the solitons using a lattice Boltzmann model

We formulate a lattice Boltzmann model which simulates Korteweg-de Vries equation by using a method of higher moments of lattice Boltzmann equation. Using a series of lattice Boltzmann equations in different time scales and the conservation law in time scale to, we obtain equilibrium distribution function. The numerical examples show that the method can be used to simulate soliton.

Anisotropic Deformation Of Polystyrene Particles By Mev Au Ion Irradiation

MeV An irradiation leads to a shape change of polystyrene (PS) and SiO2 particles from spherical to ellipsoidal, with an aspect ratio that can be precisely controlled by the ion fluence. Sub-micrometer PS and SiO2 particles were deposited on copper substrates and irradiated with Au ions at 230 K, using an ion energy and fluence ranging from 2 to 10 MeV and 1 x 10(14) ions/cm(2) to 1 x 10(15) ions/cm(2). The mechanisms of anisotropic deformation of PS and SiO2 particles are different because of their distinct physical and chemical properties. At the start of irradiation, the volume of PS particles decrease, then the aspect ratio increases with fluence, whereas for SiO2 particles the volume remains constant. (C) 2008 Elsevier B.V. All rights reserved.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

The magnetostatic relations including fluctuation fields

The influences of the fluctuation fields are important in many astrophysical environments as shown by the observations, and can not be neglected. On the basis of the first-order smoothing approximation, in the present paper, we demonstrate the magnetostatic equations for both the cases of the conventional turbulence aud the random waves, and discuss the consistent conditions of the equations. In the static problem, the fluctuation Lorentz force(▽×δB)×δB influences the large-scale configurations of magnetic field. To study this influence in detail is quite necessary for the explanations of the observation features, especially for the astrophysical environments where the magnetic fields, including the fluctuation fields, are the dominant factors in the equilibrium of momentum and energy.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<In the light of this reason there exists a great disparity of the grid Reynolds numbers RDxi between different cells in Navier-Stokes (NS) equations computations for high Reynolds number flows, an idea of solving the conservation equations for discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

参变过程中的群速延迟对晶体匹配长度的影响

基于能量守恒和三波耦合波方程, 建立了超短脉冲在参变过程中二次谐波产生时的I类和II类相位匹配条件、基波与谐波之间的群速延迟时间、以及群速失配对晶体长度限制的理论基础。以负单轴非线性光学晶体CsLiB6O10为例, 分析和数值计算了超短脉冲宽度为100 fs时, 谐波的群速匹配长度随基波波长变化的规律。研究结果表明在I类相位匹配条件下, 基波波长为642 nm时, 群速延迟最小, 相应的群速匹配晶体长度最长为19.1 mm;在II类相位匹配条件下, 基波波长为767 nm, 群速延迟最小, 群速匹配长度最

Experimental study of nonlinear switching characteristics of conventional 2×2 fused tapered couplers

The nonlinear switching characteristics of fused fiber directional couplers were studied experimentally. By using femtosecond laser pulses with pulse width of 100 fs and wavelength of about 1550 nm from a system of Ti:sapphire laser and optical parametric amplifier (OPA), the nonlinear switching properties of a null coupler and a 100% coupler were measured. The experimental results were coincident with the simulations based on nonlinear propagation equations in fiber by using super-mode theory. Nonlinear loss in fiber was also measured to get the injected power at the coupler. After deducting the nonlinear loss and input efficiency, the nonlinear switching critical peak powers for a 100% and a null fused couplers were calculated to be 9410 and 9440 W, respectively. The nonlinear loss parameter P_(N) in an expression of α_(NL)=αP/P_(N) was obtained to be P_(N)=0.23 W.

平面交叉型微透镜阵列的制作及成像特性研究

介绍了采用光刻离子交换工艺制作平面交叉型微透镜阵列的方法。利用积分形式的光线方程式讨论了平面交叉型微透镜的近轴光学特性，研究了微透镜的光线轨迹方程式和一些重要的近轴成像特性，利用ABCD定理得到了平面交叉型微透镜像距、焦距、像高、横向放大率和主平面位置的数学表达式，焦距的理论计算结果和实验数据吻合得很好。