Wavelength conversion based on degenerate-four-wave-mixing with continuous-wave pumping in silicon nanowire waveguide

We present a comprehensive numerical study on the all-optical wavelength conversion based on the degenerate four-wave-mixing with continuous-wave pumping in the silicon nanowire waveguide. It is well known that the conversion efficiency and the 3-dB bandwidth can be greatly affected by the phase-matching condition. Through proper design of the waveguide cross-section, its dispersion property can be adjusted to satisfy the phase-matching condition and therefore effective wavelength conversion can be achieved in a large wavelength range. Generally, the group velocity dispersion plays a dominant role in the wavelength conversion. However, the fourth-order dispersion takes an important effect on the wavelength conversion when the group velocity dispersion is near the zero-point. Furthermore, the conversion efficiency and the 3-dB bandwidth can also be affected by the interactive length and the initial pump power. Through the numerical simulation, the optimal values for the interactive length and the initial pump power, which are functions of the propagation loss, are obtained to realize the maximum conversion efficiency. (C) 2008 Elsevier B.V. All rights reserved.

Investigation of spectral bandwidth of optical parametric amplification

The spectral bandwidth of three-wave-mixing optical parametric amplification has been investigated. A general mathematical model for evaluating the spectral bandwidth of optical parametric amplification is developed with parametric bandwidth and gain bandwidth via three-wave noncollinear interactions. The spectral bandwidth is determined by expanding the wave-vector mismatch in a Taylor series and retaining terms through second order. The model takes into account the effects of crystal length, noncollinear angle, group velocity, group-velocity dispersion and gain coefficient. The relation between parametric bandwidth and gain bandwidth is clearly defined. The model is applied to a BBO OPA, a LBO OPA and a CLBO OPA.

Zero-broadening SBS slow light propagation in an optical fiber using two broadband pump beams

A new method of tailoring stimulated Brillouin scattering (SBS) gain spectrum for slow light propagation is proposed by use of two Gaussian-shaped broadband pump beams with different powers and spectral widths. The central frequency interval between the two pump beams are carefully set to be two inherent Brillouin frequency shift, ensuring that the gain spectrum of one pump has the same central frequency with the loss spectrum of the other one. Different gain profiles are obtained and analyzed. Among them a special gain profile is found that ensures a zero-broadening of the signal pulse independent of the Brillouin gain. This is owing to the compensation between the positive gain-dependent broadening and the negative GVD (group velocity dispersion) dependent broadening. The relationship of two pump beams is also found for constructing such a gain profile. It provides us a new idea of managing the broadening of SBS-based slow pulse by artificially constructing and optimizing the profile of gain spectrum. (c) 2008 Optical Society of America.

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.

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.

Newly designed multilayer thin film mirror for dispersion compensation in Ti:sapphire femtosecond lasers

There are two different effects to generate group delay dispersion by multilayer thin film mirrors: chirper effect and Gires-Tournois effect. Both effects are employed to introduce desired dispersion in the designed mirror. Thus the designed mirror provides large dispersion throughout broad waveband. Such mirror can be used for dispersion compensation in Ti:sapphire femtosecond lasers. Most group delay dispersion of a 5-mm Ti:sapphire crystal can be compensated perfectly with only four bounces of the designed mirror.

Slow wave effect of 2-D photonic crystal coupled cavity array

We designed a two-dimensional coupled photonic crystal resonator array with hexagonal lattice. The calculation by plane-wave-expansion method shows that the dispersion curve of coupled cavity modes in the bandgap are much flattened in all directions in the reciprocal space. We simulated the transmission spectra of transverse electric (TE) wave along the Gamma K direction. Compared with the PC single cavity structure, the transmission ratio of the coupled cavity array increases more than three orders of magnitude, while the group velocity decreases to below 1/10, reaching 0.007c. The slow wave effect has potential application in the field of miniaturized tunable optical delay components and low-threshold photonic crystal lasers.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

8阶群速度控制格式及其应用

构造了8阶精度的群速度控制型差分格式.根据激波捕捉能力及对小尺度分辨能力对格式系数进行了优化.采用该格式对可压缩均匀各向同性湍流进行了直接数值模拟,所计算的最大湍流马赫数达到0.95.相同计算条件下的结果与他人结果吻合较好,说明了格式在捕捉激波的同时对湍流小尺度有较好的分辨能力.

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

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

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

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.