Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗
NASA NTRS1986

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗
NASA NTRS1987

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗
NASA NTRS2014

Conceptual Design Optimization of an Augmented Stability Aircraft Incorporating Dynamic Response Performance Constraints

This research focused on incorporating stability and control into a multidisciplinary de- sign optimization on a Boeing 737-class advanced concept called the D8.2b. A new method of evaluating the aircraft handling performance using quantitative evaluation of the sys- tem to disturbances, including perturbations, continuous turbulence, and discrete gusts, is presented. A multidisciplinary design optimization was performed using the D8.2b transport air- craft concept. The con guration was optimized for minimum fuel burn using a design range of 3,000 nautical miles. Optimization cases were run using xed tail volume coecients, static trim constraints, and static trim and dynamic response constraints. A Cessna 182T model was used to test the various dynamic analysis components, ensuring the analysis was behaving as expected. Results of the optimizations show that including stability and con- trol in the design process drastically alters the optimal design, indicating that stability and control should be included in conceptual design to avoid system level penalties later in the design process.

Welstead, Jason↗
NASA NTRS1982

ISEE-1 and -2 observations of laminar bow shocks - Velocity and thickness

The ISEE-1 and -2 spacecraft have provided the opportunity for studying the thickness and velocity of the bow shock on a regular basis. There are two models for the thickness of quasi-perpendicular laminar shocks. Galeev (1976) has calculated the thickness of the nearly perpendicular laminar bow shock by means of a weak ion sound turbulence (WIST) approach. Morse and Greenstadt (1976), on the other hand, assume that the shock thickness adjust itself to stay close to the drift velocity for marginal stability of current driven electrostatic waves. The present investigation is concerned with an evaluation of the competing theories, taking into account ISEE-1 and -2 observations of laminar bow shocks. It is found that the formula used by Morse and Greenstadt does not predict the observed thicknesses as well as the formula given by Galeev.

Russell, C. T.↗
NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗
NASA NTRS1992

Boundary layer structure over the ocean observed by LEANDRE 1 during a tramontane event

A new airborne backscatter lidar, has been developed by CNRS (Service d'Aeronomie, (SA) Laboratoire de Meteorologie Dynamique (LMD) and the Institut des Sciences de 1'Univers) in the frame of the LEANDRE research program. It has been qualified on board the ARAT in autumn 1989 and spring 1990 and was involved in its first cooperative campaign during PYREX in October and November 1990. During this campaign, lidar observations of the perturbations induced on tropospheric flow and boundary layer structure were performed, and results are presented. A large number of experiments were performed, for synoptic situation description (meteorological radiosoundings, constant level balloons) and local flow analysis (aircrafts, radars, sodars). For the first time in such an experiment, a lidar has been flown on a research aircraft to perform altitude resolved observations of these perturbations, and we will present here results obtained for deflected flow structure. In the presence of a synoptic northerly flow, part of it is deflected to the east by the Pyrenees, and accelerated over the Mediterranean by the mountain surroundings. In this case, a low level wind is generated (the Tramontane) bringing cold and dry air over the Mediterranean Sea. As the sea is still at a warm temperature in November (around 17 degrees), an Internal Marine Boundary Layer rapidly grows over the first tens of kilometers and stabilizes at about 1 km depth, corresponding to an altitude just below the Lifting Condensation Level. The whole Marine Atmospheric Boundary Layer (MABL) is characterized by highly turbulent motions bringing large humid particles from the surface up to its top. The lidar signal due to scattering by these particles is then representative of the turbulent kinetic energy in this layer.

Flamant, C.↗
NASA NTRS1993

Effect of micron-sized roughness on transition in swept-wing flows

Boundary-layer transition-to-turbulence studies are conducted in the Arizona State University Unsteady Wind Tunnel on a 45-degree swept airfoil. The pressure gradient is designed so that the initial stability characteristics are purely crossflow-dominated. Flow visualization and hot-wire measurements show that the development of the crossflow vortices is influenced by roughness near the attachment-line. Comparisons of transition location are made between a painted surface, a machine-polished surface, and a hand-polished surface. Then, isolated 6 micron roughness elements are placed near the attachment line on the airfoil surface under conditions of the final polish (0.25 micron rms). These elements amplify a centered stationary crossflow vortex and its neighbors, resulting in localized early transition. The diameter, height, and location of these roughness elements are varied in a systematic manner. Spanwise hot-wire measurements are taken behind the roughness element to document the enhanced vortices. These scans are made at several different chord locations to examine vortex growth.

Radeztsky, Ronald H., Jr.↗
NASA NTRS2023

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally,Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗
NASA NTRS2023

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally, Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗
NASA NTRS1981

Forecasting of Clear Air Turbulence using a Diagnostic Richardson Number Tendency formulation

