SUMMARYWe describe the implementation of a computational fluid dynamics solver for the simulation of high-speed flows. It comprises a finite volume (FV) discretization using semi-discrete, non-staggered central schemes for colocated variables prescribed on a mesh of polyhedral cells that have an arbitrary number of faces. We describe the solver in detail, explaining the choice of variables whose face interpolation is limited, the choice of limiter, and a method for limiting the interpolation of a vector field that is independent of the coordinate system. The solution of momentum and energy transport in the Navier-Stokes equations uses an operator-splitting approach: first, we solve an explicit predictor equation for the convection of conserved variables, then an implicit corrector equation for the diffusion of primitive variables. Our solver is validated against four sets of data: (1) an analytical solution of the one-dimensional shock tube case; (2) a numerical solution of two dimensional, transient, supersonic flow over a forward-facing step; (3) interferogram density measurements of a supersonic jet from a circular nozzle; and (4) pressure and heat transfer measurements in hypersonic flow over a 25 • -55 • biconic. Our results indicate that the central-upwind scheme of Kurganov, Noelle and Petrova (SIAM J. Sci. Comput. 2001; 23:707-740) is competitive with the best methods previously published (e.g. piecewise parabolic method, Roe solver with van Leer limiting) and that it is inherently simple and well suited to a colocated, polyhedral FV framework.

Brenner has recently proposed modifications to the Navier-Stokes equations that are based on theoretical arguments but supported only by experiments having a fairly limited range [1,2]. These modifications relate to a diffusion of fluid volume that would be significant for flows with high density gradients. So the viscous structure of shock waves in gases should provide an excellent test case for this new model. In this paper we detail the shock structure problem and propose exponents for the gas viscosity-temperature relation based on empirical viscosity data that is independent of shock experiments. We then simulate shocks in the range Mach 1.0-12.0 using the Navier-Stokes equations, both with and without Brenner's modifications. Initial simulations showed Brenner's modifications display unphysical behaviour when the coefficient of volume diffusion exceeds the kinematic viscosity. Our subsequent analyses attribute this behaviour to both an instability to temporal disturbances and a spurious phase velocity-frequency relationship. On equating the volume diffusivity to the kinematic viscosity, however, we find the results with Brenner's modifications are significantly better than those of the standard Navier-Stokes equations, and broadly similar to those from the family of extended hydrodynamic models that includes the Burnett equations. Brenner's modifications add only two terms to the Navier-Stokes equations, and the numerical implementation is much simpler than conventional extended hydrodynamic models, particularly in respect of boundary conditions. We recommend further investigation and testing on a number of different benchmark non-equilibrium flow cases.

SUMMARYThis paper outlines the development of a new procedure for analysing continuum mechanics problems with a particular focus on fluid-structure interaction in flexible tubes. A review of current methods of fluid-structure coupling highlights common limitations of high computational cost and solution instability. It is proposed that these limitations can be overcome by an alternative approach in which both fluid and solid components are solved within a single discretized continuum domain. A single system of momentum and continuity equations is therefore derived that governs both fluids and solids and which are solved with a single mesh using finite volume discretization schemes. The method is validated first by simulating dynamic oscillation of a clamped elastic beam. It is then applied to study the case of interest-wave propagation in highly flexible tubes-in which a predicted wave speed of 8.58 m/s falls within 2% of an approximate analytical solution. The method shows further good agreement with analytical solutions for tubes of increasing rigidity, covering a range of wave speeds from those found in arteries to that in the undisturbed fluid.

A new Computational Fluid Dynamics (CFD) solver for high-speed viscous §ows in the OpenFOAM code is validated against published experimental data and Direct Simulation Monte Carlo (DSMC) results. The laminar §at plate and circular cylinder cases are studied for Mach numbers, Ma, ranging from 6 to 12.7, and with argon and nitrogen as working gases. Simulation results for the laminar §at plate cases show that the combination of accommodation coe©cient values σ u = 0.7 and σ T = 1.0 in the Maxwell/Smoluchowski conditions, and the coe©cient values A 1 = 1.5 and A 2 = 1.0 in the second-order velocity slip condition, give best agreement with experimental data of surface pressure. The values σ u = 0.7 and σ T = 1.0 also give good agreement with DSMC data of surface pressure at the stagnation point in the circular cylinder case at Kn = 0.25. The Langmuir surface adsorption condition is also tested for the laminar §at plate case, but initial results were not as good as the Maxwell/Smoluchowski boundary conditions.

This paper investigates the use of Navier-Stokes-Fourier equations with non-equilibrium boundary conditions (BCs) for simulation of rarefied hypersonic flows. It revisits a largely forgotten derivation of velocity slip and temperature jump by Patterson, based on Grad's moment method. Mach 10 flow around a cylinder and Mach 12.7 flow over a flat plate are simulated using both computational fluid dynamics using the temperature jump BCs of Patterson and Smoluchowski and the direct simulation Monte-Carlo (DMSC) method. These flow exhibit such strongly non-equilibrium behaviour that, following Patterson's analysis, they are strictly beyond the range of applicability of the BCs. Nevertheless, the results using Patterson's temperature jump BC compare quite well with the DSMC and are consistently better than those using the standard Smoluchowski temperature jump BC. One explanation for this better performance is that an assumption made by Patterson, based on the flow being only slightly non-equilibrium, introduces an additional constraint to the resulting BC model in the case of highly non-equilibrium flows.