Numerical computations of supersonic inlet flow

Air‐breathing propulsion systems for high‐speed space travel applications are studied. Ramjets and scramjets have been identified as potential candidates. The flow inlets of such systems are modelled with a simulation that can predict all complex inlet flow features, including shock due to forebody, multiple shock reflections, normal shock, shock–boundary layer interaction and associated separation for two‐dimensional and axisymmetric inlets. Computed values are in good agreement with experimental data. Copyright © 2001 John Wiley & Sons, Ltd.     

Investigation of two‐dimensional channel flow with a partially compliant wall using finite volume–finite difference approach

This paper presents a numerical simulation of steady two‐dimensional channel flow with a partially compliant wall. Navier–Stokes equation is solved using an unstructured finite volume method (FVM). The deformation of the compliant wall is determined by solving a membrane equation using finite difference method (FDM). The membrane equation and Navier–Stokes equation are coupled iteratively to determine the shape of the membrane and the flow field. A spring analogy smoothing technique is applied to the deformed mesh to ensure good mesh quality throughout the computing procedure. Numerical results obtained in the present simulation match well with that in the literature. Copyright © 2005 John Wiley & Sons, Ltd.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems

Embedded boundary methods for CFD (computational fluid dynamics) simplify a number of issues. These range from meshing the fluid domain, to designing and implementing Eulerian‐based algorithms for fluid–structure applications featuring large structural motions and/or deformations. Unfortunately, embedded boundary methods also complicate other issues such as the treatment of the wall boundary conditions in general, and fluid–structure transmission conditions in particular. This paper focuses on this aspect of the problem in the context of compressible flows, the finite volume method for the fluid, and the finite element method for the structure. First, it presents a numerical method for treating simultaneously the fluid pressure and velocity conditions on static and dynamic embedded interfaces. This method is based on the exact solution of local, one‐dimensional, fluid–structure Riemann problems. Next, it describes two consistent and conservative approaches for computing the flow‐induced loads on rigid and flexible embedded structures. The first approach reconstructs the interfaces within the CFD solver. The second one represents them as zero level sets, and works instead with surrogate fluid/structure interfaces. For example, the surrogate interfaces obtained simply by joining contiguous segments of the boundary surfaces of the fluid control volumes that are the closest to the zero level sets are explored in this work. All numerical algorithms presented in this paper are applicable with any embedding CFD mesh, whether it is structured or unstructured. Their performance is illustrated by their application to the solution of three‐dimensional fluid–structure interaction problems associated with the fields of aeronautics and underwater implosion. Copyright © 2011 John Wiley & Sons, Ltd.     

An improved Rhie–Chow interpolation scheme for the smoothed‐interface immersed boundary method

Rhie–Chow interpolation is a commonly used method in CFD calculations on a co‐located mesh in order to suppress non‐physical pressure oscillations arising from chequerboard effects. A fully parallelized smoothed‐interface immersed boundary method on a co‐located grid is described in this paper. We discuss the necessity of modifications to the original Rhie–Chow interpolation in order to deal with a locally refined mesh. Numerical simulation with the modified scheme of Choi shows that numerical dissipation due to Rhie–Chow interpolation introduces significant errors at the immersed boundary. To address this issue, we develop an improved Rhie–Chow interpolation scheme that is shown to increase the accuracy in resolving the flow near the immersed boundary. We compare our improved scheme with the modified scheme of Choi by parallel simulations of benchmark flows: (i) flow past a stationary cylinder; (ii) flow past an oscillating cylinder; and (iii) flow past a stationary elliptical cylinder, where Reynolds numbers are tested in the range 10–200. Our improved scheme is significantly more accurate and compares favourably with a staggered grid algorithm. We also develop a scheme to compute the boundary force for the direct‐forcing immersed boundary method efficiently. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of multiple shock–shock interference using implicit anti‐diffusive WENO schemes

Accurate computations of two‐dimensional turbulent hypersonic shock–shock interactions that arise when single and dual shocks impinge on the bow shock in front of a cylinder are presented. The simulation methods used are a class of lower–upper symmetric‐Gauss–Seidel implicit anti‐diffusive weighted essentially non‐oscillatory (WENO) schemes for solving the compressible Navier–Stokes equations with Spalart–Allmaras one‐equation turbulence model. A numerical flux of WENO scheme with anti‐diffusive flux correction is adopted, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of first‐order dissipative methods. Experimental flow fields of type IV shock–shock interactions with single and dual incident shocks by Wieting are computed. By using the WENO scheme with anti‐diffusive flux corrections, the present solution indicates that good accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Computed surface pressure distribution and heat transfer rate are also compared with experimental data and other computational results and good agreement is found. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of two‐phase flow with moving immersed boundaries

