Search NASASearch

Engineering topics

Murman, Scott

Publications and source records attributed to Murman, Scott.

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined.

Franciolini, Matteo

Scale-Resolving Simulations of Laminar-to-Turbulent Transition in Swept-Wing Flow

Aerodynamic shape optimization of swept Natural Laminar Flow (NLF) wing design plays an important role for the development of future energy-efficient commercial transport aircraft. Lack of predictive laminar-to-turbulent transition modeling capabilities is still limiting the extent to which design can rely on computational tools, while making costly wind tunnel and flight tests a necessity. Prediction of transition on practical 3D wing geometries using the e (sup N) method coupled with laminar boundary layer and linear stability codes or local statistical correlation CFD (Computational Fluid Dynamics) models such as the gamma - Re (sub theta) Langtry-Menter model have proven challenging, whereas transitional hybrid RANS-LES solvers have shown encouraging but still mitigated results.

Scale-resolving

Preliminary Scale-Resolving Simulations of Laminar-To-Turbulent Transition in Swept-Wing Flow

Swept wings and control surfaces are common elements of modern aircraft, and it has been shown both experimentally and theoretically that laminar-to-turbulent transition of the three-dimensional boundary layer that develops over them is highly sensitive to surface roughness. Numerous studies have been conducted on the effect of discrete roughness elements or distributed roughness elements on swept flow transition, however so far limited computational effort has been dedicated to the study of transition over swept wings with randomly distributed micron-sized roughness. In the present work, we set up to reproduce the extensive experimental data base generated by Dagenhart et al for the infinite swept wing NLF(2)-0415. To this purpose, we perform scale-resolving simulations of flow transition over smooth and rough surfaces using a high-order space-time spectral-element Discontinuous-Galerkin solver. Different types of surface roughnesses are implemented by elastically deforming the original mesh. The study shows that the experimental results cannot be accounted for by a perfectly smooth wing and reveals a strong sensitivity of the transition process to the representation of the surface roughness. The crossflow patterns and transition location approach those measured for some of the surface profiles, however a correlation between the wavenumber spectrum of the surface, grid resolution and boundary layer stability is yet to be established.

Denison, Marie

Robust Metric-Aligned Quad-Dominant Meshing Using L(sub p) Centroidal Voronoi Tessellation

We introduce a meshing algorithm that can be used to both generate and adapt meshes for bounded domains in an anisotropic manner. This is particularly beneficial when anisotropic flow features like shock waves or contact discontinuities are present in the computational domain. The algorithm presented in this paper is based upon meshing under the imposed Riemannian metric tensor, which controls the orientation and size of the mesh elements. In this way there is no need for user intervention to recognize these features. We demonstrate that the method indeed aligns the elements with the underlying metric and produces right-angled simplices that can be recombined into quadrilateral elements. The aim is to eventually incorporate this meshing strategy in the monolithic high-order spectral element solver that is currently being developed at NASA Ames. This paper has two main contributions: First, we demonstrate that we can generate quad-dominant metric-aligned meshes for bounded domains using a generalized form of L(sub p)-Centroidal Voronoi Tessellation (L(sub p)-CVT). Unlike previous works, we do not rely on a background mesh and discretize the bounded domain in a hierarchical way by first discretizing the boundaries and then the volume using the underlying metric. Second, we present an alternative for clipping the Voronoi cells on the boundary, which is common practice in CVT-based meshing algorithms, by reconstructing the Voronoi cells using the defined metric field. In this way we avoid the geometrical complexity of the clipping procedure and we show that we evaluate the energy and its gradient correctly. We show that the reconstruction of the computational domain is consistent with the Lloyds’ algorithm that is used to compute the L(sub p)-CVT.

Centroidal

DNS of Flow in a Low-Pressure Turbine Cascade Using a Discontinuous-Galerkin Spectral-Element Method

A new computational capability under development for accurate and efficient high-fidelity direct numerical simulation (DNS) and large eddy simulation (LES) of turbomachinery is described. This capability is based on an entropy-stable Discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy and is implemented in a computationally efficient manner on a modern high performance computer architecture. A validation study using this method to perform DNS of flow in a low-pressure turbine airfoil cascade are presented. Preliminary results indicate that the method captures the main features of the flow. Discrepancies between the predicted results and the experiments are likely due to the effects of freestream turbulence not being included in the simulation and will be addressed in the final paper.

Low-Pressure Turbine Cascade

Simulations of SSLV Ascent and Debris Transport

A viewgraph presentation on Computational Fluid Dynamic (CFD) Simulation of Space Shuttle Launch Vehicle (SSLV) ascent and debris transport analysis is shown. The topics include: 1) CFD simulations of the Space Shuttle Launch Vehicle ascent; 2) Debris transport analysis; 3) Debris aerodynamic modeling; and 4) Other applications.

Rogers, Stuart

Linearity-Preserving Limiters on Irregular Grids

This paper examines the behavior of flux and slope limiters on non-uniform grids in multiple dimensions. We note that on non-uniform grids the scalar formulation in standard use today sacrifices k-exactness, even for linear solutions, impacting both accuracy and convergence. We rewrite some well-known limiters in a n way to highlight their underlying symmetry, and use this to examine both traditional and novel limiter formulations. A consistent method of handling stretched meshes is developed, as is a new directional formulation in multiple dimensions for irregular grids. Results are presented demonstrating improved accuracy and convergence using a combination of model problems and complex three-dimensional examples.

