Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Numerical analyses for the material bifurcation in plane sheet under tension

The various patterns (shear banding, surface wrinkling and necking) of material bifurcation in plane sheet under tension are investigated in this paper by means of a numerical method. It is found that numerical analysis can provide better ground for searching for the lowest critical loads. The inhomogeneity caused by void damage and the nonuniformity in the stress distribution across sheet thickness are proved to have detrimental effects on the material bifurcation. Nevertheless, material stability can be promoted by any means of depressing void damage or alleviating stress, even locally across the thickness. Besides, the peculiar behaviour of material bifurcation under slight biaxiality state is demonstrated. Copyright (C) 1996 Elsevier Science Ltd

Micromechanics analyses of interface crack

A new mechanics model based on Peierls concept is presented in this paper, which can clearly characterize the intrinsic features near a tip of an interfacial crack. The stress and displacement fields are calculated under general combined tensile and shear loadings. The near tip stress fields show some oscillatory behaviors but without any singularity and the crack faces open completely without any overlapping when remote tensile loading is comparable with remote shear loading. A fracture criterion for predicting interface toughness has been also proposed, which takes into account for the shielding effects of emitted dislocations. The theoretical toughness curve gives excellent prediction, as compared with the existing experiment data.

Analyses of the mechanical properties and microstructure of bamboo epoxy composites

Bamboo reinforced epoxy possesses reasonably good properties to waarrant its use as a structural material, and is fabricated by utilizing bamboo, an abundant material resource, in the technology of fibre composites. Literature on bamboo-plastics composites is rare. This work is an experimental study of unidirectional bamboo-epoxy laminates of varying laminae number, in which tensile, compressive, flexural and interlaminar shear properties are evaluated. Further, the disposition of bamboo fibre, the parenchymatous tissue, and the resin matrix under different loading conditions are examined. Our results show that the specific strength and specific modulus of bamboo-epoxy laminates are adequate, the former being 3 to 4 times that of mild steel. Its mechanical properties are generally comparable to those of ordinary glass-fibre composites. The fracture behaviour of bamboo-epoxy under different loading conditions were observed using both acoustic emission techniques and scanning electron microscopy. The fracture mode varied with load, the fracture mechanism being similar to glass and carbon reinforced composites. Microstructural analyses revealed that natural bamboo is eligibly a fibre composite in itself; its inclusion in a plastic matrix will help solve the problems of cracking due to desiccation and bioerosion caused by insect pests. Furthermore, the thickness and shape of the composite can be tailored during fabrication to meet specific requirements, thereby enabling a wide spectrum of applications.

Spectral-line profile of turbulent gas

phases should be specified when the particle Reynolds number is higher than

Elastic-plastic finite element analyses of short cracks

Stress and strain distributions and crack opening displacement characteristics of short cracks have been studied in single edge notch bend and centre cracked panel specimens using elastic–plastic finite element analyses incorporating both a non strain hardening and a power law hardening behaviour. J contour integral solutions to describe stress strain conditions at crack tips for short cracks differ from those for long cracks. The analyses show that (i) short cracks can propagate at stress levels lower than those required for long cracks and (ii) a two-parameter description of crack tip fields is necessary for crack propagation.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Space-resolved spectral diagnosis of line-shaped laser-produced magnesium plasmas

Space-resolved spectra of line-shaped laser-produced magnesium plasmas in the normal direction of the target have been obtained using a pinhole crystal spectrograph. These spectra are treated by a spectrum analyzing code for obtaining the true spectra and fine structures of overlapped lines. The spatial distributions of electron temperature and density along the normal direction of the target surface have been obtained with different spectral diagnostic techniques. Especially, the electron density plateaus beyond the critical surface in line-shaped magnesium plasmas have been obtained with a fitting technique applied to the Stark-broadened Ly-alpha wings of hydrogenic ions. The difference of plasma parameters between those obtained by different diagnostic techniques is discussed. Other phenomena, such as plasma satellites, population inversion, etc., which are observed in magnesium plasmas, are also presented.

Spectral effects for an ultrashort pulse train propagating in a two-level atom medium

We investigate the spectra of a femtosecond pulse train propagating in a resonant two-level atom (TLA) medium. it is found that higher spectral components can be produced even for a 2 pi femtosecond pulse train. Furthermore, the spectral effects depend crucially on both the relative shift phi and the delay time tau between the successive pulses of the femtosecond pulse train.

Constrained high-order statistical blind deconvolution of spectral data

A constrained high-order statistical algorithm is proposed to blindly deconvolute the measured spectral data and estimate the response function of the instruments simultaneously. In this algorithm, no prior-knowledge is necessary except a proper length of the unit-impulse response. This length can be easily set to be the width of the narrowest spectral line by observing the measured data. The feasibility of this method has been demonstrated experimentally by the measured Raman and absorption spectral data.

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.

Conical Emission Patterns by Femtosecond Pulses with Different Spectral Bandwidths

Different conical emission (CE) patterns are obtained experimentally at various incident powers and beam sizes of pump laser pulses with pulse durations of 7 fs, 44 fs and 100 fs. The results show that it is the incident power but not the incident power density that determines a certain CE pattern. In addition, the critical powers for similar CE patterns are nearly the same for the laser pulses with the same spectral bandwidth. Furthermore, as far as a certain CE pattern is concerned, the wider the spectral bandwidth of pump laser pulse is, the higher the critical power is. This will hopefully provide new insights for the generation of CE pattern in optical medium.

Spectral evolution of angularly resolved third-order harmonic generation by infrared femtosecond laser-pulse filamentation in air

We experimentally investigate the evolution of an angularly resolved spectrum of third harmonic generated by infrared femtosecond laser pulse filamentation in air. We show that at low pump intensity, phase matching between the fundamental and third-harmonic waves dominates the nonlinear optical effect and induces a ring structure of the third-harmonic beam, whereas at high pump intensity, the dispersion properties of air begin to affect the angular spectrum, leading to the formation of a nonlinear X wave at third harmonic.