六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非特殊振荡的直接原因是在于非均一的波群的群速度，因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式，以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例，证明了文中所构造的格式具有方法简便，物理含义清楚，精度高，网格基架点小，和捕捉激波能力较强的优点。

Elastic Waves in a Hybrid Multilayered Piezoelectric Plate

An analytical-numerical method is presented for analyzing dispersion and characteristic surface of waves in a hybrid multilayered piezoelectric plate. In this method, the multilayered piezoelectric plate is divided into a number of layered elements with three-nodal-lines in the wall thickness, the coupling between the elastic field and the electric field is considered in each element. The associated frequency dispersion equation is developed and the phase velocity and slowness, as well as the group velocity and slowness are established in terms of the Rayleigh quotient. Six characteristic wave surfaces are introduced to visualize the effects of anisotropy and piezoelectricity on wave propagation. Examples provide a full understanding for the complex phenomena of elastic waves in hybrid multilayered piezoelectric media.

Behaviour analysis of numerical solutions and simulation of coherent structures in compressible mixing layers

For understanding the correctness of simulations the behaviour of numerical solutions is analysed, Tn order to improve the accuracy of solutions three methods are presented. The method with GVC (group velocity control) is used to simulate coherent structures in compressible mixing layers. The effect of initial conditions for the mixing layer with convective Mach number 0.8 on coherent structures is discussed. For the given initial conditions two types of coherent structures in the mixing layer are obtained.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Space-frequency coupling, conical waves, and small-scale filamentation in water

Numerical simulations of fs laser propagation in water have been made to explain the small-scale filaments in water we have observed by a nonlinear fluorescence technique. Some analytical descriptions combined with numerical simulations show that a space-frequency coupling mainly from the interplay among self-phase modulation, dispersion and phase mismatching will reshape the laser beam into a conical wave which plays a major role of energy redistribution and can prevent laser beam from self-guiding over a long distance. An effective group velocity dispersion is introduced to explain the pulse broadening and compression in the filamentation. (c) 2005 American Institute of Physics.

Few-cycle spatiotemporal soliton wave excited by filamentation of a femtosecond laser pulse in materials with anomalous dispersion

The nonlinear dynamics of 1.6-mu m fs laser pulses propagating in fused silica is investigated by employing a full-order dispersion model. Different from the x-wave generation in normally dispersive media, a few-cycle spatiotemporally compressed soliton wave is generated with the contrary contributions of anomalous group velocity dispersion (GVD) and self-phase-modulation. However, at the tailing edge of the pulse forms a shock wave which generates separate and strong supercontinuum peaked at 670 nm. It is also the origin of conical emission formed both in time and frequency domain with the contribution of normal GVD at visible light.

Nonlinear copropagation of two optical pulses of different frequencies in photonic crystal fibers

A theoretical investigation of the nonlinear copropagation of two optical pulses of different frequencies in a photonic crystal fiber is presented. Different phenomena are observed depending on whether the wavelength of the signal pulse is located in the normal or the anomalous dispersion region. In particular, it is found that the phenomenon of pulse trapping occurs when the signal wavelength is located in the normal dispersion region while the pump wavelength is located in the anomalous dispersion region. The signal pulse suffers cross-phase modulation by the Raman shifted soliton pulse and it is trapped and copropagates with the Raman soliton pulse along the fiber. As the input peak power of the pump pulse is increased, the red-shift of the Raman soliton is considerably enhanced with the simultaneous further blue-shift of the trapped pulse to satisfy the condition of group velocity matching.

不同功率下的飞秒激光在水中形成的脉冲压缩和超连续谱

本文研究了飞秒激光脉冲在水中的传输情况.通过改变不同的激光输入功率进行模拟，我们发现从输入功率略高于到远远高于发生自聚焦的临界功率,分别是群速度色散和多光子电离多光子吸收阻止了自聚焦导致的脉冲塌陷，当多光子电离和多光子吸收主导传输时，脉冲能被压缩到几个光学周期.在频域，多光子电离能引起很强的蓝移,而多光子吸收能对这种蓝移起到抑制作用。

超短激光脉冲在不同色散参量光子晶体光纤中传输的数值模拟

为了使得数值模拟更为精确, 采用广义非线性薛定谔方程(GNSE)描述超短激光脉冲在光子晶体光纤中的传输演化过程, 并利用二阶分步傅里叶方法通过求解方程, 数值计算了相同脉宽和能量的超短脉冲在不同色散参量的光子晶体光纤中非线性传输和超连续谱的产生。比较了超短脉冲在光纤不同色散区传输时, 高阶色散和非线性效应对超连续谱的产生以及对脉冲波形演化的影响。结果表明, 相对于超短脉冲中心波长位于光子晶体光纤的正常和反常色散区, 可以相应获得短波波段和长波波段的超连续谱输出, 当超短脉冲中心波长位于零色散波长点时, 通

Conical emission by femtosecond pulses with different spectral bandwidths

Conical emission (CE) has been investigated experimentally by laser pulses with different pulse durations and spectral bandwidths. The results show that the overall CE curve will shift as the varying of spectral bandwidth of pump laser pulse. But for pump laser pulses which have same spectral bandwidth but different pulse duration, the CE angles will be same at the spectral region close to the pump wavelength while will be different at the spectral region far away from the pump wavelength. We have also fitted the measured CE angles with X-wave model. The calculated curves and the measured CE curves match reasonably well. The best fits indicate that the group velocity of the filament pulse may be greatly controlled by controlling the spectral bandwidth of pump laser pulse. (C) 2008 Elsevier B.V. All rights reserved.

Formation of X-waves at fundamental and harmonics by infrared femtosecond pulse filamentation in air

We experimentally observe the formation of X-waves at fundamental, third harmonic, and fifth harmonic wavelengths by infrared (central wavelength at similar to 1500 nm) femtosecond laser pulse filamentation in air. By fitting the angularly resolved spectra of the fundamental and harmonic waves using X-wave relations, we confirm that all the X-waves have nearly the same group velocity, indicating that they are locked in space and time during their propagation in filament.