Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM: A comparative study of implementation strategies

Cardiovascular diseases are one of the major causes of long‐term morbidity and mortality in human beings. The nearly epidemic increase in prevalence of such diseases poses a serious threat to public health and calls for efficient methods of diagnosis and treatment. Non‐invasive diagnostic procedures such as MRI are often used in this context; however, these are limited in terms of spatial and temporal resolution and do not provide information on time‐dependent pressures and wall shear stresses—key quantities considered to be partially responsible for the formation and development of related pathologies. The present study is concerned with the numerical simulation of oscillatory flow through the abdominal aortic bifurcation. Computational fluid dynamics simulation of oscillatory flow in a branched geometry at high Reynolds numbers poses considerable challenges. The present study reports a detailed comparison of simulations performed with a finite volume and a finite element method, two approaches with significant differences in their discretization strategy, treatment of boundary conditions and other numerical aspects. Both solvers were parallelized, using loop parallelization of the BiCGStab linear solver for the finite volume and domain decomposition based on the Schur complement method for the finite element technique. The experience gained with these two approaches for the solution of flow in a bifurcation forms the focus of this study. Although similar results were obtained for both methods, the computation time required for convergence was found to be significantly smaller for the finite element approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Some exact solutions of the oscillatory motion of a generalized second grade fluid in an annular region of two cylinders

The velocity field and the associated shear stress corresponding to the longitudinal oscillatory flow of a generalized second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. Initially, the fluid and cylinders are at rest and at t = 0＋ both cylinders suddenly begin to oscillate along their common axis with simple harmonic motions having angular frequencies Ω1 and Ω2. The solutions that have been obtained are presented under integral and series forms in terms of the generalized G and R functions and satisfy the governing differential equation and all imposed initial and boundary conditions. The respective solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for the similar flow of ordinary second grade fluid and Newtonian fluid are also obtained as limiting cases of our general solutions. At the end, the effect of different parameters on the flow of ordinary second grade and generalized second grade fluid are investigated graphically by plotting velocity profiles.     

Implementation of physical boundary conditions into computational domain in modelling of oscillatory bottom boundary layers

This paper discusses the importance of realistic implementation of the physical boundary conditions into computational domain for the simulation of the oscillatory turbulent boundary layer flow over smooth and rough flat beds. A mathematical model composed of the Reynolds averaged Navier–Stokes equation, turbulent kinetic energy (k) and dissipation rate of the turbulent kinetic energy (ε) has been developed. Control‐volume approach is used to discretize the governing equations to facilitate the numerical solution. Non‐slip condition is imposed on the bottom surface, and irrotational main flow properties are applied to the upper boundary. The turbulent kinetic energy is zero at the bottom, whereas the dissipation rate is approaching to a constant value, which is proportional to the kinematic viscosity times the second derivative of the turbulent kinetic energy. The output of the model is compared with the available experimental studies conducted in oscillatory tunnels and wave flume. It is observed that the irrotational flow assumption at the upper boundary is not realistic in case of water tunnels. Therefore, new upper boundary conditions are proposed for oscillatory tunnels. The data of wave flume show good agreement with the proposed numerical model. Additionally, several factors such as grid aspect ratio, staggered grid arrangement, time‐marching scheme and convergence criteria that are important to obtain a robust, realistic and stable code are discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

Fermi-edge singularities in the photoluminescence spectrum of modulation-doped GaAs v-groove quantum wires

We report the observation of strong Fermi-edge singularities in the photoluminescence spectrum of strongly-confined, modulation-doped GaAs v-groove quantum wires. The behaviour of the singularity has been investigated at high excitation intensity, and both lattice and electrical heating. The latter produces a strong reduction of the singularity due to Fermi surface smearing, whereas, increased photoexcitation produces complex electron–hole correlation effects.     

OSCILLATING CRITERIA FOR CERTAIN EVEN ORDER NEUTRAL DIFFERENTIAL EQUATIONS WITH DEVIATING ARGUMENTS

By using the integral averaging technique, oscillating criteria for certain even order neutral differential equation with deviating arguments are established. These results extend some known oscillation criteria due to Xu and Xia.     

Isolated complex singularities and their CR links

In this paper, we discuss some recent studies on the complex structure of an isolated normal singularity by using the information from its link. We also give some open problems to be further pursued. Dedicated to Professor Sheng GONG on the occasion of his 75th birthday     

Characterization of isolated homogeneous hypersurface singularities in C~4 Dedicated to Professor Sheng GONG on the occasion of his 75th birthday

Let V be a hypersurface with an isolated singularity at the origin in Cn 1. It is a natural question to ask when V is defined by weighted homogeneous polynomial or homogeneous polynomial up to biholomorphic change of coordinates. In 1971, a beautiful theorem of Saito gives a necessary and sufficient condition for V to be defined by a weighted homogeneous polynomial. For a two-dimensional isolated hypersurface singularity V, Xu and Yau found a coordinate free characterization for V to be defined by a homogeneous polynomial. Recently Lin and Yau gave necessary and sufficient conditions for a 3-dimensional isolated hypersurface singularity with geometric genus bigger than zero to be defined by a homogeneous polynomial. The purpose of this paper is to prove that Lin-Yau's theorem remains true for singularities with geometric genus equal to zero.     

三圆盘扭振系统主共振的理论与实验研究

本文应用非线性振动的平均法，分析了一个具有立方非线性并受简谐激励作用的三圆盘扭振系统的第二阶主共振，求得了稳态响应的分岔方程，并进行了奇异性分析，理论结果与实验结果相符·