Berger, Marsha

CFD Support for STS-107 Ascent Investigation

The research described in this viewgraph presentation investigates the ascent of STS-107 and foam-debris impact, and contributes to understanding of the STS-107 accident using CFD tools. The goals of the research are to: 1) Quantify loads on foam bipod ramp during ascent; 2) Provide steady-state flow-fields to debris-transport simulations; 3) Simulate flight of foam debris using unsteady six-degree-of-freedom calculations; 4) Provide estimates of foam mass, velocity, and impact angle which correlate with video and film evidence.

Rogers, Stuart E.

High Energy Boundary Conditions for a Cartesian Mesh Euler Solver

Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.

Pandya, Shishir

Parametric Powered-Lift Navier-Stokes Computations

The goal of this work is to enable the computation of large numbers of unsteady high-fidelity flow simulations for a YAV-8B Harrier aircraft in ground effect by improving the solution process and taking advantage of NASA parallel supercomputers. The YAV-8B Harrier aircraft can take off and land vertically, or utilize short runways by directing its four exhaust nozzles toward the ground. Transition to forward flight is achieved by rotating these nozzles into a horizontal position.

Chaderjian, Neal M.

Parametric Time-Dependent Navier-Stokes Computations for a YAV-8B Harrier in Ground Effect

The Harrier Jump Jet has the distinction of being the only powered-lift aircraft in the free world to achieve operational status and to have flown in combat. This V/STOL aircraft can take-off and land vertically or utilize very short runways by directing its four exhaust nozzles towards the ground. Transition to forward flight is achieved by rotating these nozzles into a horizontal position. Powered-lift vehicles have certain advantages over conventional strike fighters. Their V/STOL capabilities allow for safer carrier operations, smaller carrier size, and quick reaction time for troop support. Moreover, they are not dependent on vulnerable land-based runways. The AV-8A Harrier first entered service in the British Royal Air Force (RAF) during 1969, and the U.S. Marine Corps (USMC) in 1971. The AV-8B was a redesign to achieve improved payload capacity, range, and accuracy. This modified design first entered service with the USMC and RAF in 1985. The success and unique capabilities of the Harrier has prompted the design of a powered-lift version of the Joint Strike Fighter (JSF). The flowfield for the Harrier near the ground during low-speed or hover flight operations is very complex and time-dependent. A sketch of this flowfield is shown. Warm air from the fan is exhausted from the front nozzles, while a hot air/fuel mixture from the engine is exhausted from the rear nozzles. These jets strike the ground and move out radially forming a ground jet-flow. The ambient freestream, due to low-speed forward flight or - headwind during hover, opposes the jet-flow. This interaction causes the flow to separate and form a ground vortex. The multiple jets also interact with each other near the ground and form an upwash or jet fountain, which strikes the underside of the fuselage. If the aircraft is sufficiently close to the ground, the inlet can ingest ground debris and hot gases from the fountain and ground vortex. This Hot Gas Ingestion (HGI) can cause a sudden loss of thrust (powered lift), and the vehicle may crash. The high-speed jet flow along the ground can also entrain the ambient flow, resulting in a low pressure region underneath the vehicle. The accompanied loss of lift is referred to as the suckdown effect. The ground vortex may also be highly unsteady, dramatically changing its size and position with time at low frequencies, e.g., 1 Hz.

Chaderjian, Neal M.

Unsteady Computations of a Jet in a Crossflow with Ground Effect

A numerical study of a jet in crossflow with ground effect is conducted using OVERFLOW with dual time-stepping and low Mach number preconditioning. The results of the numerical study are compared to an experiment to show that the numerical methods are capable of capturing the dominant features of the flow field as well as the unsteadiness associated with the ground vortex.

Pandya, Shishir

A PTC Optimization and Control Surface Interference Study HSR Airframe Technical Review

This report considers the effect of canard and horizontal tail vertical position on the aerodynamic characteristics of the PTC configuration without nacelles and diverters. This analysis is followed by three optimization studies using canard and tail incidence as design variables in the first problem followed by an optimization run with canard and tail incidence and wing camber design variables and finally an optimization run with canard incidence and wing camber. The first problem was run at fixed lift while the other two problems were run at fixed angle of attack. The final investigation reported here will show data from a component buildup study using the PTC configuration. This final study will show the aerodynamic interference between the canard, wing and horizontal tail.

Hicks, Raymond

CFD Data Generation Process for Nonlinear Loads

This paper discusses the development of a process to generate a CFD database for the non-linear loads process capability for critical loads evaluation at Boeing Long Beach. The CFD simulations were performed for wing/body configurations at high angles of attack and Reynolds numbers with transonic and elastic deflection effects. Convergence criteria had to be tailored for loads applications rather than the usual drag performance. The time-accurate approach was subsequently adopted in order to improve convergence and model possible unsteadiness in the flowfield. In addition, uncertainty issues relating to the turbulence model and grid resolution in areas of high vortical flows were addressed and investigated for one of the cases.

Arslan, Alan