Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Edge-Based Viscous Method for Mixed-Element Node-Centered Finite-Volume Solvers

A novel, efficient, edge-based viscous (EBV) discretization method has been recently developed, implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver, and applied to viscous-kernel computations that include evaluations of meanflow viscous fluxes, turbulence-model and chemistry-model diffusion terms, and the corresponding Jacobian contributions. Initially, the EBV method had been implemented for tetrahedral grids and demonstrated multifold acceleration of all viscous-kernel computations. This paper presents an extension of the EBV method for mixed-element grids. In addition to the primal edges of a given mixed-element grid, virtual edges are introduced to connect cell nodes that are not connected by a primal edge. The EBV method uses an efficient loop over all (primal and virtual) edges and features a compact discretization stencil based on the nearest neighbors. This study verifies the EBV method and assesses its efficiency on mixed-element grids by comparing the EBV solution accuracy and iterative convergence with those of well-established solutions obtained using a cell-based viscous (CBV) discretization method. The EBV solver’s memory footprint is optimized and often smaller than the memory footprint of the CBV solver. A multifold speedup is demonstrated for all viscous-kernel computations resulting in significant reduction of the time to solutions for several benchmark mixed-element-grid computations, including simulations of a flow around NASA’s juncture-flow model and a hypersonic, chemically reacting flow around a blunt body.

Edge-based viscous method↗
NASA NTRS1999

NCC: A Multidisciplinary Design/Analysis Tool for Combustion Systems

A multi-disciplinary design/analysis tool for combustion systems is critical for optimizing the low-emission, high-performance combustor design process. Based on discussions between NASA Lewis Research Center and the jet engine companies, an industry-government team was formed in early 1995 to develop the National Combustion Code (NCC), which is an integrated system of computer codes for the design and analysis of combustion systems. NCC has advanced features that address the need to meet designer's requirements such as "assured accuracy", "fast turnaround", and "acceptable cost". The NCC development team is comprised of Allison Engine Company (Allison), CFD Research Corporation (CFDRC), GE Aircraft Engines (GEAE), NASA Lewis Research Center (LeRC), and Pratt & Whitney (P&W). This development team operates under the guidance of the NCC steering committee. The "unstructured mesh" capability and "parallel computing" are fundamental features of NCC from its inception. The NCC system is composed of a set of "elements" which includes grid generator, main flow solver, turbulence module, turbulence and chemistry interaction module, chemistry module, spray module, radiation heat transfer module, data visualization module, and a post-processor for evaluating engine performance parameters. Each element may have contributions from several team members. Such a multi-source multi-element system needs to be integrated in a way that facilitates inter-module data communication, flexibility in module selection, and ease of integration.

Liu, Nan-Suey↗
NASA NTRS2012

A Full-Envelope Air Data Calibration and Three-Dimensional Wind Estimation Method Using Global Output-Error Optimization and Flight-Test Techniques

A novel, efficient air data calibration method is proposed for aircraft with limited envelopes. This method uses output-error optimization on three-dimensional inertial velocities to estimate calibration and wind parameters. Calibration parameters are based on assumed calibration models for static pressure, angle of attack, and flank angle. Estimated wind parameters are the north, east, and down components. The only assumptions needed for this method are that the inertial velocities and Euler angles are accurate, the calibration models are correct, and that the steady-state component of wind is constant throughout the maneuver. A two-minute maneuver was designed to excite the aircraft over the range of air data calibration parameters and de-correlate the angle-of-attack bias from the vertical component of wind. Simulation of the X-48B (The Boeing Company, Chicago, Illinois) aircraft was used to validate the method, ultimately using data derived from wind-tunnel testing to simulate the un-calibrated air data measurements. Results from the simulation were accurate and robust to turbulence levels comparable to those observed in flight. Future experiments are planned to evaluate the proposed air data calibration in a flight environment.

Taylor, Brian R.↗
NASA NTRS2022

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗
NASA NTRS2016

Development of a High-Order Space-Time Matrix-Free Adjoint Solver