The results of several case studies of Clear Air Turbulence (CAT) using the Diagnostic Richardson Number Tendency (DRT) formulation are highlighted. The performance of this technique in resolving regions of documented CAT encounters is encouraging. Its operational adaptability seems particularly attractive in that the input data can be supplied by the currently operational Rawinsonde system. A CAT index is calculated deterministically, sensing synoptic-scale changes in static stability and vertical wind shear conducive for supporting meso-scale CAT layers. This index reveals volumes of the troposphere which act as 'source regions' of Kelvin-Helmholtz instabilities. The results suggest that these regions are particularly efficient with respect to the synoptic/meso-scale energetical coupling necessary for supporting significantly turbulent layers. The output highlights specific regions of the atmosphere which can be interpreted operationally in terms of CAT-encounter probabilities.

Keller, J. L.↗
NASA NTRS1991

Hi-alpha forebody design. Part 2: Determination of body shapes for positive directional stability

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low speed high angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows of the F-5A forebody and Erickson forebody were recomputed with better and refined grids. The results were obtained using a modified version of cfl3d to solve either the Euler equations or the Reynolds equations employing a form of the Baldwin-Lomax turbulence model. Based on those results, we conclude that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a systematic study of forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. A novel way of presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, R.↗
NASA NTRS2001

A Multiyear Dataset of SSM/I-Derived Global Ocean Surface Turbulent Fluxes

The surface turbulent fluxes of momentum, latent heat, and sensible heat over global oceans are essential to weather, climate and ocean problems. Evaporation is a key component of the hydrological cycle and the surface heat budget, while the wind stress is the major forcing for driving the oceanic circulation. The global air-sea fluxes of momentum, latent and sensible heat, radiation, and freshwater (precipitation-evaporation) are the forcing for driving oceanic circulation and, hence, are essential for understanding the general circulation of global oceans. The global air-sea fluxes are required for driving ocean models and validating coupled ocean-atmosphere global models. We have produced a 7.5-year (July 1987-December 1994) dataset of daily surface turbulent fluxes over the global oceans from the Special Sensor microwave/Imager (SSM/I) data. Daily turbulent fluxes were derived from daily data of SSM/I surface winds and specific humidity, National Centers for Environmental Prediction (NCEP) sea surface temperatures, and European Centre for Medium-Range Weather Forecasts (ECMWF) air-sea temperature differences, using a stability-dependent bulk scheme. The retrieved instantaneous surface air humidity (with a 25-km resolution) validated well with that of the collocated radiosonde observations over the global oceans. Furthermore, the retrieved daily wind stresses and latent heat fluxes were found to agree well with that of the in situ measurements (IMET buoy, RV Moana Wave, and RV Wecoma) in the western Pacific warm pool during the TOGA COARE intensive observing period (November 1992-February 1993). The global distributions of 1988-94 seasonal-mean turbulent fluxes will be presented. In addition, the global distributions of 1990-93 annual-means turbulent fluxes and input variables will be compared with those of UWM/COADS covering the same period. The latter is based on the COADS (comprehensive ocean-atmosphere data set) and is recognized to be one of the best climatological analyses of fluxes derived from ship observations.

Chou, Shu-Hsien↗
NASA NTRS2004

COMSAC: Computational Methods for Stability and Control

The unprecedented advances being made in computational fluid dynamic (CFD) technology have demonstrated the powerful capabilities of codes in applications to civil and military aircraft. Used in conjunction with wind-tunnel and flight investigations, many codes are now routinely used by designers in diverse applications such as aerodynamic performance predictions and propulsion integration. Typically, these codes are most reliable for attached, steady, and predominantly turbulent flows. As a result of increasing reliability and confidence in CFD, wind-tunnel testing for some new configurations has been substantially reduced in key areas, such as wing trade studies for mission performance guarantees. Interest is now growing in the application of computational methods to other critical design challenges. One of the most important disciplinary elements for civil and military aircraft is prediction of stability and control characteristics. CFD offers the potential for significantly increasing the basic understanding, prediction, and control of flow phenomena associated with requirements for satisfactory aircraft handling characteristics.

Fremaux, C. Michael↗
NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗
NASA NTRS1982

Exhaust gas measurements in a propane fueled swirl stabilized combustor

Exhaust gas temperature, velocity, and composition are measured and combustor efficiencies are calculated in a lean premixed swirl stabilized laboratory combustor. The radial profiles of the data between the co- and the counter swirl cases show significant differences. Co-swirl cases show evidence of poor turbulent mixing across the combustor in comparison to the counter-swirl cases. NO sub x levels are low in the combustor but substantial amounts of CO are present. Combustion efficiencies are low and surprisingly constant with varying outer swirl in contradiction to previous results under a slightly different inner swirl condition. This difference in the efficiency trends is expected to be a result of the high sensitivity of the combustor to changes in the inner swirl. Combustor operation is found to be the same for propane and methane fuels. A mechanism is proposed to explain the combustor operation and a few important characteristics determining combustor efficiency are identified.

Aanad, M. S.↗