On air-chemistry reduction for hypersonic external flow applications

In external hypersonic flows, viscous and compressibility effects generate very high temperatures leading to significant chemical reactions among air constituents. Therefore, hypersonic flow computations require coupled calculations of flow and chemistry. Accurate and efficient computations of air-chemistry kinetics are of much importance for many practical applications but calculations accounting for detailed chemical kinetics can be prohibitively expensive. In this paper, we investigate the possibility of applying chemical kinetics reduction schemes for hypersonic air-chemistry. We consider two chemical kinetics sets appropriate for three different temperature ranges: 2500 K to 4500 K; 4500 K to 9000 K; and above 9000 K. By demonstrating the existence of the so-called the slow manifold in each of the chemistry sets, we show that judicious chemical kinetics reduction leading to significant computational savings is possible without much loss in accuracy.     

Efficient approximation of Sparse Jacobians for time‐implicit reduced order models

This paper introduces a sparse matrix discrete interpolation method to effectively compute matrix approximations in the reduced order modeling framework. The sparse algorithm developed herein relies on the discrete empirical interpolation method and uses only samples of the nonzero entries of the matrix series. The proposed approach can approximate very large matrices, unlike the current matrix discrete empirical interpolation method, which is limited by its large computational memory requirements. The empirical interpolation indices obtained by the sparse algorithm slightly differ from the ones computed by the matrix discrete empirical interpolation method as a consequence of the singular vectors round‐off errors introduced by the economy or full singular value decomposition (SVD) algorithms when applied to the full matrix snapshots. When appropriately padded with zeros, the economy SVD factorization of the nonzero elements of the snapshots matrix is a valid economy SVD for the full snapshots matrix. Numerical experiments are performed with the 1D Burgers and 2D shallow water equations test problems where the quadratic reduced nonlinearities are computed via tensorial calculus. The sparse matrix approximation strategy is compared against five existing methods for computing reduced Jacobians: (i) matrix discrete empirical interpolation method, (ii) discrete empirical interpolation method, (iii) tensorial calculus, (iv) full Jacobian projection onto the reduced basis subspace, and (v) directional derivatives of the model along the reduced basis functions. The sparse matrix method outperforms all other algorithms. The use of traditional matrix discrete empirical interpolation method is not possible for very large dimensions because of its excessive memory requirements. Copyright © 2016 John Wiley & Sons, Ltd.     

High-order discontinuous Galerkin method for applications to multicomponent and chemically reacting flows

This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chem-ically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic proper-ties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicompo-nent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.     

Development of a lattice Boltzmann method for two‐layered shallow‐water flow

In this paper the dynamics of a two‐layered liquid, made of two immiscible shallow‐layers of different density, has been investigated within the framework of the lattice Boltzmann method (LBM). The LBM developed in this paper for the two‐layered, shallow‐water flow has been obtained considering two separate sets of LBM equations, one for each layer. The coupling terms between the two sets have been defined as external forces, acted on one layer by the other. Results obtained from the LBM developed in this paper are compared with numerical results obtained solving the two‐layered, shallow‐water equations, with experimental and other numerical results published in literature. The results are interesting. First, the numerical results obtained by the LBM and by the shallow‐water model can be considered as equivalent. Second, the LBM developed in this paper is able to simulate motion conditions on nonflat topography. Third, the agreement between the LBM (and also shallow‐water model) numerical results and the experimental results is good when the evolution of the flow does not depend on the viscosity, that is, during the initial phase of the flow, dominated by gravity and inertia forces. When the viscous forces dominate the evolution of the flow the agreement between numerical and experimental results depends strongly on the viscosity; it is good if the numerical LBM viscosity has the same order of magnitude of the liquid's kinematic viscosity. Copyright © 2012 John Wiley & Sons, Ltd.     

Reduced‐order controllers for control of flow past an airfoil

Reduced‐order controller design by means of reduced‐order model for control of a wake flow is presented. Reduced‐order model is derived by combining the Galerkin projection with proper orthogonal decomposition (POD) or with other related reduced‐order approaches such as singular value decomposition or reduced‐basis method. In the present investigation, we discuss the applicability of the reduced‐order approaches for fast computation of the optimal control for control of vortex shedding behind a thin airfoil through unsteady blowing on the airfoil surface. Accuracy of the reduced‐order model is quantified by comparing flow fields obtained from the reduced‐order models with those from the full‐order simulations under the same free‐stream conditions. A control of vortex shedding is demonstrated for Reynolds number 100. It is found that downstream directed blowing on the upper surface of the airfoil near the leading edge is more efficient in mitigating flow separation and suppressing the vortex shedding. Copyright © 2005 John Wiley & Sons, Ltd.     

