A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

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.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Laser intensity effects in photoemission and photoelectron spectra produced in laser-induced autoionisation

Strong laser-field-induced autoionisation in the presence of both photoionising and radiative decay of the autoionising state (AS) is investigated, focusing on the laser intensity dependence of the photoemission and photoelectron spectra. In contrast to previous predictions, power broadening and increasing reduction of the doublet peak heights with field strength are found in the photoemission spectrum. Similar effects leading to considerable suppression and even complete disappearance of the lowest-order peaks in the photoelectron spectrum, together with peak switching, are also demonstrated, which are closely related to above-threshold ionisation. In addition, it is suggested that the total number of energetic photoelectrons may serve as an alternative to measuring the atomic parameters of the AS. All these effects are attributed to the presence of the strong `probe': laser-induced decay of the AS.

Diffraction integral formulas of the pulsed wave field in the temporal domain

To resolve the diffraction problems of the pulsed wave field directly in the temporal domain, we extend the Rayleigh diffraction integrals to the temporal domain and then discuss the approximation condition of this diffraction formula. (C) 1997 Optical Society of America.

The intensity distributions of collected signals in coherent anti-Stokes Raman scattering microscopy

Coherent anti-Stokes Raman scattering (CARS) microscopy with the combining of confocal and CARS techniques is a remarkable alternative for imaging chemical or biological specimens that neither fluoresce nor tolerate labeling. The CARS is a nonlinear optical process, the imaging properties of CARS microscopy will be very different from the conventional confocal microscopy. In this paper, we calculated the propagation of CARS signals by using the wave equation in medium and the slowly varying envelope approximation (SVEA), and find that the intensity angular distributions vary considerably with the different experimental configurations and the different specimen shapes. So the conventional description of microscopy (e.g.. the point spread function) will fail to descript the imaging properties of CARS microscopy. (c) 2004 Elsevier B.V. All rights reserved.

Generating monoenergetic heavy-ion bunches with laser-induced electrostatic shocks

A method for efficient laser acceleration of heavy ions by electrostatic shock is investigated using particle-in-cell (PIC) simulation and analytical modeling. When a small number of heavy ions are mixed with light ions, the heavy ions can be accelerated to the same velocity as the light ions so that they gain much higher energy because of their large mass. Accordingly, a sandwich target design with a thin compound ion layer between two light-ion layers and a micro-structured target design are proposed for obtaining monoenergetic heavy-ion beams.

Laser intensity dependence of high-order harmonic generation from aligned CO2 molecules

We investigate experimentally the high-order harmonic generation from aligned CO2 molecules and demonstrate that the modulation inversion of the harmonic yield with respect to molecular alignment can be altered dramatically by fine-tuning the intensity of the driving laser pulse for harmonic generation. The results can be modeled by employing the strong field approximation including a ground state depletion factor. The laser intensity is thus proved to be a parameter that can control the high-harmonic emission from aligned molecules.

Intensity-dependent asymmetric photoionization in few-cycle laser pulses

The asymmetric photoionization of atoms irradiated by intense, few-cycle laser pulses is studied numerically. The results show that the pulse intensity affects the asymmetric photoionization in three aspects. First, at higher intensities, the asymmetry becomes distinctive for few-cycle pulses of longer durations. Second, as the laser intensity increases, the maximal asymmetry first decreases then increases after it has reached a minimal value. Last, the value of the carrier-envelope phase corresponding to the maximal asymmetry varies with the pulse intensity. This study reveals that the increasing of pulse intensity is helpful for observing the asymmetric photoionization.

Short-pulse laser absorption via J x B heating in ultrahigh intensity laser plasma interaction

An analytical fluid model for JxB heating during the normal incidence by a short ultraintense linearly polarized laser on a solid-density plasma is proposed. The steepening of an originally smooth electron density profile as the electrons are pushed inward by the laser is included self-consistently. It is shown that the JxB heating includes two distinct coupling processes depending on the initial laser and plasma conditions: for a moderate intensity (a <= 1), the ponderomotive force of the laser light can drive a large plasma wave at the point n(e)=4 gamma(0)n(c) resonantly. When this plasma wave is damped, the energy is transferred to the plasma. At higher intensity, the electron density is steepened to a high level by the time-independent ponderomotive force, n(e)> 4 gamma(0)n(c), so that no 2 omega resonance will occur, but the longitudinal component of the oscillating ponderomotive field can lead to an absorption mechanism similar to "vacuum heating." (c) 2006 American Institute of Physics.