A two‐dimensional multi‐phase model for immiscible binary fluid flow including moving immersed objects is presented. The fluid motion is described by the incompressible Navier–Stokes equation coupled with a phase‐field model based on van der Waals' free energy density and the Cahn–Hilliard equation. A new phase‐field boundary condition was implemented with minimization of the free energy in a direct way, to specifically improve the physical behavior of the contact line dynamics for moving immersed objects. Numerical stability and execution time were significantly improved by the use of the new boundary condition. Convergence toward the analytical solution was demonstrated for equilibrium contact angle, the Lucas–Washburn theory and Stefan's problem. The proposed model may be used for multi‐phase flow problems with moving boundaries of complex geometry, such as the penetration of fluid into a deformable, porous medium. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite volume solution of the Navier–Stokes equations in velocity–vorticity formulation

For the incompressible Navier–Stokes equations, vorticity‐based formulations have many attractive features over primitive‐variable velocity–pressure formulations. However, some features interfere with the use of the numerical methods based on the vorticity formulations, one of them being the lack of a boundary conditions on vorticity. In this paper, a novel approach is presented to solve the velocity–vorticity integro‐differential formulations. The general numerical method is based on standard finite volume scheme. The velocities needed at the vertexes of each control volume are calculated by a so‐called generalized Biot–Savart formula combined with a fast summation algorithm, which makes the velocity boundary conditions implicitly satisfied by maintaining the kinematic compatibility of the velocity and vorticity fields. The well‐known fractional step approaches are used to solve the vorticity transport equation. The paper describes in detail how we accurately impose no normal‐flow and no tangential‐flow boundary conditions. We impose a no‐flux boundary condition on solid objects by the introduction of a proper amount of vorticity at wall. The diffusion term in the transport equation is treated implicitly using a conservative finite update. The diffusive fluxes of vorticity into flow domain from solid boundaries are determined by an iterative process in order to satisfy the no tangential‐flow boundary condition. As application examples, the impulsively started flows through a flat plate and a circular cylinder are computed using the method. The present results are compared with the analytical solution and other numerical results and show good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical solution of single and multi‐phase internal transonic flow problems

This paper presents numerical methods for solving turbulent and two‐phase transonic flow problems. The Navier–Stokes equations are solved using cell‐vertex Lax–Wendroff method with artificial dissipation or cell‐centred upwind method with Roe's Riemann solver and linear reconstruction. Due to a big difference of time scales in two‐phase flow of condensing steam a fractional step method is used. Test cases including 2D condensing flow in a nozzle and one‐phase transonic flow in a turbine cascade with transition to turbulence are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

A high‐resolution method for compressible two‐fluid flows and simulation of three‐dimensional shock–bubble interactions

A high‐resolution method is developed to capture the material interfaces of compressible two‐fluid flows in multiple dimensions. A fluid mixture model system with single velocity and pressure is used, and viscous effect can also be taken into account. A consistent thermodynamic law based on the assumption of pressure equilibrium is employed to describe the thermodynamic behaviors of the pure fluids and mixture of two components. The splitting and unsplit Eulerian formulations of piecewise parabolic method are extended to numerically integrate the hyperbolic part of the model system, whereas the system of diffusion equations is solved using an explicit, central difference scheme. The block‐structured adaptive mesh refinement (AMR) capability is built in the hydrodynamic code to locally improve grid resolution. The resulting method is verified to be at least second‐order accurate in space. Numerical results show that the discontinuities, particularly contact discontinuities, can be resolved sharply. The use of AMR allows flow features at disparate scales to be resolved sufficiently. In addition, three‐dimensional shock–bubble interactions are simulated to investigate effects of Mach number on bubble evolution. The flow structures including those peculiar to three‐dimensional bubble are resolved correctly, and some physical phenomena with increasing Mach number are reported. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of vortex shedding past a square cylinder with different turbulence models

This paper presents the results of numerical simulations of vortex shedding past a free-standing square cylinder at Re_D=22 000, obtained with different turbulence models. Using wall functions, the standard k–ε model is compared with a modification suggested by Kato and Launder (Proc. 9th Symp. Turbulent Shear Flows, Kyoto, 10-4-1 (1993)). In addition, both versions are used in a two-layer approach, in which the flow close to the cylinder is computed with a locally more suitable one-equation turbulence model and only outside the viscous near-wall layer with the two mentioned high-Re model versions. To allow a comparison, the simulations are performed first using the same computational domain and boundary conditions as in previous investigations. Then results are presented that were obtained on a computational domain and with boundary conditions more suitable for a comparison with the experiments. © 1998 John Wiley & Sons, Ltd.     

Numerical simulation of viscous flow interaction with an elastic membrane