Immersed boundary conditions method for unsteady flow problems described by the Laplace operator

An implicit, spectral algorithm for the analysis of unsteady flow problems governed by the Laplace operator in corrugated geometries is described. The algorithm treats the physical boundary conditions as constraints along lines internal to the solution domain. The method eliminates the need for coordinate generation and can be quickly adapted to changing geometries. Various tests confirm the spectral accuracy in space and the first‐ and second‐order accuracies in time. Copyright © 2007 John Wiley & Sons, Ltd.     

The numerical solution of the thin film flow surrounding a horizontal cylinder resulting from a vertical cylindrical jet

The numerical solution of the thin film flow surrounding a horizontal cylinder resulting from a single vertical cylindrical jet is obtained. This is effected by transforming the domain of the flow, which contains a free surface, onto a rectangular parallelepiped and using a marching strategy to solve the ensuing parabolic equations. The flow terminates at a finite distance along the cylinder, its position depending on the velocity and mass flux of the jet. A comparison with the usual two-dimensional model in which the jet is replaced by a vertical sheet shows that such a representation is valid provided the overall width of the flow is not too large. In particular, the differences in heat transfer characteristics amount to a few per cent, thus validating the use of the two-dimensional model when applied to heat exchanger tubes. A comparison with the more usual multicolumn case is also considered.     

Parallel solution of a coupled flow and transport model for shallow water

We present an integrated approach for the concurrent solution of a 3D hydrodynamical model coupled with a 3D transport model. Since both models are quite similar in nature, the same numerical method has been employed. This leads to a code that is more efficient than when two existing codes would have been combined. Discretization of the spatial differential operators, and the boundary conditions, results in a stiff initial value problem. To cope with the stiffness, we select an implicit time‐integration formula, viz. the second‐order, L‐stable BDF method because of its excellent stability properties. To reduce the huge amount of linear algebra involved in solving the implicit relations, an Approximate Factorization technique has been used. Essentially, this technique replaces a ‘multi‐dimensional’ system by a series of ‘one‐dimensional’ systems. Since the output of the hydrodynamical model (i.e., the flow field) serves as input for the transport model, we solve the hydrodynamical model one time step ahead in time. This allows us to solve the models in parallel, using two different groups of processors. By a little tuning of the parameters in the algorithm, a load‐balancing has been obtained that is close to optimal. As a result, both models require roughly the same amount of CPU time, so that one of them, effectively, can be considered as obtained ‘for free’. Copyright © 2002 John Wiley & Sons, Ltd.     

White in Time Scalar Advection Model as a Tool for Solving Joint Composition PDF Equations

A rapidly decorrelating velocity field model is used to derive stochastic partial differential equations (SPDE) allowing one to compute the modeled one-point joint probability density function of turbulent reactive scalars. Those SPDEs are shown to be hyperbolic advection/reaction equations. They are dealt with in a generalized sense, so that discontinuities in the scalar fields can be treated. The Eulerian Monte Carlo (EMC) method thus defined is coupled with a RANS solver and applied to the computation of a turbulent premixed methane flame over a backward facing step.     

An approach formulated in terms of conserved variables for the characterisation of propellant combustion in internal ballistics

A model formulated in terms of conserved variables is proposed for its use in the study of internal ballistic problems of pyrotechnical mixtures and propellants. It is a transient two‐phase flow model adapted from the non‐conservative Gough model. This conversion is mathematically attractive because of the wide range of numerical methods for this kind of systems that may be applied. We propose the use of the AUSM+, AUSM + up and Rusanov schemes as an efficient alternative for this type of two‐phase problem. A splitting technique is applied, which solves the system of equations in several steps. A second‐order approach based on Monotonic Upstream‐Centred Scheme for Conservation Laws (MUSCL) is also used. Some tests are used to validate the code, namely a shock wave test, a contact discontinuity problem and an internal ballistics problem. In this last case, one‐dimensional numerical results are compared with experimental data of 155‐mm gunshots. Copyright © 2015 John Wiley & Sons, Ltd.     

Nonisothermal two‐ and three‐dimensional flow simulations of inelastic and viscoelastic fluids by a finite‐volume method

This paper applies the finite‐volume method to computations of steady flows of viscous and viscoelastic incompressible fluids in complex two and three‐dimensional geometries. The materials adopted in the study obey different constitutive laws: Newtonian, purely viscous Carreau–Yasuda as also Upper‐Convected Maxwell and Phan‐Thien/Tanner differential models, with a Williams–Landel–Ferry (WLF) equation for temperature dependence. Specific analyses are made depending on the rheological model. A staggered grid is used for discretizing the equations and unknowns. Stockage possibilities allow us to solve problems involving a great number of degrees of freedom, up to 1 500 000 unknowns with a desk computer. In relation to the fluid properties, our numerical simulations provide flow characteristics for various 2D and 3D configurations and demonstrate the possibilities of the code to solve problems involving complex nonlinear constitutive equations with thermal effects. Copyright © 2009 John Wiley & Sons, Ltd.     

