ALBERT

Description: The effect of turbulence models in the off-body grids on the accuracy of solutions for rotor flows in hover has been investigated. Results from the Reynolds-Averaged Navier-Stokes and Laminar Off-Body models are compared. Advection of turbulent eddy viscosity has been studied to find the mechanism leading to inaccurate solutions. A coaxial rotor result is also included.

Description: Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Description: An improved solution adaption capability has been implemented in the OVERFLOW overset grid CFD code. Building on the Cartesian off-body approach inherent in OVERFLOW and the original adaptive refinement method developed by Meakin, the new scheme provides for automated creation of multiple levels of finer Cartesian grids. Refinement can be based on the undivided second-difference of the flow solution variables, or on a specific flow quantity such as vorticity. Coupled with load-balancing and an inmemory solution interpolation procedure, the adaption process provides very good performance for time-accurate simulations on parallel compute platforms. A method of using refined, thin body-fitted grids combined with adaption in the off-body grids is presented, which maximizes the part of the domain subject to adaption. Two- and three-dimensional examples are used to illustrate the effectiveness and performance of the adaption scheme.

Description: An analysis is presented for obtaining two-dimensional temperature profiles and heat transfer in a radiation-absorbing gray gas of uniform absorptivity under the combined influence of thermal radiation, conduction, and gas flow. The gas is enclosed in a channel of infinite width and finite length formed by two semi-infinite parallel flat plates. These plates are black emitting surfaces, and the ends of the channel are formed by porous black surfaces through which the gas can flow into or out of the channel. These porous black end surfaces are used to simulate the radiation environment external to the channel. First, results are obtained for heat transfer between the plates in the absence of both conduction and flow. These results are found to be in good agreement with those obtained for the same conditions by previous workers. Results are then presented for heat transfer between the plates for the case of a radiating and conducting, but stagnant, gas separating the plates. The effects of the interactions between radiation and conduction are discussed. It was found that the heat transfer for combined radiation and conduction in an absorbing gas is slightly greater than the sum for each process taken separately. Finally, results are given for heat transfer from the plates to a flowing, radiating gas in the absence of conduction. The two plates are at the same temperature, and the gas enters the channel with uniform velocity and temperature. The results obtained for this case indicate that the heat transferred to the flowing gas from the constant temperature surfaces goes through a maximum as the absorptivity of the gas increases. This is in qualitative agreement with earlier results obtained by other investigators. All the results are presented in terms of dimensionless parameters, for the sake of generality, and the derivation of the dimensionless parameters, which are indicative of the effects of conduction and flow is presented.

Description: A thermal protection system (TPS) comprising a mixture of silicon carbide and SiOx that has been converted from Si that is present in a collection of diatom frustules and at least one diatom has quasi-periodic pore-to-pore separation distance d(p-p) in a selected range. Where a heat shield comprising the converted SiC/SiOx frustules receives radiation, associated with atmospheric (re)entry, a portion of this radiation is reflected so that radiation loading of the heat shield is reduced.

Description: Preconditioning methods have become the method of choice for the solution of flowfields involving the simultaneous presence of low Mach and transonic regions. It is well known that these methods are important for insuring accurate numerical discretization as well as convergence efficiency over various operating conditions such as low Mach number, low Reynolds number and high Strouhal numbers. For unsteady problems, the preconditioning is introduced within a dual-time framework wherein the physical time-derivatives are used to march the unsteady equations and the preconditioned time-derivatives are used for purposes of numerical discretization and iterative solution. In this paper, we describe the implementation of the preconditioned dual-time methodology in the OVERFLOW code. To demonstrate the performance of the method, we employ both simple and practical unsteady flowfields, including vortex propagation in a low Mach number flow, flowfield of an impulsively started plate (Stokes' first problem) arid a cylindrical jet in a low Mach number crossflow with ground effect. All the results demonstrate that the preconditioning algorithm is responsible for improvements to both numerical accuracy and convergence efficiency and, thereby, enables low Mach number unsteady computations to be performed at a fraction of the cost of traditional time-marching methods.

Description: Many thermal protection system materials used for spacecraft heatshields have anisotropic thermal properties, causing them to display significantly different thermal characteristics in different directions, when subjected to a heating environment during flight or arcjet tests. This paper investigates the effects of sidewall heating coupled with anisotropic thermal properties of thermal protection materials in the arcjet environment. Phenolic Impregnated Carbon Ablator (PICA) and LI-2200 materials (the insulation material of Shuttle tiles) were used for this study. First, conduction-based thermal response simulations were carried out, using the Marc.Mentat finite element solver, to study the effects of sidewall heating on PICA arcjet coupons. The simulation showed that sidewall heating plays a significant role in thermal response of these models. Arcjet tests at the Aerodynamic Heating Facility (AHF) at NASA Ames Research Center were performed later on instrumented coupons to obtain temperature history at sidewall and various radial locations. The details of instrumentation and experimental technique are the prime focus of this paper. The results obtained from testing confirmed that sidewall heating plays a significant role in thermal response of these models. The test results were later used to verify the two-dimensional ablation, thermal response, and sizing program, TITAN. The test data and model predictions were found to be in excellent agreement

Description: An investigation to determine the steady-state and surge characteristics of the J57-P-1 two-spool turbojet engine with various inlet air-flow distortions was conducted in the altitude wind tunnel at the NACA Lewis laboratory. Along with a uniform inlet total-pressure distribution, one circumferential and three radial pressure distortions were investigated. Data were obtained over a complete range of compressor speeds both with and without intercompressor air bleed at a flight Mach number of 0.8 and at altitudes of 35,000 and 50,000 feet. Total-pressure distortions of the magnitudes investigated had very little effect on the steady-state operating line for either the outer or inner compressor. The small radial distortions investigated also had engine over that obtained with the uniform inlet pressure distribution. The circumferential distortion, however, raised the minimum speed at which the engine could operate without encountering surge when the intercompressor bleeds were closed. This increase in minimum speed resulted in a substantial reduction in the operable speed range accompanied by a reduction in the altitude operating limit.

Description: A process for preparing a silicon based substrate with a protective coating having improved thermal resistance at temperature up to at least 1500.degree. C., and the resulting article.

Description: This research paper investigates current experimental designs and regression models for calibrating internal wind tunnel force balances of non-monolithic design. Such calibration methods are necessary for this class of balance because it has an electrical response that is dependent upon the sign of the applied forces and moments. This dependency gives rise to discontinuities in the response surfaces that are not easily modeled using traditional response surface methodologies. An analysis of current recommended calibration models is shown to lead to correlated response model terms. Alternative modeling methods are explored which feature orthogonal or near-orthogonal terms.