Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid

A solution methodology has been developed for incompressible flow in general curvilinear co‐ordinates. Two staggered grids are used to discretize the physical domain. The first grid is a MAC quadrilateral mesh with pressure arranged at the centre and the Cartesian velocity components located at the middle of the sides of the mesh. The second grid is so displaced that its corners correspond to the centre of the first grid. In the second grid the pressure is placed at the corner of the first grid. The discretized mass and momentum conservation equations are derived on a control volume. The two pressure grid functions are coupled explicitly through the boundary conditions and implicitly through the velocity of the field. The introduction of these two grid functions avoids an averaging of pressure and velocity components when calculating terms that are generated in general curvilinear co‐ordinates. The SIMPLE calculation procedure is extended to the present curvilinear co‐ordinates with double grids. Application of the methodology is illustrated by calculation of well‐known external and internal problems: viscous flow over a circular cylinder, with Reynolds numbers ranging from 10 to 40, and lid‐driven flow in a cavity with inclined walls are examined. The numerical results are in close agreement with experimental results and other numerical data. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

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.     

Quantum-Chemical Evaluation of the Protolytic Properties of Thiophenols

The enthalpies and free energies of proton affinity (PA) were calculated by the PM3 method for 27 thiophenoxide anions. For thiophenols, linear correlation has been found between the PAs and the values of pK _a measured by different authors in aqueous acetone, aqueous ethanol, and methanol media. The dependences found permit one to predict a priori the protolytic properties of thiophenols. The differential effects of aqueous ethanol media with different ethanol contents on thiophenol acidity were evaluated based on quantum-chemical data.Original Russian Text Copyright © 2004 by A. N. Pankratov and A. V. Shalabai__________Translated from Zhurnal Strukturnoi Khimii, Vol. 45, No. 5, pp. 800–806, September–October, 2004.     

Stress and heat flux stabilization for viscous compressible medium equations with a nonmonotone state function

We study a system of quasilinear equations describing one-dimensional flow of a viscous compressible heat-conducting medium with a nonmonotone state function and mass force. The large-time behavior of solutions is considered for arbitrarily large initial data. In spite of possible nonuniqueness and discontinuity of the stationary solution, we prove L²-stabilization for the stress and heat flux as t → ∞ along with corresponding global energy estimates for them. The new method of proof utilizes a combination of energy type equalities for the stress and heat flux. Consequently, H¹-stabilization of the velocity and temperature along with global estimates for their derivatives are valid as well.     

Radiation Transfer in a Binary, Stochastic Medium with Strongly Fluctuating Parameters

The propagation of radiation in a three-dimensional, inhomogeneous, stochastic, scattering medium and in an equivalent homogeneous medium has been investigated by the small-angle iteration method. The statistical structure of the scattering parameters of a three-dimensional medium has been described within the framework of the two-dimensional Poisson process and the binary Markov process. The analytical data obtained point to the fact that the three-dimensional inhomogeneity of a stochastic medium significantly influences the radiation transfer in it.     

Green chemical synthesis of fluorinated 1,3,5-triaryl-s-triazines in aqueous medium under microwaves as potential antifungal agents

Reasonable pure fluorinated s-triazines were synthesized in quantitative yield (96-99%) in 2-3 min in aqueous medium under microwaves, by reaction of fluorinated anilines and aqueous formaldehyde. All synthesized compounds have been screened in vitro for their antifungal activity against Rhizoctonia solani, Fusarium oxysporum, and Collectotrichum capsici.     

Pattern forming precipitation in gels due to coupling of chemical reactions with diffusion

Precipitate-forming chemical reactions have been studied in chemically cross-linked poly(vinyl alcohol) gel medium. One of the reactive components was incorporated into the gel, the other was allowed to diffuse into it. Depending on the experimental conditions the reaction-diffusion process often results in patterns of different type. Experiments performed in tubes and in thin layers were carried out in order to investigate the effects of various factors (cross-linking density, swelling degree as well as the concentrations of the outer and inner electrolytes) on the morphologies of the precipitate patterns. It was found that precipitation occurs not only in the Liesegang bands, but also between bands. Beside Liesegang-type structures, tree-like patterns have been observed, showing a characteristic periodicity in the density profile obtained by digitalized image analyses.