A numerical fluid–structure interaction model is developed for the analysis of viscous flow over elastic membrane structures. The Navier–Stokes equations are discretized on a moving body‐fitted unstructured triangular grid using the finite volume method, taking into account grid non‐orthogonality, and implementing the SIMPLE algorithm for pressure solution, power law implicit differencing and Rhie–Chow explicit mass flux interpolations. The membrane is discretized as a set of links that coincide with a subset of the fluid mesh edges. A new model is introduced to distribute local and global elastic effects to aid stability of the structure model and damping effects are also included. A pseudo‐structural approach using a balance of mesh edge spring tensions and cell internal pressures controls the motion of fluid mesh nodes based on the displacements of the membrane. Following initial validation, the model is applied to the case of a two‐dimensional membrane pinned at both ends at an angle of attack of 4° to the oncoming flow, at a Reynolds number based on the chord length of 4 × 10³. A series of tests on membranes of different elastic stiffness investigates their unsteady movements over time. The membranes of higher elastic stiffness adopt a stable equilibrium shape, while the membrane of lowest elastic stiffness demonstrates unstable interactions between its inflated shape and the resulting unsteady wake. These unstable effects are shown to be significantly magnified by the flexible nature of the membrane compared with a rigid surface of the same average shape. Copyright © 2007 John Wiley & Sons, Ltd.     

Computation of head–disk interface gap micro flowfields using DSMC and continuum–atomistic hybrid methods

This paper discusses computational modeling of micro flow in the head–disk interface (HDI) gap using the direct simulation Monte Carlo (DSMC) method. Modeling considerations are discussed in detail both for a stand‐alone DSMC computation and for the case of a hybrid continuum–atomistic simulation that couples the Navier–Stokes (NS) equation to a DSMC solver. The impact of the number of particles and number of cells on the accuracy of a DSMC simulation of the HDI gap is investigated both for two‐ and three‐dimensional configurations. An appropriate implicit boundary treatment method for modeling inflow and outflow boundaries is used in this work for a three‐dimensional DSMC micro flow simulation. As the flow outside the slider is in the continuum regime, a hybrid continuum–atomistic method based on the Schwarz alternating method is used to couple the DSMC model in the slider bearing region to the flow outside the slider modeled by NS equation. Schwarz coupling is done in two dimensions by taking overlap regions along two directions and the Chapman–Enskog distribution is employed for imposing the boundary condition from the continuum region to the DSMC region. Converged hybrid flow solutions are obtained in about five iterations and the hybrid DSMC–NS solutions show good agreement with the exact solutions in the entire domain considered. An investigation on the impact of the size of the overlap region on the convergence behavior of the Schwarz method indicates that the hybrid coupling by the Schwarz method is weakly dependent on the size of the overlap region. However, the use of a finite overlap region will facilitate the exchange of boundary conditions as the hybrid solution has been found to diverge in the absence of an overlap region for coupling the two models. Copyright © 2009 John Wiley & Sons, Ltd.     

Role of the inclination of an inverted‐V upper plate on the heat and flow behavior of trapped gases in a modified Rayleigh–Bénard cavity

Thermal buoyant air inside a modified Rayleigh–Bénard (RB) cavity bounded by a lower flat plate and an inverted‐V upper plate has been investigated numerically using the finite‐volume method. The second‐order‐accurate QUICK and SIMPLE schemes were used for the discretization of the convective terms and the pressure–velocity coupling in the set of conservation equations, respectively. The problem under study is controlled by two parameters: (1) the Rayleigh number ranging from 10³ to 10⁶ and (2) the relative height of the vertical sidewalls d. In reference to the latter, it varies from one limiting case corresponding to the standard RB cavity (a rectangle with d = 1) to another limiting case represented by an isosceles triangular cavity where d = 0. The numerical results for the velocity and temperature fields are presented in terms of streamlines, isotherms, local and mean heat fluxes. An additional effort was devoted to determine the critical Ra values characterizing the transition from symmetrical to asymmetrical buoyant airflow responsive to incremental changes in Ra. For purposes of engineering design, a general correlation equation for the Nusselt number in terms of the pertinent Ra and d was constructed using nonlinear multiple regression theory. Copyright © 2007 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Computation of unsteady transonic flow over a fighter wing using a zonal Navier–Stokes/full‐potential method

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the inner zone, the Navier–Stokes equations are solved using a diagonal form of an alternating‐direction implicit (ADI) approximate factorisation procedure. In the outer zone, the unsteady full‐potential equation (FPE) is solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic‐based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are about 60% of the full Navier–Stokes CPU times for unsteady flows in non‐vector processing machines. Applications of the method are presented for a F‐5 wing in steady and unsteady transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to the full Navier–Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier–Stokes equations can be retained with significant savings in computational time. Copyright © 1999 John Wiley & Sons, Ltd.