Collisional Relaxation of Luminescence Anisotropy in Rarefied Gases in the Condensed Phase

The orientational relaxation of optically induced anisotropy in rarefied gases and at a damped rotation has been investigated. It has been found that the anisotropy relaxation in rarefied gases is described by a reduced kinetic equation depending only on free rotation integrals. The behavior of the integral anisotropy of luminescence for free symmetric and asymmetric top molecules has been elucidated. The law of luminescence depolarization has been obtained for asymmetric top molecules in the Gordon J-diffusion model. It represents the sum of two Stern–Volmer-type dependences, whose relative contribution is determined by the orientation of the dipole moments of transitions with absorption and emission of light in the molecular coordinate system and by the principal moments of inertia of the molecular top. It has been established that in the limit of a strongly damped rotation, kinetic equations of the general form reduce to equations of rotational diffusion. A number of modified diffusion equations correctly describing the contribution of inertial effects to the orientational relaxation of anisotropy have been obtained.     

Eigenvalues of the antiplane-shear crack problem for a power-law material

This paper discusses the problem of finding the eigenvalue spectrum in determining the stress and strain fields at the tip of an antiplane-shear crack in a power-law material. It is shown that the perturbation method provides an analytical dependence of the eigenvalue on the material nonlinearity parameter and the eigenvalue of the linear problem. Thus, it is possible to find the entire spectrum of eigenvalues and not only the eigenvalue of the Hutchinson-Rice-Rosengren problem. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 173–180, January–February, 2008.     

Numerical simulation on an interaction of a vortex street with an elliptical leading edge using an unstructured grid

An algorithm for a time accurate incompressible Navier–Stokes solver on an unstructured grid is presented. The algorithm uses a second order, three‐point, backward difference formula for the physical time marching. For each time step, a divergence free flow field is obtained based on an artificial compressibility method. An implicit method with a local time step is used to accelerate the convergence for the pseudotime iteration. To validate the code, an unsteady laminar flow over a circular cylinder at a Reynolds number of 200 is calculated. The results are compared with available experimental and numerical data and good agreements are achieved. Using the developed unsteady code, an interaction of a Karman vortex street with an elliptical leading edge is simulated. The incident Karman vortex street is generated by a circular cylinder located upstream. A clustering to the path of the vortices is achieved easily due to flexibility of an unstructured grid. Details of the interaction mechanism are analysed by investigating evolutions of vortices. Characteristics of the interactions are compared for large‐ and small‐scale vortex streets. Different patterns of the interaction are observed for those two vortex streets and the observation is in agreement with experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Calculation of impulsively started incompressible viscous flows

The paper's focus is the calculation of unsteady incompressible 2D flows past airfoils. In the framework of the primitive variable Navier–Stokes equations, the initial and boundary conditions must be assigned so as to be compatible, to assure the correct prediction of the flow evolution. This requirement, typical of all incompressible flows, viscous or inviscid, is often violated when modelling the flow past immersed bodies impulsively started from rest. Its fulfillment can however be restored by means of a procedure enforcing compatibility, consisting in a pre‐processing of the initial velocity field, here described in detail. Numerical solutions for an impulsively started multiple airfoil have been obtained using a finite element incremental projection method. The spatial discretization chosen for the velocity and pressure are of different order to satisfy the inf–sup condition and obtain a smooth pressure field. Results are provided to illustrate the effect of employing or not the compatibility procedure, and are found in good agreement with those obtained with a non‐primitive variable solver. In addition, we introduce a post‐processing procedure to evaluate an alternative pressure field which is found to be more accurate than the one resulting from the projection method. This is achieved by considering an appropriate ‘unsplit’ version of the momentum equation, where the velocity solution of the projection method is substituted. Copyright © 2004 John Wiley & Sons, Ltd.     

Passive control of the flow around a square cylinder using porous media

The passive control of bluff body flows using porous media is investigated by means of the penalization method. This method is used to create intermediate porous media between solid obstacles and the fluid in order to modify the boundary layer behaviour. The study covers a wide range of two‐dimensional flows from low transitional flow to fully established turbulence by direct numerical simulation of incompressible Navier–Stokes equations. A parametric study is performed to illustrate the effect of the porous layer permeability and thickness on the passive control. The numerical results reveal the ability of porous media to both regularize the flow and to reduce the drag forces up to 30%. Copyright © 2004 John Wiley & Sons, Ltd.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

延迟微分方程单支方法的非线性稳定性

本文讨论延迟微分方程单支方法的非线性稳定性 .对于 Kα,β,γ类非线性延迟微分方程 ,我们证明带有线性插值的 G( c,p,q) -代数稳定的单支方法当 c≤ 1时是 GR( p/2 ,q/2 ) -稳定及弱 GAR( p/2 ,q/2 ) -稳定的 ,当 c<1时是 GAR( p/2 ,q/2 ) -稳定的 .最后的数值试验表明了上述结论的正确性 .     

The Calogero–Bogoyavlenskii–Schiff Equation in 2+1 Dimensions

We use the classical and nonclassical methods to obtain symmetry reductions and exact solutions of the (2+1)-dimensional integrable Calogero–Bogoyavlenskii–Schiff equation. Although this (2+1)-dimensional equation arises in a nonlocal form, it can be written as a system of differential equations and, in potential form, as a fourth-order partial differential equation. The classical and nonclassical methods yield some exact solutions of the (2+1)-dimensional equation that involve several arbitrary functions and hence exhibit a rich variety of qualitative behavior.