The growth in computational power and algorithm development in the past few decades has granted the science and engineering community the ability to simulate flows over complex geometries, thus making Computational Fluid Dynamics (CFD) tools indispensable in analysis and design. Currently, one of the pacing items limiting the utility of CFD for general problems is the prediction of unsteady turbulent ows.1{3 Reynolds-averaged Navier-Stokes (RANS) methods, which predict a time-invariant mean flowfield, struggle to provide consistent predictions when encountering even mild separation, such as the side-of-body separation at a wing-body junction. NASA's Transformative Tools and Technologies project is developing both numerical methods and physical modeling approaches to improve the prediction of separated flows. A major focus of this e ort is efficient methods for resolving the unsteady fluctuations occurring in these flows to provide valuable engineering data of the time-accurate flow field for buffet analysis, vortex shedding, etc. This approach encompasses unsteady RANS (URANS), large-eddy simulations (LES), and hybrid LES-RANS approaches such as Detached Eddy Simulations (DES). These unsteady approaches are inherently more expensive than traditional engineering RANS approaches, hence every e ort to mitigate this cost must be leveraged. Arguably, the most cost-effective approach to improve the efficiency of unsteady methods is the optimal placement of the spatial and temporal degrees of freedom (DOF) using solution-adaptive methods.

Adjoint↗
NASA NTRS1991

A unified viscous theory of lift and drag of 2-D thin airfoils and 3-D thin wings

A unified viscous theory of 2-D thin airfoils and 3-D thin wings is developed with numerical examples. The viscous theory of the load distribution is unique and tends to the classical inviscid result with Kutta condition in the high Reynolds number limit. A new theory of 2-D section induced drag is introduced with specific applications to three cases of interest: (1) constant angle of attack; (2) parabolic camber; and (3) a flapped airfoil. The first case is also extended to a profiled leading edge foil. The well-known drag due to absence of leading edge suction is derived from the viscous theory. It is independent of Reynolds number for zero thickness and varies inversely with the square root of the Reynolds number based on the leading edge radius for profiled sections. The role of turbulence in the section induced drag problem is discussed. A theory of minimum section induced drag is derived and applied. For low Reynolds number the minimum drag load tends to the constant angle of attack solution and for high Reynolds number to an approximation of the parabolic camber solution. The parabolic camber section induced drag is about 4 percent greater than the ideal minimum at high Reynolds number. Two new concepts, the viscous induced drag angle and the viscous induced separation potential are introduced. The separation potential is calculated for three 2-D cases and for a 3-D rectangular wing. The potential is calculated with input from a standard doublet lattice wing code without recourse to any boundary layer calculations. Separation is indicated in regions where it is observed experimentally. The classical induced drag is recovered in the 3-D high Reynolds number limit with an additional contribution that is Reynold number dependent. The 3-D viscous theory of minimum induced drag yields an equation for the optimal spanwise and chordwise load distribution. The design of optimal wing tip planforms and camber distributions is possible with the viscous 3-D wing theory.

Yates, John E.↗
NASA NTRSDate not supplied

Uranus Probe Entry and Descent Mission Concept

Introduction: Uranus was identified as the third highest priority flagship mission in the 2012-2022 Planetary Science Decadal Survey. This latest concept study was requested by the Decadal Survey panel to determine NASA’s planetary science priorities from 2022-2032. This study focused on the probe’s entry and descent aspects and associated trades for viable trajectory options. Uranus Mission and Descent Probe: The proposed Uranus Orbiter and Probe (UOP) Flagship mission will investigate Uranus and its surrounding moons using an orbiting spacecraft with a Uranus descent probe. Unlike previous studies the probe release will occur after orbit insertion allowing sufficient separation of critical events during the orbit insertion burn. The probe will be released at an altitude that allows one hour of in situ atmospheric readings that will be relayed to the orbiter. Afterwards the orbiter will transition to the moon tour phase of the mission. The configuration chosen for the entry aeroshell was a 45° sphere-cone. This shape has been used in the past in the Pioneer Venus Galileo missions. However, the nose radius considered in the present study differed from the values used in either of the previous configurations, primarily to reduce the heat flux at the stagnation point. A two-step approach was used in the development of flight trajectories for the chosen configuration. In the first step, the trajectory code POST2 was used to screen the thousands of entry states provided by interplanetary trajectory simulations, which were terminated at an altitude of 2000 km from the reference surface (1 bar) of Uranus. The screening criteria were: (i) optimization of the communication geometry between the entry probe and orbiter to ensure at least 1 hour of science measurements, (ii) peak stagnation point pressures to be less than six bar, (iii) peak heat fluxes to be less than 5 kW/cm2; the latter two constraints being the limits of ground-test capabilities of the arc jets at NASA Ames Research Center. The entry team investigated two trajectories that met the criteria above, a shallow entry (high heat load ~44 kJ/cm2) and a steeper entry (high heat rate ~1950 W/cm2). In the second step, the two bounding candidate entry states from the POST2 screening process were used in developing flight trajectories using TRAJ coupled with FIAT (a materials thermal response and sizing code) and a margins policy to determine a margined uniform thickness (hence mass) of the forward heatshield material based on the aerothermal environments at the stagnation point. Since the combination of TRAJ and FIAT size the TPS based on stagnation point environments only, flow field computations using DPLR were necessary to determine turbulent aerothermal environments on the conical flank, and the augmentation of these environments due to surface roughness. The environments at select locations on the forward heatshield were then used to size the thermal protection material, with the largest thickness value then used to estimate the mass. The newly developed woven thermal protection material –HEEET (Heatshield for Extreme Entry Environments Technology) –was considered for the forward heatshield and PICA (Phenolic-Impregnated Carbon Ablator) was considered for the backshell. These NASA-developed materials are at TRL 6 and TRL 9, respectively. Furthermore, two options were considered for the HEEET material: (i) a dual-layer option with a denser recession layer on top and an insulative layer underneath it, and (ii) a single layer option consisting of the insulative layer alone. Results: It is clear that probe entry states are feasible and the selected TPS options are able to perform in the predicted aerothermal environments thus enabling the mission to meet of the descent probe portion of this flagship mission

Entry Descent and Landing↗
NASA NTRS2004

Usability Evaluation of a Flight-Deck Airflow Hazard Visualization System

Many aircraft accidents each year are caused by encounters with unseen airflow hazards near the ground, such as vortices, downdrafts, low level wind shear, microbursts, or turbulence from surrounding vegetation or structures near the landing site. These hazards can be dangerous even to airliners; there have been hundreds of fatalities in the United States in the last two decades attributable to airliner encounters with microbursts and low level wind shear alone. However, helicopters are especially vulnerable to airflow hazards because they often have to operate in confined spaces and under operationally stressful conditions (such as emergency search and rescue, military or shipboard operations). Providing helicopter pilots with an augmented-reality display visualizing local airflow hazards may be of significant benefit. However, the form such a visualization might take, and whether it does indeed provide a benefit, had not been studied before our experiment. We recruited experienced military and civilian helicopter pilots for a preliminary usability study to evaluate a prototype augmented-reality visualization system. The study had two goals: first, to assess the efficacy of presenting airflow data in flight; and second, to obtain expert feedback on sample presentations of hazard indicators to refine our design choices. The study addressed the optimal way to provide critical safety information to the pilot, what level of detail to provide, whether to display specific aerodynamic causes or potential effects only, and how to safely and effectively shift the locus of attention during a high-workload task. Three-dimensional visual cues, with varying shape, color, transparency, texture, depth cueing, and use of motion, depicting regions of hazardous airflow, were developed and presented to the pilots. The study results indicated that such a visualization system could be of significant value in improving safety during critical takeoff and landing operations, and also gave clear indications of the best design choices in producing the hazard visual cues.

Aragon, Cecilia R.↗
NASA NTRS2010

HYTHIRM Radiance Modeling and Image Analyses in Support of STS-119, STS-125 and STS-128 Space Shuttle Hypersonic Re-entries

We provide the first geometrically accurate (i.e., 3-D) temperature maps of the entire windward surface of the Space Shuttle during hypersonic reentry. To accomplish this task we began with estimated surface temperatures derived from CFD models at integral high Mach numbers and used them, the Shuttle's surface properties and reasonable estimates of the sensor-to-target geometry to predict the emitted spectral radiance from the surface (in units of W sr-1 m-2 nm-1). These data were converted to sensor counts using properties of the sensor (e.g. aperture, spectral band, and various efficiencies), the expected background, and the atmosphere transmission to inform the optimal settings for the near-infrared and midwave IR cameras on the Cast Glance aircraft. Once these data were collected, calibrated, edited, registered and co-added we formed both 2-D maps of the scene in the above units and 3-D maps of the bottom surface in temperature that could be compared with not only the initial inputs but also thermocouple data from the Shuttle itself. The 3-D temperature mapping process was based on the initial radiance modeling process. Here temperatures were guessed for each node in a well-resolved 3-D framework, a radiance model was produced and compared to the processed imagery, and corrections to the temperature were estimated until the iterative process converged. This process did very well in characterizing the temperature structure of the large asymmetric boundary layer transition the covered much of the starboard bottom surface of STS-119 Discovery. Both internally estimated accuracies and differences with CFD models and thermocouple measurements are at most a few percent. The technique did less well characterizing the temperature structure of the turbulent wedge behind the trip due to limitations in understanding the true sensor resolution. (Note: Those less inclined to read the entire paper are encouraged to read an Executive Summary provided at the end.)

Gibson, David M.↗
NASA NTRSDate not supplied

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

The reliance on CFD simulations to develop, design, and optimize scramjet systems (or components) has become commonplace. This reliance inevitably hinges on the ability of the computational analyst to quantify the level of confidence in their computational results. Unfortunately, nearly all the measured data available for this assessment comes from antiquated experimental datasets, or from tests that focused on the extraction of scramjet system (or component) performance. The objective of a CFD validation experiment is to quantify the predictive accuracy of one or more of the CFD physics submodels, implying that other uncertainties related to replicating the facility flow environment (e.g., knowledge of boundary conditions) must be minimized to the extent possible. This inevitably places stringent requirements on the quality and quantity of measurements taken to accurately specify inflow, outflow, and surface conditions for the CFD simulations; in addition to the measurements taken for the validation of physics submodels. This places additional demands on the experimental process above and beyond those for test article performance assessment. A recent high speed code credibility workshop series sponsored by AFRL identified a gap in existing validation data for fundamental assessments of turbulent mixing and combustion CFD closure models at scramjet engine relevant conditions. To address this gap, engineers at AFRL have designed a coaxial jet flame configuration that will be tested at two facilities (Research Cell 19 at the Air Force Research Lab, and at Purdue University). The effort described here documents pretest simulations of this validation experiment with the goal of fleshing out the extent of the facility flowpath that must be included to adequately reproduce the facility test section flow environment. The findings indicate that the flow around the support structure for the fuel injection centerbody upstream of the facility nozzle generates disturbances that persist throughout the nozzle expansion process; corrupting the azimuthal symmetry that was desired in the fuel/air mixing region of the test section. Simulations without this support structure maintained a high degree of azimuthal symmetry up until the fuel injection plane. However, even in this scenario the azimuthal symmetry was not maintained once the centerbody boundary layer transitioned to a wake flow downstream of the fuel injection plane.

CFD↗
NASA NTRSDate not supplied

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

The reliance on CFD simulations to develop, design, and optimize scramjet systems (or components) has become commonplace. This reliance inevitably hinges on the ability of the computational analyst to quantify the level of confidence in their computational results. Unfortunately, nearly all the measured data available for this assessment comes from antiquated experimental datasets, or from tests that focused on the extraction of scramjet system (or component) performance. The objective of a CFD validation experiment is to quantify the predictive accuracy of one or more of the CFD physics submodels, implying that other uncertainties related to replicating the facility flow environment (e.g., knowledge of boundary conditions) must be minimized to the extent possible. This inevitably places stringent requirements on the quality and quantity of measurements taken to accurately specify inflow, outflow, and surface conditions for the CFD simulations; in addition to the measurements taken for the validation of physics submodels. This places additional demands on the experimental process above and beyond those for test article performance assessment. A recent high speed code credibility workshop series sponsored by AFRL identified a gap in existing validation data for fundamental assessments of turbulent mixing and combustion CFD closure models at scramjet engine relevant conditions. To address this gap, engineers at AFRL have designed a coaxial jet flame configuration that will be tested at two facilities (Research Cell 19 at the Air Force Research Lab, and at Purdue University). The effort described here documents pretest simulations of this validation experiment with the goal of fleshing out the extent of the facility flowpath that must be included to adequately reproduce the facility test section flow environment. The findings indicate that the flow around the support structure for the fuel injection centerbody upstream of the facility nozzle generates disturbances that persist throughout the nozzle expansion process; corrupting the azimuthal symmetry that was desired in the fuel/air mixing region of the test section. Simulations without this support structure maintained a high degree of azimuthal symmetry up until the fuel injection plane. However, even in this scenario the azimuthal symmetry was not maintained once the centerbody boundary layer transitioned to a wake flow downstream of the fuel injection plane.

CFD↗
NASA NTRS2010

Atmospheric Compensation for Uplink Arrays via Radiometry

Uplink arrays for communications applications are gaining increased visibility within the NASA and military community due to the enhanced flexibility and reliability they provide. When compared with the conventional large, single aperture antennas currently comprising the Deep Space Network (DSN), for example, smaller aperture antenna arrays have the benefits of providing fault tolerance (reduced single-point failure), reduced maintenance cost, and enhanced capabilities such as electronic beam-steering and multi-beam operation. However, signal combining of antenna array elements spaced many wavelengths apart becomes problematic due to the inherent instability of earth's turbulent atmosphere, particularly at the frequencies of interest to the DSN (i.e., Ka-band). Degradation in the power combining of the individual elements comprising the array arises due to uncorrelated phase errors introduced as the signals propagate through the troposphere. It is well known that the fundamental source of this error is due to the inhomogeneous distribution of water vapor in the atmosphere [1]. Several techniques have been proposed to circumvent this issue, including the use of phase calibration towers and a moon bounce to generate a feedback loop which would provide a means of intermittent calibration of the system phase errors (thermal drifts, atmosphere) [2,3]. However, these techniques require repositioning of the antenna elements to perform this operation which ultimately results in reduced system availability. And, though they are sufficient for compensating for slow varying phase drifts, they are insufficient to compensate for faster varying phase errors, such as those introduced by the atmosphere. In this paper, preliminary radiometry and interferometry measurements collected by the NASA Glenn Research Center are analyzed and indicate that the use of optimized water vapor radiometers as a feedback system in a communications platform could provide the necessary atmospheric compensation technique to enhance the beamforming of uplink arrays.

Nessel, James A.↗
NASA NTRS1993

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗
NASA NTRS2016

Nonlinear Transient Growth and Boundary Layer Transition

Parabolized stability equations (PSE) are used in a variational approach to study the optimal, non-modal disturbance growth in a Mach 3 at plate boundary layer and a Mach 6 circular cone boundary layer. As noted in previous works, the optimal initial disturbances correspond to steady counter-rotating streamwise vortices, which subsequently lead to the formation of streamwise-elongated structures, i.e., streaks, via a lift-up effect. The nonlinear evolution of the linearly optimal stationary perturbations is computed using the nonlinear plane-marching PSE for stationary perturbations. A fully implicit marching technique is used to facilitate the computation of nonlinear streaks with large amplitudes. To assess the effect of the finite-amplitude streaks on transition, the linear form of plane- marching PSE is used to investigate the instability of the boundary layer flow modified by spanwise periodic streaks. The onset of bypass transition is estimated by using an N- factor criterion based on the amplification of the streak instabilities. Results show that, for both flow configurations of interest, streaks of sufficiently large amplitude can lead to significantly earlier onset of transition than that in an unperturbed boundary layer without any streaks.

Laminar turbulent transition↗
NASA NTRS1994

Optimized design of a hypersonic nozzle

Conventional procedures for designing nozzles involve the design of an inviscid contour (using the method of characteristics) that is corrected with a displacement thickness calculated from boundary-layer theory. However, nozzles designed using this classical procedure have been shown to exhibit poor flow quality at Mach numbers characteristic of hypersonic applications. The nozzle to be designed will be a part of the NASA HYPULSE facility which is being used for hypervelocity flight research. Thus, the flow quality of the nozzle is a critical question that needs to be addressed. Design of nozzles for hypersonic applications requires a proper assessment of the effects of the thick boundary layer on the inviscid flowfield. Since the flow field is largely supersonic, the parabolized form of the Navier-Stokes (PNS) equations can be used. The requirement of a uniform flow at the exit plane of the nozzle can be used to define an objective function as part of an optimization procedure. The design procedure used in this study involves the coupling of a nonlinear (least-squares) optimization algorithm with an efficient, explicit PNS solver. The thick boundary layers growing on the walls of the nozzle limit the extent of the usable core region (region with uniform flow) for testing models (especially rectangular). In order to maximize the region of uniform flow, it was decided to have the exit plane of this nozzle to be (nearly) rectangular. Thus, an additional constraint on the nozzle shape resulted, namely the nozzle will have a shape transitioning from a circular one at the inlet to that of a rectangle at the exit. In order to provide for a smooth shape transition, the cross sectional contour of the nozzle is defined by a superellipse. The nozzle is taken to be a meter in length. The axial variations of the major and minor radii of the superellipse are governed by cubic splines. The design parameters are the coefficients of the splines associated with the local nozzle wall slopes. Extensive calculations have been made (with a three-dimensional Euler code) to understand the effects of various parameters such as location of the knot points of the spline function, different ways of characterizing the uniformity of the flow in the exit plane, as well as the effect of constraining the area of the nozzle to be invariant. Turbulent flow (measurements indicate that the flow at the nozzle inlet is turbulent) calculations are now being performed (with the inviscidly designed nozzle contours) to assess the flow quality.

Krishnamurthy, Ramesh↗
NASA NTRSDate not supplied

A Data-Driven Reduced Order Model of an Isolated Rotor

There are numerous conceptual design stage rotorcraft analysis tasks which demand a high-fidelity and low cost method for rotor load distribution predictions. Considering Urban Air Mobility (UAM) vehicles aim to operate in close proximity to buildings and with unique rotor configurations, there is a significant challenge in quickly and accurately modeling rotors operating in complex, turbulent flow fields. One potential path for deriving a high-fidelity, low cost rotor model is with data-driven surrogate modeling. In this study, an initial investigation is taken to apply a proper orthogonal decomposition (POD) based reduced order model (ROM) for the purpose of pressure distribution prediction. In this study, a POD ROM was derived to produce distributed pressure predictions on rotor blades subjected to topology change due to variation in twist and taper ratio. Rotor twist was varied between 0◦, 10◦, 20◦, and 30◦ while taper ratio was varied between 1.0, 0.9, 0.8, and 0.7. All rotors consisted of a single blade. The POD ROM was validated for three demonstration cases; a high thrust rotor in hover, a low thrust rotor in hover, and a rotor in forward flight with a flight speed of M = 0.1. Results showed highly accurate distributed load predictions could be achieved at minimal computational cost. Computational cost for hovering blade surface pressure modeling was reduced from 12 hours on 440 cores to 10−5 seconds on a single core. For blade in forward flight cost was reduced from 20 hours on 440 cores to 0.6 seconds on a single core. For cases of high thrust and low thrust rotors, POD ROM was used to undergo a design optimization of the rotor such that figure of merit was maximized. Total optimization time for each case was 1 minute.

Data-Driven↗
NASA NTRS2006

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗
NASA NTRS2005

The Impact of Model Uncertainty on Spatial Compensation in Structural Acoustic Control

Turbulent boundary layer (TBL) noise is considered a primary contribution to the interior noise present in commercial airliners. There are numerous investigations of interior noise control devoted to aircraft panels; however, practical realization is a potential challenge since physical boundary conditions are uncertain at best. In most prior studies, pinned or clamped boundary conditions were assumed; however, realistic panels likely display a range of boundary conditions between these two limits. Uncertainty in boundary conditions is a challenge for control system designers, both in terms of the compensator implemented and the location of transducers required to achieve the desired control. The impact of model uncertainties, specifically uncertain boundaries, on the selection of transducer locations for structural acoustic control is considered herein. The final goal of this work is the design of an aircraft panel structure that can reduce TBL noise transmission through the use of a completely adaptive, single-input, single-output control system. The feasibility of this goal is demonstrated through the creation of a detailed analytical solution, followed by the implementation of a test model in a transmission loss apparatus. Successfully realizing a control system robust to variations in boundary conditions can lead to the design and implementation of practical adaptive structures that could be used to control the transmission of sound to the interior of aircraft. Results from this research effort indicate it is possible to optimize the design of actuator and sensor location and aperture, minimizing the impact of boundary conditions on the desired structural acoustic control.

Clark, Robert L.↗