An inhomogeneous three-phase model for the flow in aluminium reduction cells

The exact prediction of both the melts flow and metal–bath interface deformation is critical to the design of commercial aluminium reduction cells which operate steadily at high current efficiency (CE). An inhomogeneous flow model of three phases, metal, bath and gas bubbles, was established in this study. The flow in a 300 kA cell was numerically computed with the model by using the finite volume scheme, and the results reveal the relation among the motion of the three phases and prove that the globality of the flow in cells could not be neglected. An improved calculation model of CE, based on the surface renewal theory, was also developed by connecting the melts flow with the mass transfer of reduced entity at the metal–bath interface. Using the model, the local CE in the 300 kA cell was predicted, and the calculated CE value of the whole cell is reasonable.     

A note on probabilistic interpretation for quasilinear mixed boundary problems

I.Intr0ducti0nWekn0wthattheprobabilisticrepresentationofsolutionsofpartialdifferentia1equationswiththemixedboundaryconditi0nshasmanyimp0rtantapplicationsbothintheoryofpartialdifferentialequationsandinthat0fstochasticdifferentialequationswithreflectingboun…     

Differential characteristic set algorithm for the complete symmetry classification of partial differential equations

In this paper, we present a differential polynomial characteristic set algorithm for the complete symmetry classification of partial differential equations （PDEs） with some parameters. It can make the solution to the complete symmetry classification problem for PDEs become direct and systematic. As an illustrative example, the complete potential symmetry classifications of nonlinear and linear wave equations with an arbitrary function parameter are presented. This is a new application of the differential form characteristic set algorithm, i.e., Wu＇s method, in differential equations.     

Flooding in urban drainage systems: coupling hyperbolic conservation laws for sewer systems and surface flow

In this paper, we propose a model for a sewer network coupled to surface flow and investigate it numerically. In particular, we present a new model for the manholes in storm sewer systems. It is derived using the balance of the total energy in the complete network. The resulting system of equations contains, aside from hyperbolic conservation laws for the sewer network and algebraic relations for the coupling conditions, a system of ODEs governing the flow in the manholes. The manholes provide natural points for the interaction of the sewer system and the runoff on the urban surface modeled by shallow‐water equations. Finally, a numerical method for the coupled system is presented. In several numerical tests, we study the influence of the manhole model on the sewer system and the coupling with 2D surface flow. Copyright © 2014 John Wiley & Sons, Ltd.     

A two‐dimensional numerical scheme of dry/wet fronts for the Saint‐Venant system of shallow water equations

We propose a new two‐dimensional numerical scheme to solve the Saint‐Venant system of shallow water equations in the presence of partially flooded cells. Our method is well balanced, positivity preserving, and handles dry states. The latter is ensured by using the draining time step technique in the time integration process, which guarantees non‐negative water depths. Unlike previous schemes, our technique does not generate high velocities at the dry/wet boundaries, which are responsible for small time step sizes and slow simulation runs. We prove that the new scheme preserves ‘lake at rest’ steady states and guarantees the positivity of the computed fluid depth in the partially flooded cells. We test the new scheme, along with another recent scheme from the literature, against the analytical solution for a parabolic basin and show the improved simulation performance of the new scheme for two real‐world scenarios. Copyright © 2014 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

Accuracy and efficiency of two numerical methods of solving the potential flow problem for highly nonlinear and dispersive water waves

The accuracy and efficiency of two methods of resolving the exact potential flow problem for nonlinear waves are compared using three different one horizontal dimension (1DH) test cases. The two model approaches use high‐order finite difference schemes in the horizontal dimension and differ in the resolution of the vertical dimension. The first model uses high‐order finite difference schemes also in the vertical, while the second model applies a spectral approach. The convergence, accuracy, and efficiency of the two models are demonstrated as a function of the temporal, horizontal, and vertical resolutions for the following: (1) the propagation of regular nonlinear waves in a periodic domain; (2) the motion of nonlinear standing waves in a domain with fully reflective boundaries; and (3) the propagation and shoaling of a train of waves on a slope. The spectral model approach converges more rapidly as a function of the vertical resolution. In addition, with equivalent vertical resolution, the spectral model approach shows enhanced accuracy and efficiency in the parameter range used for practical model applications. Copyright © 2015 John Wiley & Sons, Ltd.