International Journal for Numerical Methods in Fluids

Prasanth, T. K.; Mittal, S.

doi: 10.1002/fld.1771pmid: N/A

The effect of the blockage on vortex‐induced vibrations of a circular cylinder of low non‐dimensional mass (m*=10) in the laminar flow regime is investigated in detail. A stabilized space–time finite element formulation is utilized to solve the incompressible flow equations in primitive variables form in two dimensions. The transverse response of the cylinder is found to be hysteretic at both ends of synchronization/lock‐in region for 5% blockage. However, for the 1% blockage hysteresis occurs only at the higher Re end of synchronization/lock‐in region. Computations are carried out at other blockages to understand its effect on the hysteretic behavior. The hysteresis loop at the lower Re end of the synchronization decreases with decrease in blockage and is completely eliminated for blockage of 2.5% and less. On the other hand, hysteresis persists for all values of blockage at the higher Re end of synchronization/lock‐in. Although the peak transverse oscillation amplitude is found to be same for all blockage (∼0.6D), the peak value of the aerodynamic coefficients vary significantly with blockage. The r.m.s. values show lesser variation with blockage. The effect of streamwise extent of computational domain on hysteretic behavior is also studied. The phase between the lift force and transverse displacement shows a jump of almost 180° at, approximately, the middle of the synchronization region. This jump is not hysteretic and is independent of blockage. Copyright © 2008 John Wiley & Sons, Ltd.

Valencia, Alvaro; Ledermann, Darren; Rivera, Rodrigo; Bravo, Eduardo; Galvez, Marcelo

doi: 10.1002/fld.1786pmid: N/A

Hemodynamics plays an important role in the progression and rupture of cerebral aneurysms. The current work describes the blood flow dynamics and fluid–structure interaction in seven patient‐specific models of bifurcating cerebral aneurysms located in the anterior and posterior circulation regions of the circle of Willis. The models were obtained from 3D rotational angiography image data, and blood flow dynamics and fluid–structure interaction were studied under physiologically representative waveform of inflow. The arterial wall was assumed to be elastic, isotropic and homogeneous. The flow was assumed to be laminar, non‐Newtonian and incompressible. In one case, the effects of different model suppositions and boundary conditions were reported in detail. The fully coupled fluid and structure models were solved with the finite elements package ADINA. The vortex structure, pressure, wall shear stress (WSS), effective stress and displacement of the aneurysm wall showed large variations, depending on the morphology of the artery, aneurysm size and position. The time‐averaged WSS, effective stress and displacement at the aneurysm fundus vary between 0.17 and 4.86 Pa, 4.35 and 170.2 kPa and 0.16 and 0.74 mm, respectively, for the seven patient‐specific models of bifurcating cerebral aneurysms. Copyright © 2008 John Wiley & Sons, Ltd.

Liu, G. R.; Xu, George X.

doi: 10.1002/fld.1788pmid: N/A

A novel gradient smoothing method (GSM) based on irregular cells and strong form of governing equations is presented for fluid dynamics problems with arbitrary geometries. Upon the analyses about the compactness and the positivity of coefficients of influence of their stencils for approximating a derivative, four favorable schemes (II, VI, VII and VIII) with second‐order accuracy are selected among the total eight proposed discretization schemes. These four schemes are successively verified and carefully examined in solving Poisson's equations, subjected to changes in the number of nodes, the shapes of cells and the irregularity of triangular cells, respectively. Numerical results imply us that all the four schemes give very good results: Schemes VI and VIII produce a slightly better accuracy than the other two schemes on irregular cells, but at a higher cost in computation. Schemes VII and VIII that consistently rely on gradient smoothing operations are more accurate than Schemes II and VI in which directional correction is imposed. It is interestingly found that GSM is insensitive to the irregularity of meshes, indicating the robustness of the presented GSM. Among the four schemes of GSM, Scheme VII outperforms the other three schemes, for its outstanding overall performance in terms of numerical accuracy, stability and efficiency. Finally, GSM solutions with Scheme VII to some benchmarked compressible flows including inviscid flow over NACA0012 airfoil, laminar flow over flat plate and turbulent flow over an RAE2822 airfoil are presented, respectively. Copyright © 2008 John Wiley & Sons, Ltd.

Iskandarani, M.

doi: 10.1002/fld.1791pmid: N/A

The thin aspect ratio of oceanic basins is simultaneously a complication to contend with when developing ocean models and an opportunity to simplify the equations of motion. Here we discuss these two aspects of this geometric feature in the context of hydrostatic and non‐hydrostatic ocean models. A simple analysis shows that the horizontal viscous operator in the hydrostatic primitive equations plays a central role in the specification of boundary conditions on the lateral vertical surfaces bounding the domain. The asymptotic analysis shows that for very thin aspect ratios the leading‐order flow cannot be closed unless additional terms in the equations are considered, namely either the horizontal viscous forces or the non‐hydrostatic pressure forces. In either case, narrow boundary layers must be resolved in order to close the circulation properly. The computational cost increases substantially when non‐hydrostatic effects are taken into account. Copyright © 2008 John Wiley & Sons, Ltd.

Fernández Oro, J. M.; Argüelles Díaz, K. M.; Santolaria Morros, C.; Cobo Hedilla, A. F.; Lemaille, M.

doi: 10.1002/fld.1793pmid: N/A

This paper describes the numerical modelling of the three‐phase flow that is established when pouring molten glass onto a tin bath (float) under a non‐reactive atmosphere. Discharged from the furnace, the glass is spouted over a lip onto the bath, where it floats and spreads out in the form of a long ribbon, leading to the flat glass forming process. Numerically, the physical modelling must consider a multiphase problem with three immiscible phases: glass, tin and air. The simulation of this complex free‐surface flow involves the solution of the Navier–Stokes set of equations for all the phases simultaneously, using a volume‐of‐fluid formulation that introduces a marker function convected by the flow to identify each phase. The evolution of the interphases is tracked over time with the implementation of a continuous surface force algorithm. A general purpose, well‐tested commercial code, FLUENT, is employed for the computations. Firstly, a two‐dimensional model considering the symmetry plane of the float is developed to fit accurate physical and numerical parameters. The high complexity of the interphases in addition to great differences between the physical properties of the phases has required extensive tests to ensure the consistency and accuracy of the solutions. Afterwards, a complete three‐dimensional model is built to simulate the pouring process in a real geometry and predict the behaviour of the industrial facility when the operating conditions are modified. Typical flow phenomena inside the tin bath, such as tin currents, wet back flow or even the equilibrium thickness, are also obtained as a promising result of the numerical modelling. Copyright © 2008 John Wiley & Sons, Ltd.

van Loon, R.; van de Vosse, F. N.

doi: 10.1002/fld.1891pmid: N/A