Two-dimensional simulation of high-order laterally-coupled GaAs-AlGaAs DFB laser diodes

Laterally-coupled distributed feedback (LC-DFB) laser diodes made without an epitaxial re-growth process have the advantage of a simple fabrication process. In this paper, two-dimensional optical field distribution of the fundamental quasi TE (transverse electric) mode is calculated by means of a semivectorial finite-difference method (SV-FDM). The dependence of the effective coupling coefficient (kappa(eff)) on the dutycycle of first-, second- and third-order LC-DFB LDs is investigated using modified coupled wave equations.

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.

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.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

High-Resolution Finite Compact Difference Schemes for Hyperbolic Conservation Laws

A finite compact (FC) difference scheme requiring only bi-diagonal matrix inversion is proposed by using the known high-resolution flux. Introducing TVD or ENO limiters in the numerical flux, several high-resolution FC-schemes of hyperbolic conservation law are developed, including the FC-TVD, third-order FC-ENO and fifth-order FC-ENO schemes. Boundary conditions formulated need only one unknown variable for third-order FC-ENO scheme and two unknown variables for fifth-order FC-ENO scheme. Numerical test results of the proposed FC-scheme were compared with traditional TVD, ENO and WENO schemes to demonstrate its high-order accuracy and high-resolution.

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

Quasi-Dammann Grating With Proportional Intensity Array Spots

A quasi-Dammann grating is proposed to generate array spots with proportional-intensity orders in the far field. To describe the performance of the grating, the uniformities of the array spots are redefined. A two-dimensional even-sampling encode scheme is adopted to design the quasi-Dammann grating. Numerical solutions of the binary-phase quasi-Dammann grating with proportional-intensity orders are given. The experimental results with a third-order quasi-Dammann grating, which has an intensity proportion of 3:2:1 from zero order to second order, are presented. (C) 2008 Optical Society of America

A micromechanical model of influence of particle fracture and particle cluster on mechanical properties of metal matrix composites

A general analytical model for a composite with an isotropic matrix and two populations of spherical inclusions is proposed. The method is based on the second order moment of stress for evaluating the homogenised effective stress in the matrix and on the secant moduli concept for the plastic deformation. With Webull's statistical law for the strength of SiCp particles, the model can quantitatively predict the influence of particle fracture on the mechanical properties of PMMCs. Application of the proposed model to the particle cluster shows that the particle cluster has neglected influence on the strain and stress curves of the composite. (C) 1998 Elsevier Science B.V.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

较高阶非线性抽运-探测过程中的空间效应

在飞秒抽运一探测光谱技术中，空间分辨的探测光信号反映了在不同空间位置的材料的非线性效应。当抽运光强度增大时，探测光信号中会出现明显的高阶特别是五阶非线性效应。利用劈裂算子方法直接解决了一维非线性传播方程的问题。在数值模拟中，研究了在不同抽运强度和位置下的抽运一探测过程中的五阶非线性效应。在足够高的抽运场下，探测信号出现清晰的振荡，显示了三阶和五阶非线性效应之间的干涉。当空间位置离抽运场中心足够远时，五阶比三阶非线性效应的衰减快得多，对其物理机制和趋势进行了定性的讨论。

XPM-induced pulse compression in birefringent dispersion decreasing fiber

The pulse compression induced by cross-phase modulation in birefringent dispersion decreasing fiber is discussed theoretically by solving the coupled Schrodinger equations which include the contribution of the high-order non-linear effects, and third-order dispersion. In particular, it is found that a high quality compressed signal pulse can be obtained by a pump pulse of low intense through the technique. The dependence of optimum compression on the non-linear factor N, time delay tau(d) and the dispersive ratio f is also discussed in detail. (c) 2005 Elsevier B.V. All rights reserved.

Ultrafast non-linear optical properties of Ge20As25Se55 chalcogenide films

The non-resonant third-order non-linear optical properties of amorphous Ge20As25Se55 films were studied experimentally by the method of the femtosecond optical heterodyne detection of optical Kerr effect. The real and imaginary parts of complex third-order optical non-linearity could be effectively separated and their values and signs could be also determined, which were 6.6 x 10(-12) and -2.4 x 10(-12) esu, respectively. Amorphous Ge20As25Se55 films showed a very fast response in the range of 200 fs under ultrafast excitation. The ultrafast response and large third-order non-linearity are attributed to the ultrafast distortion of the electron orbitals surrounding the average positions of the nucleus of Ge, As and Se atoms. The high third-order susceptibility and a fast response time of amorphous Ge20As25Se55 films makes it a promising material for application in advanced techniques especially in optical switching. (c) 2005 Elsevier B.V. All rights reserved.

Femtosecond optical Kerr effect study of Ge(10)AS(40)S(30)Se(20) film

The results of the femtosecond optical heterodyne detection of optical Kerr effect at 805 nm with the 80 fs ultrafast pulses in amorphous Ge10As40S30Se20 film is reported in this paper. The film shows an optical non-linear response of: 200 fs under ultrafast 80 fs-pulse excitation and the values of real and imaginary parts of non-linear susceptibility chi((3)) were 9.0 X 10(-12) and -4.0 X 10(-12) esu, respectively. The large third-order non-linearity and ultrafast response are attributed to the ultrafast distortion of the electron orbits surrounding the average positions of the nucleus of Ge, As, S and Se atoms. This Ge10As40S30Se20 chalcogenide glass would be expected as a promising material for optical switching technique. (c) 2005 Elsevier Ltd. All rights reserved.