Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

High order asymptotic field for plane stress crack problem

This paper presents an exact analysis for high order asymptotic field of the plane stress crack problem. It has been shown that the second order asymptotic field is not an independent eigen field and should be matched with the elastic strain term of the first order asymptotic field. The second order stress field ahead of the crack tip is quite small compared with the first order stress field. The stress field ahead of crack tip is characterized by the HRR field. Hence the J integral can be used as a criterion for crack initiation.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Asymptotic solutions to a class of nonlinear oscillation systems

In this paper, we present an asymptotic method for the analysis of a class of strongly nonlinear oscillators, derive second-order approximate solutions to them expressed in terms of their amplitudes and phases, and obtain the equations governing the amplitudes and phases, by which the amplitudes of the corresponding limit cycles and their behaviour can be determined. As an example, we investigate the modified van der Pol oscillator and give the second-order approximate analytical solution of its limit cycle. The comparison with the numerical solutions shows that the two results agree well with each other.

The interaction of two pairs of solitary waves in a two-fluid system

This paper deals with the interaction of solitary waves in a two-fluid system which consistsof two superimposed incompressible inviscid fluids with a free surface and a horizontal rigidbottom. Under the assumption of shallow water wave, we first derive the basic equationssuitable for the model considered, a generalized form of the Boussinesq equations, then usingthe PLK method and the reductive perturbation method, obtain the second-order approximatesolution for the head-on collision between two pairs of interface and surface solitary waves,and give their maximum amplitudes during the collision and the nonuniform phase shiftsafter the collision which lead to the distortion of the wave profiles.

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.

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.

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

Displacive To Order-Disorder Two-Step Phase Transition Model For Para-Ferroelectric Transition

Two-step phase transition model, displacive to order-disorder, is proposed. The driving forces for these two transitions are fundamentally different. The displacive phase transition is one type of the structural phase transitions. We clearly define the structural phase transition as the symmetry broking of the unit cell and the electric dipole starts to form in the unit cell. Then the dipole-dipole interaction takes place as soon as the dipoles in unit cells are formed. We believe that the dipole-dipole interaction may cause an order-disorder phase transition following the displacive phase transition. Both structural and order-disorder phase transition can be first-order or second-order or in between. We found that the structural transition temperatures can be lower or equal or higher than the order-disorder transition temperature. The para-ferroelectric phase transition is the combination of the displacive and order-disorder phase transitions. It generates a variety of transition configurations along with confusions. In this paper, we discuss all these configurations using our displacive to order-disorder two-step phase transition model and clarified all the confusions.

Ejection of energetic photoelectrons in laser-induced autoionisation: effects of high-order ionisation processes

Resonant interaction of an autoionising state with a strong laser field is considered and effects of second-order ionisation processes are investigated. The authors show that these processes play a very important role in laser-induced autoionisation (LIA). They drastically affect the lowest-order peaks in the photoelectron spectrum. In addition to these peaks, high-order peaks due to ejection of energetic photoelectrons appear. For the laser intensities of current interest, second-order peaks are much stronger than the original ones, an important result that, they believe, can be observed experimentally. Moreover, `peak switching', a general feature of above-threshold ionisation, is also manifest in the electron spectrum of LIA.

Ghost imaging with thermal light by third-order correlation

Ghost imaging with classical incoherent light by third-order correlation is investigated. We discuss the similarities and the differences between ghost imaging by third-order correlation and by second-order correlation, and analyze the effect from each correlation part of the third-order correlation function on the imaging process. It is shown that the third-order correlated imaging includes richer correlated imaging effects than the second-order correlated one, while the imaging information originates mainly from the correlation of the intensity fluctuations between the test detector and each reference detector, as does ghost imaging by second-order correlation.

Corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media

The coupled differential recurrence equations for the corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media are established in terms of the perturbation method. All the corrections (including the longitudinal field corrections) to the paraxial approximation solutions are presented in the weak-guidance approximation. As a concrete application, the first-order longitudinal field correction and the second-order transverse field correction to the paraxial approximation of a Gaussian beam propagating in a transversely quadratic refractive index medium are analytically investigated. (C) 1999 Optical Society of America [S0740-3232(99)00310-5].

Second-harmonic generation in Ge20As25S55 glass irradiated by an electron beam

Second-harmonic generation was observed in Ge(20)AS(25)S(55) chalcogenide glass irradiated by an electron beam. The second-harmonic intensity increased with increasing electron-beam current and accelerating voltage. The second-harmonic generation in Ge20As25S55 glass was caused by the space-charge electrostatic field that was generated by irradiation of an electron beam. Second-order nonlinearity chi ((2)) as great as 0.8 pm/V was obtained. The results of measurements of thermally stimulated depolarization current indicated that the glass was poled in the thin layers of its surface (several micrometers) and that the nonlinearity was stable. (C) 2001 Optical Society of America.

Photoinduced stable second-harmonic generation in chalcogenide glasses

We report on photoinduced second-harmonic generation (SHG) in chalcogenide glasses. Fundamental and second-harmonic waves from a nanosecond pulsed Nd:YAG laser were used to induce second-order nonlinearity in chalcogenide glasses. The magnitude of SHG in 20Ge . 20As . 60S glass was 10(4) larger than that of tellurite glass with a composition of 15Nb(2)O(5) . 85TeO(2) (mol.%). Moreover, no apparent decay of photoinduced SHG in 20Ge . 20As . 60S glass was observed after optical poling at room temperature. We suggest that the large and stable value of X-(2) is due to the induced defect structures and large X-(3) of the chalcogenide glasses. (C) 2001 Optical Society of America

Lowered laser threshold in photonic crystals investigated by semiclassical theory

Using the response formula of the photonic crystals and the Bloch equations, the lasing threshold in arbitrary 2D photonic crystals was obtained by an investigation of steady-state laser behavior. The lasing threshold is expressed by the population inversion. It shows that the population inversion threshold is proportional to the second order of the group velocity, and to the relaxation coefficient. (c) 2004 Elsevier B.V. All rights reserved.