Face- and Cell-Averaged Nodal-Gradient Approach to Cell-Centered Finite-Volume Method on Mixed Grids

In this paper, the averaged nodal-gradient approach previously developed for triangular grids is extended to mixed triangular-quadrilateral grids. It is shown that the face- averaged approach leads to deteriorated iterative convergence on quadrilateral grids. To develop a convergent solver, we consider cell-averaging instead of face-averaging for quadri- lateral cells. We show that the cell-averaged approach leads to a convergent solver and can be efficiently combined with the face-averaged approach on mixed grids. The method is demonstrated for various inviscid and viscous problems from low to high Mach numbers on two-dimensional mixed grids

Efficient and Robust Weighted Least-Squares Cell-Average Gradient Construction Methods for the Simulation of Scramjet Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches designed to minimize the number of cells used to augment/stabilize the least-squares stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cellaverage gradients is described. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet engine flowpath is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the preferred gradient method for a realistic 3-D geometry on a non-hex-dominant grid

An Efficient Quadratic Interpolation Scheme for a Third-Order Cell-Centered Finite-Volume Method on Tetrahedral Grids

In this paper, we propose an efficient quadratic interpolation formula utilizing solution gradients computed and stored at nodes and demonstrate its application to a third-order cell-centered finite-volume discretization on tetrahedral grids. The proposed quadratic formula is constructed based on an efficient formula of computing a projected derivative. It is efficient in that it completely eliminates the need to compute and store second derivatives of solution variables or any other quantities, which are typically required in upgrading a second-order cell-centered unstructured-grid finite-volume discretization to third-order accuracy. Moreover, a high-order flux quadrature formula, as required for third-order accuracy, can also be simplified by utilizing the efficient projected-derivative formula, resulting in a numerical flux at a face centroid plus a curvature correction not involving second derivatives of the flux. Similarly, a source term can be integrated over a cell to high-order in the form of a source term evaluated at the cell centroid plus a curvature correction, again, not requiring second derivatives of the source term. The discretization is defined as an approximation to an integral form of a conservation law but the numerical solution is defined as a point value at a cell center, leading to another feature that there is no need to compute and store geometric moments for a quadratic polynomial to preserve a cell average. Third-order accuracy and improved second-order accuracy are demonstrated and investigated for simple but illustrative test cases in three dimensions

Boundary Layer Transition over Blunt Hypersonic Vehicles Including Effects of Ablation-Induced Out-Gassing

Computations are performed to study the boundary layer instability mechanisms pertaining to hypersonic flow over blunt capsules. For capsules with ablative heat shields, transition may be influenced both by out-gassing associated with surface pyrolysis and the resulting modification of surface geometry including the formation of micro-roughness. To isolate the effects of out-gassing, this paper examines the stability of canonical boundary layer flows over a smooth surface in the presence of gas injection into the boundary layer. For a slender cone, the effects of out-gassing on the predominantly second mode instability are found to be stabilizing. In contrast, for a blunt capsule flow dominated by first mode instability, out-gassing is shown to be destabilizing. Analogous destabilizing effects of outgassing are also noted for both stationary and traveling modes of crossflow instability over a blunt sphere-cone configuration at angle of attack

Very unusual case of a primary sinonasal germ cell tumour

Sinonasal malignancies are a very rare diagnosis. We present a unique case of a 32-year-old man who presented with symptoms of worsening sinusitis and periorbital cellulitis. Investigation found a sinonasal malignancy and pathology confirmed this to be a primary germ cell tumour. The patient was managed with chemotherapy, surgery and consolidation radiotherapy and has remained well to date. This case report outlines an unusual presentation and diagnostic challenge for the primary care physician, ear, nose and throat surgeon, pathologist and oncologist with review of the surrounding literature

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved

Analysis of Instabilities in Non-Axisymmetric Hypersonic Boundary Layers Over Cones

Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above

Secondary Instability of Second Modes in Hypersonic Boundary Layers

Second mode disturbances dominate the primary instability stage of transition in a number of hypersonic flow configurations. The highest amplification rates of second mode disturbances are usually associated with 2D (or axisymmetric) perturbations and, therefore, a likely scenario for the onset of the three-dimensionality required for laminar-turbulent transition corresponds to the parametric amplification of 3D secondary instabilities in the presence of 2D, finite amplitude second mode disturbances. The secondary instability of second mode disturbances is studied for selected canonical flow configurations. The basic state for the secondary instability analysis is obtained by tracking the linear and nonlinear evolution of 2D, second mode disturbances using nonlinear parabolized stability equations. Unlike in previous studies, the selection of primary disturbances used for the secondary instability analysis was based on their potential relevance to transition in a low disturbance environment and the effects of nonlinearity on the evolution of primary disturbances was accounted for. Strongly nonlinear effects related to the self-interaction of second mode disturbances lead to an upstream shift in the upper branch neutral location. Secondary instability computations confirm the previously known dominance of subharmonic modes at relatively small primary amplitudes. However, for the Purdue Mach 6 compression cone configuration, it was shown that a strong fundamental secondary instability can exist for a range of initial amplitudes of the most amplified second mode disturbance, indicating that the exclusive focus on subharmonic modes in the previous applications of secondary instability theory to second mode primary instability may not have been fully justified

Investor Perceptions of Board Performance: Evidence From Uncontested Director Elections

This paper provides evidence that uncontested director elections provide informative polls of investor perceptions regarding board performance. We find that higher (lower) vote approval is associated with lower (higher) stock price reactions to subsequent announcements of management turnovers. In addition, firms with low vote approval are more likely to experience CEO turnover, greater board turnover, lower CEO compensation, fewer and better-received acquisitions, and more and better-received divestitures in the future. These findings hold after controlling for other variables reflecting or determining investor perceptions, suggesting that elections not only inform as a summary statistic, but incrementally inform as well

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling