Turbulent passive scalar field of a small prandtl number

The statistical-mechanics theory of the passive scalar field convected by turbulence, developed in an earlier paper [Phys. Fluids 28, 1299 (1985)], is extended to the case of a small molecular Prandtl number. The set of governing integral equations is solved by the equation-error method. The resultant scalar-variance spectrum for the inertial range is F(k)~x−5/3/[1+1.21x1.67(1+0.353x2.32)], where x is the wavenumber scaled by Corrsin's dissipation wavenumber. This result reduces to the − (5)/(3) law in the inertial-convective range. It also approximately reduces to the − (17)/(3) law in the inertial-diffusive range, but the proportionality constant differs from Batchelor's by a factor of 3.6.

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<numbers RDxi between different cells in Navier-Stokes (NS) equations computations for high Reynolds number flows, an idea of solving the conservation equations for discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-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.

Smaller Deborah number inducing more serrated plastic flow of metallic glass

Spherical nano-indentations of Cu46Zr54 bulk metallic glass (BMG) model systems were performed using molecular dynamics (MD) computer simulations, focusing specifically on the physical origin of serrated plastic flow. The results demonstrate that there is a direct correlation between macroscopic flow serration and underlying irreversible rearrangement of atoms, which is strongly dependent on the loading (strain) rate and the temperature. The serrated plastic flow is, therefore, determined by the magnitude of such irreversible rearrangement that is inhomogeneous temporally. A dimensionless Deborah number is introduced to characterize the effects of strain rate and temperature on serrations. Our simulations are shown to compare favorably with the available experimental observations.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

A universal Biot number determining the susceptibility of ceramics to quenching

A universal Biot number of ceramics, which not only determines the susceptibility of the ceramics to quenching but also indicates the duration that the ceramics fail during thermal shock, is theoretically obtained. The present analysis shows that the thermal shock failure of the ceramics with a Biot number greater than this universal value is a very rapid process that just occurs in the initial regime of the heat conduction of the ceramics. This universal Biot number provides a guide to the selection of the ceramics applying to the thermostructural engineering including thermal shock.

Number of phase levels of a Talbot array illuminator

The number of phase levels of a Talbot array illuminator is an important factor in the estimation of practical fabrication complexity and cost. We show that the number it) of phase levels of a Talbot array illuminator has a simple relationship to the prime number. When there is an alternative pi -phase modulation in the output array, the relations are similar. (C) 2001 Optical Society of America OCIS codes: 070.6760, 050.1950, 050.1980.

A New theory for evaluating the number density of inclusions in films

A new formulation derived from thermal characters of inclusions and host films for estimating laser induced damage threshold has been deduced. This formulation is applicable for dielectric films when they are irradiated by laser beam with pulse width longer than tens picoseconds. This formulation can interpret the relationship between pulse-width and damage threshold energy density of laser pulse obtained experimentally. Using this formulation, we can analyze which kind of inclusion is the most harmful inclusion. Combining it with fractal distribution of inclusions, we have obtained an equation which describes relationship between number density of inclusions and damage probability. Using this equation, according to damage probability and corresponding laser energy density, we can evaluate the number density and distribution in size dimension of the most harmful inclusions. (c) 2005 Elsevier B.V. All rights reserved.